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Interactions of the Heavy Nuclei of the Cosmic Radiation af 


By P. H. Fowter, R. R. Hirer and C. J. WADDINGTON 
H. H. Wills Physical Laboratory, Bristol 


[Received August 10, 1956] 


ABSTRACT 


‘ We have observed in a systematic manner 317 nuclear interactions 
produced by the heavy nuclei of the primary cosmic radiation in nuclear 
emulsions. From a study of these interactions we have found the inter- 
action mean free paths and fragmentation probabilities. The interaction 
mean free paths are in good agreement with those reported previously. 
The fragmentation probabilities, particularly those for the production of 
light nuclei, are considerably smaller than those currently accepted. In 
addition, from a study of the angular distribution of the nuclei emitted 
from the interactions we have shown that the energies of emission in the 
C-system of the incident nucleus are greater than those previously accepted 

-and depend on the charge of the primary nucleus. 


§ 1. INTRODUCTION 


Tue heavy nuclei in the cosmic radiation, those with charges equal to or 
greater than that of lithium, are rapidly absorbed in the atmosphere by 
nuclear interactions. As a result, in order to study these nuclei it is 
necessary to place the detecting equipment as high in the atmosphere as 
possible. However, with present day techniques it is impossible to keep 
such equipment above the atmosphere for any length of time. It is there- 
fore necessary to make measurements with matter above the detectors, 
and to correct for the effects produced by this residual matter. These 
effects are threefold : 

(i) Nuclei with comparatively low primary energies may be brought to 
rest by ionization energy losses. The correction introduced in this manner 
may be calculated. 

(ii) Nuclei may be broken up into nucleons and singly charged particles 
in nuclear interactions above the detectors, and thus not be observed. 
(iii) Nuclei may arrive at the detectors with reduced charge due to 
nuclear interactions in which they were broken into one or more heavy 
fragments. , pest” 
In this paper we are concerned with effects (11) and. (iil) > Le. with deter- 
mining the interaction mean free paths and fragmentation probabilities. 
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In the primary radiation there appear to be nuclei with all charges 
between those of lithium and iron. It is not, as yet, possible to determine 
these quantities for individual elements. Indeed it has been customary 
to separate these nuclei into only three charge groups. However, due to 
the somewhat greater statistical weight of our results, and because we 
intend to increase this weight later, we have separated our data into four 
groups. These we have defined as : light nuclei, L, 3<Z<5; medium 
nuclei, M, 6<Z<9; heavy nuclei, H, 10<Z<19 ; and very heavy nuclei, 
VA, Z<20. 

We have determined the mean free paths and fragmentation proba- 
bilities of these charge groups in nuclear emulsion. An estimate of the 
corresponding values in air has then been deduced from these experimental 
values. In addition, from the angular distribution of the fragments 
produced in these interactions we have been able to estimate the r.m.s. 
transverse momentum for emission of «-particles in the C-system of the 
incident nucleus. 


§ 2. EXPERIMENTAL PROCEDURE 


We have used a stack of 80, 20 x 15 em, 600 » Ilford G5 stripped photo- 
graphic nuclear emulsions. This stack was flown on the 14th of October 
1954 over Northern Italy, at a mean altitude of 105 000 feet. The cut- 
off energy at the top of the atmosphere for multiply charged nuclei has 
been determined previously as 1-55--0-06 Bev per nucleon (Fowler and 
Waddington 1956). We have therefore assumed that all the nuclei 
arriving at the emulsions had energies greater than | Bev per nucleon. 

In order to detect the tracks of the heavy nuclei a line 1 em below the 
top edge of the emulsions was scanned for tracks longer than a certain 
minimum length per emulsion. This length ranged between 2 and 6 mm. 
The track of every nucleus found in this manner was followed through 
the stack until the nucleus either interacted or left the stack. Ifa heavy 
fragment emerged from an interaction, this fragment was also followed. 

The charges of the nuclei were measured in two ways. For those nuclei 
with charges between those of lithium and oxygen inclusive a measure- 
ment was made of the 5-ray density by counting the number of 5-rays with 
four or more grains. For those nuclei with charges greater than or equal 
to that of oxygen the 6-ray density was measured by a photoelectric 
method. It is estimated that these charge determinations each had a 
probable error of about a charge. Fuller details of these measurements 
will be given in a later paper. 

Each interaction was classified in terms of the charge of the incident 
nucleus, the number of black and grey prongs emitted, NV ,, the number of 
minimum (g<1-49)) tracks emitted, n,, the number of fast a-particles 
emitted, n,, and the number and type of heavy nuclei emitted. In each 
interaction where helium or heavy nuclei were emitted we have measured 
the nee angles between their directions of motion and that of the primary 
nucleus. 
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§ 3. EXPERIMENTAL Resutts 


3.1. Interaction Mean Free Paths 


The interaction mean free paths for the four different charge groups are 
shown in table 1. The results are given in g/cm", assuming that the 
density of emulsion was 3-83 g/cm. Also shown in this table are the 
values found by Noon and Kaplon (1955) for two groups of nuclei 
defined by them as medium, M’, 6<Z <10; and heavy H’, Z>10. The 
mean charge of the nuclei in the M’ group of these authors is the same as 
that observed by us in our M group. 

It can be seen from this table that our results are in good agreement 
with those of Noon and Kaplon. However, we believe that our values 
can represent only upper limits to the true values, since some interactions 
may have been missed. For example, interactions of the type where one 
evaporation track only is formed, and the incident nucleus emerges 
apparently unchanged, are difficult to detect, as are those where no 
evaporation tracks are formed, but the incident nucleus loses a few charges 
in the form of minimum tracks. The detection of this last type of imter- 
action will be very dependent on the relative change of charge. The 
values found for fragmentation probabilities involving only a small change 
in charge are therefore dependent on the efficiency of detection of these 
interactions. For this reason we have not included anywhere in the 
analysis of our results interactions of the first type. The exclusion of 
these events is supported by the fact that they may be examples of the 
diffraction scattering of the hydrogen in emulsions by the incident 
nuclei, and thus represent a different physical process to the other inter- 
actions. 

In order to extrapolate the values found in emulsion to those in air we 
have used a method originally due to Bradt and Peters (1950), which was 
also used by Noon and Kaplon. The reader is referred to this last paper 
for details. Essentially this method requires that some overlap of the 
nuclear volumes has to occur before an interaction takes place, so that the 
interaction cross section o, is given by 


o=n(R, +R, —24R)? and R=r,A¥3(r,=1-45 x 10-% em) 


where R, and R, are the radii of the target and incoming nuclei respectively, 
A is the atomic weight, and 42 has the empirical value of 0:85 < 105m, 
Using this expression, and the known composition of emulsion, a mean 
free path can be calculated for the various groups of nuclei. The values 
found by us and by Noon and Kaplon are shown in table 1. It can be 
seen that they are in good agreement with the experimental values. 
Therefore, we have used the above expression to calculate the mean free 
path in dry air, and the values found are also shown in table 1. In order 
to make these calculations we have had to determine the mean charge of 
These mean charges were determined from the 


h group of nuclei. 
Se s rather than 


number of nuclei found capable of producing interaction 
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from the number that actually produced interactions, and are shown in 
table 1. 


Table 1. The interaction mean free paths in emulsion found in the 
present work and by Noon and Kaplon (1955), together with 


calculated values in emulsion and air 
ee eee 


Present work Noon and Kaplon 
Z glem’ @ Bie 
emul. (exp.) | 51:6+6:1 } 

L nuclei 4:2 | emul. (cal.) 60:1 L nuclei | — emul. (exp.) 61-74 19-4 
air (cal.) 34-2 air (cal.) 31-5 
emul. (exp.) | 51-938 ‘| emul. (exp.) | 59-6+6-0 

M nuclei 7:0 | emul, (eal.) 49-7 M’ nuclei | 7-0 emul. (cal.) 56-5 (2?) 
air (cal.) 27-6 air (cal.) 26-5 
emul. (exp.) | 37-8+5-0 

H nuclei } 12-8 | emul. (cal.) 38:3 " 
air (cal.) 20-6 emul. (exp.) 36-5+4:-8 

H’ nuclei | 15 emul. (cal.) 36-5 

emul. (exp.) | 45-1£9:-5 air (cal.) 18-0 
VA nuclei | 23-6 | emul. (cal.) 284; 
air (cal.) 14-8 


3.2. Angular Distribution of Secondary Nucler 


In those interactions where the geometrical conditions were suitable, 
the space angle was measured between the directions of the primary 
nucleus and those of the emitted helium and heavy nuclei. The angular 
distributions for various classes of interaction are shown in the figure. 

From the data shown in this figure we have estimated the root mean 
square (r.m.s.) transverse momentum of the emitted «-particles in the 
C-system of the incident heavy nuclei. In order to make these estimations 
we have used the results of a previous experiment made in these emulsions, 
(Fowler and Waddington 1956) which showed that the integral energy 
spectrum was proportional to H,-'*, where H, was the total energy ; and 
that the cut-off energy at the top of the atmosphere was 1-55 Bev per 
nucleon. In addition, we have assumed that the fragmentation proba- 
bilities and characteristics of the interactions in the C-system are indepen- 
dent of the energy of the primary nucleus. In order to allow for the loss 
of energy bv ionization in air and emulsion above the interactions we have 
calculated the ‘typical’ depth at which interactions produced by the 
nuclei of different charge groups occur. We find this typical depth for 
L nuclei to be equivalent to 11-5 em of emulsion, for M nuclei 11 cm of 
emulsion, and for H and VH nuclei 10 cm of emulsion. These values 
have been calculated from the observed positions of the interactions, 
allowing for the effects of production of secondary nuclei in the overlying 
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atmosphere, and for the slightly differing cut-off energies of those nuclei 
for which 4/Z 2. 

From the above data we have calculated numerically ‘ effective ’ values 
for the incident momentum p,,, which would result in the same r.m.s. 
angles as those observed experimentally. The values of this ‘ effective ’ 


ALPHA PARTICLES HEAVY FRAGMENTS 


NUMBER OF ALPHA PARTICLES 
FRAGMENTS 


NUMBER OF HEAVY 


3 4 Ge a8 
SPACE ANGLE IN DEGREES 
istributi i fragments emitted from 
The angular distribution of the a-particles and heavy 
= interactions of VH, H, M and L nuclei. 


Table 2. The transverse momenta of the «-particles emitted 
from interactions. For explanation see text 


«-particles Heavy fragments 


L M H VH M H* 
27 2920 
3220 3010 2740 
Fe. 5, ms " z 2 a 0-06 L1LLO-13 
ve degrees 1-1-£0-14 1-47+0-11 | 1:83+0-18 Sea a 0-48+ 
p4(Mev/c) 64-8 83.6 ae = 
Hy (Mev) 13 22 


——————— | __ a arEEEne 
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momentum are shown in table 2. In this table we have separated the 
interactions producing «-particles into those produced by L, M, H, and 
VH nuclei. In each case we give the number of «-particles n,, the r.m.s. 
angle 1/3, in degrees, the r.m.s. transverse momentum /,, in Mev/c per 
nucleon, and the kinetic energy in the C-system of the incident nucleus 
E;, in Mev, assuming isotropy of emission in this system. The errors 
shown on the values of 4/#? and p, are standard deviations calculated from 
the observed distributions. For the sake of completeness we also give 
the values of n and 4/% for interactions of M and H* (H#*=H+V#H) 
nuclei which produce secondary heavy nuclei. 

Within the limited statistical weight of these results two main features 
are apparent. Firstly the values of the mean energy of emission in the 


C-system appear to be considerably greater than those assumed previously ~ 


for such interactions, see for example, Peters (1952). Secondly the r.m.s. 
transverse momentum increases with increasing primary charge. As a 
consequence it appears that primary energies estimated from the opening 
angles of «-particles produced in such interactions have frequently been 
underestimated, and that these energy estimates depend on the charge of 
the primary particle. It may be noted that these experimentally observed 
transverse momenta are in good agreement with the results of Perkins 


(1950) for interactions with large excitation energy on silver and bromine 
nuclei. 


Table 3. The effect on the transverse momenta of separating the emitted 


a-particles into those which are emitted alone, ‘ single « ’, and those 
produced together with other «-particles, ‘ multiple « ’ 


Single « Multiple a Single a Multiple « 


45 76 22 67 


Ng 

o/#* degrees 1-82 +0-20 1-18+0-12 1-50 +0-20 2-50 +0-03 
 (Mev/c) 102411 6647 77+10 127+15 
E;, (Mev) 35 14 20 54 


We have also separated the interactions producing «-particles into those 
with one «-particle, and those with more than one. The results for M and 
H* nuclei are shown in table 3. It can be seen that for M nuclei the r.m.s. 
transverse momenta when only one «-particle is removed is considerably 
greater than when more than one is produced ; whereas for H* nuclei the 
situation is reversed. It appears, therefore, that considerable care must 


be taken in assuming what transverse momentum is appropriate to a given 
case. 
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3.3. Fragmentation Probabilities 


The characteristics of each of the 317 nuclear interactions observed so 
far are shown in table 4. These interactions are listed according to the 
charge of the primary nucleus. It should be noted that the distribution 
in charge in this table does not represent the charge distribution observed 
to enter the emulsions, even allowing for the differing mean free paths, 
since a number of the interactions have been produced by fragments from 
other interactions. In addition, in order to increase the statistical weight 
on the higher charge nuclei, the length criterion was sometimes reduced 
for these, although never below 2 mm per emulsion. 

In table 5 are shown the fragmentation probabilities obtained in this 
experiment, together with those obtained by Noon and Kaplon, and by 
Gottstein (1954). The errors shown on these values are standard deviations 
assuming a normal distribution for the fragments, and using the observed 
distributions for «-particles. Noon and Kaplon found their interactions 
in stripped emulsion stacks by following the tracks of heavy nuclei having 
lengths, in 400 emulsions, of greater than 750, or, in one stack, of 
greater than 450. Gottstein detected his disintegrations by making a 
direct scan for them in glass backed emulsions 400 1 thick. The maximum 
inclination to the plane of the emulsion accepted in this experiment is not 
specified. 

The results shown in table 5 suggest that the fragmentation proba- 
bilities into H or M nuclei found in the various experiments are not 
significantly different. However, those for fragmentation into LZ nuclei 
or «-particles appear to differ widely, with the values found in the present 
experiment being considerably less than those found previously. In 
addition, the probability that no nuclei should emerge from an inter- 
action is found to be considerably higher than suggested previously. 
The divergence between our results and those of Noon and Kaplon is 
exemplified by the frequency of occurrence of interactions in which two 
or more heavy fragments are emitted. Noon and Kaplon, in 214 inter- 
actions, found 11 where two fragments were emitted, and 3 with three 
fragments. We, on the other hand, in 249 comparable interactions find 
only 2 where two fragments were emitted. We have never in this labora- 
tory seen an interaction from which three fragments emerge. : 

This discrepancy between our results and those of the earlier workers 
can be explained in two ways. Either we have failed to detect a 
considerable proportion of the ZL nuclei and «-particles emitted from 
these interactions, or the earlier workers have wrongly identified 
evaporation particles from the target nuclei as energetic fragments. It 
cannot be explained by assuming that the charge determinations on the 
primary nuclei were seriously in error in any of the experiments. For 
example, if we make the extreme assumption that in the present experi- 
ment the fragments identified as carbon nuclei are, in fact, boron, and that 
all the boron primaries producing fragments are carbon, then the fragmen- 
tation probability for M, H and VH nuclei into L nuclei is still only 
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0-21+-0-03 as compared with the value of 0-37-L0-04 found by Noon and 
Kaplon. 

We believe that it is extremely unlikely that we could have missed an 
appreciable proportion of the fragments, since we have examined all 
possible tracks emitted within a projected angle of 10° to the direction of 
the primary nucleus, and we have not found any fragments with a solid 
angle of more than 9°, see the figure. All our primary particles had lengths 
per emulsion of more than 2 mm, and the majority had more than 4 mm, 
whereas previous workers have used considerably steeper tracks, and under 
these conditions the identification of the secondary fragments might well 
be subject to considerable uncertainties. In addition, the fact that the 
interaction mean free paths which we observe are in good agreement with 
those found by Noon and Kaplon suggests that the efficiencies of detection 
of interactions were comparable in the two experiments. We therefore 
suggest that the fragmentation probabilities found by Gottstein and by 
Noon and Kaplon for the production of Z nuclei and «-particles are too 
high. 

In order to obtain an estimate of the fragmentation probabilities in air 
from those observed in emulsion we have adopted the following procedure. 
Nuclear emulsion consists of hydrogen; carbon, nitrogen and oxygen; 
and silver and bromine; while air consists of nitrogen and oxygen. 
Unfortunately, except in special cases, it is not possible to identify the 
target nucleus involved in an interaction. However, we may divide the 
observed interactions into two classes on a basis of the number of black 
and grey prongs emitted. We have assumed that those interactions 
characterized by N,>7 are due to moderately central collisions with 
silver or bromine nuclei, while those with N,<7 are interactions with 
light nuclei and hydrogen or peripherial collisions with heavy nuclei. 
It appears likely to us that these interactions will result in somewhat 
larger fragmentation probabilities than typical interactions with nitrogen 
and oxygen. Therefore, the fragmentation probabilities found for inter- 
actions with N,<7 will presumably be upper limits to those in air. These 
probabilities are shown in table 6 for both our and Noon and Kaplon’s data. 

While interactions with hydrogen can sometimes be identified from 
their appearance we have not attempted to give the fragmentation 
probabilities, as the present statistical weight does not justify our doing so. 
We hope that further work will enable us to obtain these important figures. 


§ 4. CONCLUSIONS 


Our main conclusion from this work is that the probability of a heavy 
nucleus fragmenting into a 1 nucleus is considerably less than has been 
claimed previously on the basis of similar experiments. it 80, then those 
workers who have used the earlier fragmentation probabilities will have 
underestimated the abundance of L nuclei at the top of the atmosphere. 
For example, in an experiment by Kaplon et al. (1954) these workers 
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obtained a value of 0-46 for the ratio of the fluxes of Z nuclei to M nuclei 
at the top of the atmosphere. Using our figures for the fragmentation 
probabilities this ratio is increased to about 0-80. 
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ABSTRACT 


Occurrence of secondary flow of a viscous fluid flowing under a pressure 
gradient in curved channels has often been explained by invoking the 
existence of boundary layers before the primary flow is established. 
In this paper it is shown that instability can occur even after the full 
establishment of the primary flow. Chandrasekhar’s method is used 
in the analysis and the relationship between the stability parameter 
and the wave number of the disturbance for neutral stability is obtained. 


§ 1. INTRODUCTION 

WHEREAS the occurrence of spirals in fluid flow under a pressure gradient 
in curved channels is common knowledge, explanations of the 
phenomenon have been based upon the existence of boundary layers 
before the primary flow is fully established, and the fully developed 
two-dimensional flow has generally been believed to be free of such 
spirals (Goldstein 1938). In the present paper, this phenomenon will 
be analysed by Chandrasekhar’s (1954) method under the assumption 
that the primary flow is two-dimensional and between two concentric 
circular cylinders of small spacing, and that it is already fully established. 
The object is to obtain definite information about the stability of such 
a flow against the formation of Taylor—Gortler vortices of various 
transverse wavelengths. 


§ 2, VELocIty DISTRIBUTION IN THE PRIMARY FLOW 
Cylindrical coordinates (7, 4, 2) will be used, with the z-axis coinciding 
with the axis of the cylinders. If the radius of the inner cylinder is 
a and that of the outer one is b, the distribution of the velocity V of the 
primary flow can be found by solving the Navier-Stokes equations and 
is given (Goldstein 1938) by 


V r Ob Y - 
pia ee ae eS Sire erie. Bee oe Le it 
z-4(; s) +5 mG (1) 
in which mae hie 
a?(In a—In 6) > 0 OL 2 
Al apie 5 KeE= Du Od . . . . ( ) 


uw being the dynamic viscosity and P the pressure for the primary flow. 
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If the spacing d of the cylinders (or b—a) is very small compared with 
the radii of the cylinders it can be shown from eqn. (1) that the velocity 
distribution of the primary flow is essentially the same as that for plane 
Poiseuille flow, i.e. 

d2ncr Kd? 


Jed oe CP ee ee (Eee Meee mr Ll 
in which 
r—a 
oe : ae oe eo. 
f=" aca! 
The mean velocity is, to the same degree of approximation, 
a7 
P= ~ Gok 


so that the Reynolds number based on J,,, and d is 
tee) nee 
ye = | brat 


R (5) 


in which v is the kinematic viscosity. From eqn. (3) it follows that, 
with terms of higher order in d neglected, 


V Ka : 
~=—e (Pf gy oe 
v4 2 SE (xg—1) Ah er 


in which the prime indicates differentiation with respect to r. These 
quantities will be useful later in the investigation of stability. 


§ 3. FORMULATION OF THE STABILITY PROBLEM 


To investigate the stability of the primary flow, one may superpose 
on it a time-dependent perturbation with (small) velocity components 
(u, v, w) in the directions of r, 6, and z, respectively, and with corresponding 
pressure p. If it is kept in mind that the primary flow is such that the 
Navier-Stokes equations and the equation of continuity are satisfied, 
and if terms of higher order than the first in the perturbation quantities 
are neglected, the Navier-Stokes equations are 


Ou 2Vv_  1édp ee 

ov we Va v 

ape ae al aM (v- *): eae as OLY (9) 
Ow —ildp A 
FE pemmararr ean ees 3 KG) 


in which ¢ and p denote time and density respectively, and 
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The equation of continuity is 


O(ur)  d(wr) 
Now, following Taylor (1923), one may assume 
(wu, v, w)=exp (ct)[w4(r) cos kz, v,(r) cos kz, w,(r) sin kz]. . (12) 


If p is eliminated between eqns. (8) and (10), and if eqn. (12) is used, 
one has 


v ; o 2V i! 
EP? (V2—18— *) w4=— a Cav (vie—ie— 2 — *) Uz, . (13) 


in which 
D=d/dr, Vi=Di+ — Dp. 
Furthermore eqn. (9) can now be written as 
an a ae 
ui(V +2) =» (ve-e—5—*) Vy cline si) ots (14) 


and the equation of continuity is 
Du, +hw,=0. ieee a te oa aad at 


Although for neutral stability o may be purely imaginary, one may 
follow Taylor in taking it to be zero—at any rate if one is looking for the 
formation of Taylor—Gortler vortices. For neutral stability, therefore, 
eqn. (14) becomes 

lp) me Vite cee 3 a fey CLG) 
in which 


1 I 
= [2 a aly Gs, 
L,=D*+ -D—B—-, 


and it can be readily shown that eqns. (13) and (15) yield 
7 Ly%,= — oy otis. ie, "oa aan. (La) 
Since the spacing d is much smaller than a or 6, and since the perturba- 
tion velocities change appreciably as r varies between a and b, and the 
wavelength (27/k) is comparable with d, the part 
1 ch 
ro 
in the operator L, is small compared with the rest of the terms. If 
now one denotes kd by « and expresses w, and v, in terms of V,,, eqns. (16) 
and (17) become, in virtue of eqns. (5), (6), and (7), 


Kd 
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7 1,26 
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in which 
L=C/d2—a?, p=d/a. 
If one sets 
Uy= — 120? RBW We ge os (4) 
one has finally 
L2W=(&—£)v, 1 te A eee 
Dv,=SeX(2E—1)W ss ea) 
where 
OST ee (23) 


is the stability parameter. Equations (21) and (22) are to nee with 
the boundary conditions 

U(O)\=0, (1) =O W(0)=W(1)=0, W’'(0)=W'(1)=0 . (24) 
the last two being a consequence of eqns. (15) and (20) and of the fact 
that wu, and w, are zero at €=0 and =1. The primes in the last two 
of eqns. (24) indicate differentiation with respect to €. The differential 
system consisting of eqns. (21), (22), and (24) represents an eigenvalue 
problem, the solution of which will result in a functional relationship 
between S and « for neutral stability since c has been assumed to be zero 
to start with. 


§ 4. OUTLINE OF THE METHOD OF SOLUTION 


Since the method of solution used in this paper is that of Chandrasekhar 
and is already described in detail elsewhere (Chandrasekhar 1954), only 
an outline of the essential steps will be presented here. First, one takes 


v= Pa EEN mare er et rey (O15) 


to satisfy the first two conditions in eqn. (24). Then this expression for 
v, 18 aa into eqn. (21) which can be solved to yield the result 


Ve 2 STE ii| Bs cosh «€+C,,, sinh «€+ D,,é cosh «€+ E,,€ sinh a€ 


: 4mm(2E—1) cosmm€ | 4(M—6m?n?) sin mr 
+(@—€) sin mag “ree 
(26) 
in which M=m?r?+«? and the constants B to D can be determined to 
satisfy the last four conditions in eqns. (24). With these constants so 
determined, eqn. (26) is then substituted into eqn. (22) and the right-hand 
side expressed in a Fourier sine series. Then coefficients of the individual 
harmonics at left and right are equated. This result in an infinite number 
of homogeneous linear equations involving the A’s. Finally the eliminant 
of these equations (or the condition that the A’s are not all zero) furnishes 
the desired relationship between S and «. 
Although the eliminant is a determinant of infinitely many rows 
and columns, Chandrasekhar (1954) has shown that often it is sufficiently 
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accurate to equate the first element (in the first column and first row) 
to zero, and that the result of the second approximation (by equating a 
two-by-two determinant to zero) does not differ much from that of the 
first. The same was found to be true in a previous calculation (for a 
slightly different primary flow) by the authors whenever the second 
approximation was tried for spot checking, so that the relationship 
between S and « determined by the first approximation is believed to be 
sufficiently accurate. The ease with which this relationship was obtained 
is additional proof of the power of Chandrasekhar’s method in solving 
problems of the present type. 


§ 5. RESULTS aND ConNCcLUDING REMARKS 


The desired relationship between S and « for neutral stability of the 
primary flow is given graphically in the figure. It is seen from this 


0 1000 C000 3000 4000 5000 


Relationship between the wave number a and the stability parameter S for 
neutral stability. 


figure that established laminar flow under a pressure gradient in a narrow 
curved channel is unstable for values of the parameter S greater than the 
approximate critical value 600, which corresponds to a wavelength of the 
Taylor-Gortler vortices equal to 27/2-25 or 2-79 times the width of the 
channel. Since the flow tends to be unstable as S is increased, for the 
same Reynolds number and the same (small) width of the channel it. 
does so as the curvature is increased, as can be seen from the definition of S, 
Thus centrifugal force is as much a cause of instability for the established 
flow as it is for the unestablished one, and the occurrence of secondary 
flow cannot be considered to be exclusively a consequence of boundary- 
layer development. 

Attention might also be called to the fact that, since centrifugal force: 
can be considered as a form of gravitation (Mach, Hinstein), the non- 
uniformity of it is a kind of stratification—a stratification not of density, 
for the density is homogeneous, but of specific weight in the general 


YZ 
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sense. Secondary flow occurs if this stratification is not stable, and 
the mechanics of it is not unlike the falling of a heavier fluid in a lighter 
one, particularly if the boundary conditions are alike, and diffusion takes 
place to render the variation of density and that of centrifugal acceleration 
strictly comparable. This point of view is substantiated elsewhere by 
the equality of the Rayleigh number for thermal instability and the 
Taylor number for rotational instability (both being 1708), as first 
pointed out by Low (1929). 
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THERE is considerable experimental evidence that the lower lying levels 
of nuclei of masses 25, 29 and 31 may be rather well described as over- 
lapping systems of rotational bands (Litherland et al. 1956). The 
moments of inertia required are close to those indicated from a similar 
interpretation of the levels of neighbouring even—-even nuclei. The 
deformations required are in the range B=0-2-0-3 which are consistent 
with the deformations indicated by the few measured quadrupole moments 
in this region. One would not have thought that the collective model of 
Bohr and Mottelson (1953) which has had such marked success in the 
heavy nuclei would be useful for such light nuclei where the zero point 
oscillation amplitude might be expected to exceed the deformation value. 
This would prevent the rotational states from being developed to any 
significant extent. However when the detailed properties of these nuclei 
such as relative and absolute transition probabilities and level ordering are 
compared with a model such as that of Nilsson (1955) who has computed 
the single particle energy levels in a deformed potential, the agreement 
is remarkably good. 

Considerable success has been achieved in the region around 10 in 
describing the nuclei in terms of configuration mixing in an intermediate 
coupling shell model; it is of interest therefore to explore the two 
approaches through a direct comparison for the same nucleus. The 
shell model calculation becomes impractically laborious in its rigorous 
form for nuclei above 2°Ne having more than four nucleons outside the 
16Q closed shell. The nucleus !°F has been very successfully treated on 
the intermediate coupling shell model by Elliott and Flowers (1955). 
Some modifications arising from collective effects have been considered by 
Barker (1956). This paper is an attempt to apply the rotational model 
to this nucleus in a similar way to that used for ?°Al. No attempt has 
been made to introduce refinements to the theory which may be necessary 
for nuclei in this region. Such refinements are being considered by 
Rakavy (private communication). 


The Even Parity States of 1°F 


Figure 1 shows the known states of °F below 6 Mev together with some 
of their measured properties (Ajzenberg and Lauritsen 1955). Let us 
consider first the even parity states. Based on the ground state with 
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I=} we expect a rotational band with K=}, where K is the projection of 

the total angular momentum on the nuclear symmetry axis. The 

level order for bands with K=} is given by (Bohr and Mottelson 1954). 
HE, =Ey+h/28{1(1 +1) +a(—Y +9) 

where S is the moment of inertia and a depends on the configuration mixing 

of the last odd particle. For reasonable values of a given by Nilsson’s 

model and of the moment of inertia, the level spacing is much greater than 


that observed, as shown for example on the left-hand side of fig. 2, although 
the level order is correct. 


Fig. 1 


The known states of 19F below 6 Mev. The arrows from the left represent 
B-transitions from 1°O, that from the right the f-transition from 19Ne. 
All three are labelled with the value of log ft. 


Let us now consider the next higher rotational band which may occur. 
According to Nilsson’s model this will have K=3 and will be about 3 Mev 
higher. Kerman (1956) has shown that through the rotation—particle 
coupling there is the possibility for K-admixture between bands with 
K differing by one unit. This means that K is no longer a very good 
quantum number. The level ordering is affected since states having the 
same spin in the two series will be pushed apart. In this case the state 
with J=} is not changed since it has no corresponding member in the 
K=% band. The mixing amplitudes and the changed level positions 
can be calculated using the matrix element for the interaction between the 
two bands which can be derived from Nilsson’s wave functions. 
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The final fit shown in fig. 2 has been derived by requiring (1) the 3 state 
with K =} to be at 1-55 mev where it is observed and (2) the 3 state with 

=} to be at 4 Mev near the first two states one of which could be 
identified with it. For a value of B=0-3 the agreement shown in the 
third and fourth columns of fig. 2 is obtained. The rotational energy 
i?/23 required for both bands is 0-3 Mev. When the first excited states 
of 2+ of 18O and 2°Ne are assumed to be rotational states, values of 0-33 
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and 0:27 mev for /2/25 are obtained. The 8 state occurs " ues ae 
to be compared with the observed position of 0- 197 mev while a 3 i : 
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contain 46% of K=3. 
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Using this specification for levels, other properties may be derived from 
the wave functions of Nilsson. These are listed in table 1 along with the 
experimentally determined values. The corresponding predictions of the 
shell model calculation of Elliott and Flowers are also listed. In the 
calculation of the H2 lifetime of the 3+state on the rotational model, the 
collective as well as the particle part of the #2 transition probability has 
been included. For 19F the two parts are of comparable importance, 
unlike the heavy nuclei where the collective part predominates. 


The Beta Decay of 1°O 


To explain the f-transition probabilities for the decay of '°O to the 
states of 19F, it is necessary to interpret the states of 1°O also on the 
rotational model. This is difficult because of the meagre experimental 


Table 1. Predictions of the Two Models compared with the 
Experimental Values 


Rotational Shell 


model miodel Experimental 


Property 


y-ray branching 
Ratio $+ 3+ 
to $> 3+ 0:8% 0-6% 40/,(@) 


Magnetic moments : 
of 4+ state 2-75 n.m. 2-80 n.m. 2-63 n.m.@) 


of 3+ state 3°70 3-3 3°50 +0-240) 


Lifetime of $+ 
state 0-9x10-7 sec |5+10-7 sec] 1-25+0-025+-10-7 sec 
1+10-? sec 


Log ft value for 
19Ne decay 


(a) Ajzenberg and Lauritsen (1955). 

(b) Mack (1950). 

(c) Phillips and Jones (1956). 

(d) The first value is that obtained when core vibrations are neglected, the 
second value is obtained when such vibrations are included. 


(ec) Johnson (1956). 
information available. The spin of 1°O should be 3+, K=3 on Nilsson’s 
model but experimentally may be $ or 3+. The spins and positions of 


other states in '°O need to be determined before an unambiguous descrip- 
tion of the 6-decay is possible. 


The Odd Parity States of °F 


A state of $— has been identified at 0-110 Mev in 19F. In addition two 
states have been identified (Toppel et al. 1956, Freeman 1955) at 1:35 and 
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1-46 Mev which are not accounted for in the interpretation overleaf of the 
even parity states. These might be identified as /=1, 3 and 3 members of a 
K=}4 band resulting from a hole in the p-shell. From the position of the 
first excited states of 160 one would expect that 6 or 7 Mev would be re- 
quired to break up the p-shell. Theoretically, however, it appears 
that for large enough deformations the rotational energy may be great 
enough to cancel out this energy for the break-up of the shell and so bring 
an odd parity particle level down near the ground state. The observed 
y-ray branching from the states at 1-35 and 1-46 Mev (assuming them to 
have odd parity) is characterized by the weakness of the #1 transitions 
compared to the 1 or £2 transitions. This discouragement of the #1 
transitions would also result from a large deformation. On the other hand 
the spacing of the levels indicates a fairly small deformation. This 
inconsistency makes it difficult to draw further conclusions. However, 
this model like the shell model requires the two states at 1-35 and 1-46 Mev 
to have odd parity. 


Comparison with the Shell Model 


From table 2 we can see that the predictions of the two models are 
remarkably similar. In table 2 a further comparison is made. The 


Table 2. Configuration of the Last Odd Particle for the Even 
Parity States according to the Two Models 


State 4 Rotational model Shell model 

(%) (%) 

$ 4 29 49 
: 3 9 9 
$ 62 42 

3 = 1 15 21 
3 9 19 

5 76 60 

8 } 23 19 
3 9 10 

5 68 71 


probability of finding the last odd particle in a given configuration j is 
listed for each even parity state. Those derived from the intermediate 
shell model calculation are in the coordinate system fixed in space, while 
those derived from the rotational model are in the coordinate system fixed 
inside the deformed nucleus. The results are comparable if the period of 
the rotational motion is slow compared to the particle motion. If this is 
the case the mixtures of configurations derived from the two approaches 


show similar trends. 
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ABSTRACT 


Natural growth features observed on crystals of lead nitrate when 
growing from aqueous solution are described in detail. The modes of 
formation and behaviour of growth layers are discussed. A special 
feature is the frequent occurrence of stacking fault lines over which the 
growth layers cannot pass. These often appear at centres of growth. 
The effect of methylene blue as impurity in solution restricts the growth 
on the cube faces in a direction parallel to a diagonal in the face. This 
results in the formation of small striations or tilt faces. 


§1. INTRODUCTION 


FoLLow1nc the interest in spiral layer growth as predicted by Burton et al. 
(1951), in their theory of crystal growth, it was decided to re-examine the 
growth of some of the crystals reported by Bunn and Emmett (1949) as 
showing layer growth. Bunn and Emmett used cine-photomicrography 
with an optical microscope, a method by which layer growth can be 
studied as it occurs and a record of the variations with time made. It 
is advantageous to use a high magnification, which somewhat restricts 
the condition of growth of the crystal. It cannot be suspended in a freely 
stirred medium, and the faces which can be examined depend both on the 
habit of growth, and on which face the crystal lies. If it is hoped to see 
any layer growth, the step height and the refractive properties of the 
solution and crystal will determine whether this is possible. The speed 
of growth of the layers will also affect visibility. If growth is rapid 
even statically visible layers will become blurred. 

Bunn and Emmett describe the growth of a few crystals in detail and 
others are just mentioned as either showing or not showing layers. They 
reported frequent initiation of layers at points inside the crystal faces, 
which was puzzling in the light of theory at that time. Information on 
the nature of the centres of growth was not given. = 

Forty (1951, 1952) examined the growth of cadmium iodide (one of the 
crystals studied by Bunn and Emmett) in greater detail, showing that it 
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does grow as a result of dislocations in the crystal at later stages of 
growth. He has made determinations of step heights interferometrically. 
Korndorffer et al. (1952) have shown that spiral growth of layers may result 
from deformation of previously perfect crystal platelets of cadmium iodide. 
Both screw and edge dislocations were produced and the growth rate was 
greatly increased in the deformed crystals. 

Spiral layer growth on crystal faces has now been reported for several 
crystals and measurement of step height has been carried out using multiple 
beam interferometry which is well suited to this type of examination of 
crystal surfaces. This procedure however is only easily done with crystals 
which have ceased growing. It is mostly on the basal planes of thin 
plate-like crystals that these features have been seen. Such crystals 
have grown slowly in the direction normal to the basal plane but rapidly 
in extension sideways. The initial thin plate formed would be expected to 
undergo strain and this would give rise to dislocations as suggested by 
Frank (1951). Such behaviour would not be expected to occur so readily 
in a cube or octahedral form crystal which tends to grow more uniformly 
in all directions. Hence a study of the dynamic growth of lead nitrate 
has been made. Lead nitrate belongs to the cubic system and shows 
clear layer growth on octahedral faces when growing from pure solution. 


§ 2. EXPERIMENTAL DETAILS 


A saturated solution of lead nitrate was prepared using Analar materials 
and distilled water. The solution was warmed above room temperature 
and then allowed to cool, excess of the salt having been added, so that a 
supersaturated solution was obtained. A drop was placed on a thin cover 
slip on the stage of the microscope. The crystals grow as octahedra with 
one of these faces lying in contact with the cover slip. The face above 
this was then observed during growth. The crystals appear when the 
solution has reached a labile condition. If the solution is left untouched, 
dendritic crystals tend to form from the edge of the drop. Crystals could 
be nucleated by scratching the cover slip in the centre of the drop when 
the solution was in the labile state, a crop of small crystals forming. The 
growth features on crystals nucleated in this way were the same as those 
on naturally nucleated crystals. 


§3. MrrHop oF OBSERVATION 


Using an inverted type microscope, the crystals were studied in trans- 
mitted light, with phase contrast illumination. Light traversing a 
crystal in solution with a layer growing on its upper surface will undergo 
a differential phase change according to which side of the layer it passes 
through the crystal, corresponding to the thickness of the layer. The 
phase difference will depend on the relative refractive indices of the solu- 
tion and the solid crystal and the thickness of the step. With transmitted 
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light, positive phase contrast takes place and the top of the step appears 
dark, the lower side appears bright. 


§4. Narurat Layer GrowtH on Leap Nirrare Crystats 


The natural forms of layer growth on lead nitrate crystals may be 
described as follows. Immediately after nucleation of the crystal the 
growth is very rapid due to high supersaturation. On the octahedral 
face the visible layers are very irregular and optimum illumination and 
adjustment shows closely spaced smooth layers of very small step height, 
moving rapidly, sweeping over the face in about one second in the region 
between the thicker layers, continually building these up and breaking 
them down. It is unfortunately impossible to photograph these layers. 
Figures 1, 2 (order of magnification x 600) show successive photographs 
at an interval of about one second of a face on which this is occurring. 
The larger layers are visible terminating in the face, constantly altering 
in length and thickness. The projected film shows the effect clearly. 
Growth is very rapid at this stage and nucleation of layers occurs 
apparently at random at different points on the crystal face. The shape 
of the thick layers is generally dendritic in form. These conditions would 
probably be normal with faster growth at a higher temperature. The 
probability for misfits to occur must be high at this stage. 

When the growth conditions become steadier, after the first few 
minutes, thick layers predominate, through the small layers bunching 
into the thick layers and these are now continuous over the surface. 
They may be found either to move out from centres of growth, figs. 3, 4 
(x 600), or they may move inwards forming a pit of encirclement as 
described by Bunn and Emmett, figs. 5, 6 (300). The lower side of a 
step in these photographs is the bright edge, and the dark side is the top 
of the step. The time taken for a layer to cross the crystal surface was 
of the order of ten seconds. 

From several observations on fresh crops of crystals the natural growth 
appears to be that of the Kossel-Stranski model of layers nucleating 
preferentially at the corners and edges and moving inwards. The 
layers advance more rapidly at the edges of the crystal face, and result in 
pits of encirclement in the centre. If dislocations were responsible for 
the nucleation of these layers, the centres would eventually move within 
the face when the crystal extends further outwards. This did not happen. 
Bunn suggested that these pits might be related to the trigons observed 
on the octahedral faces of diamond, but these pits never show a precise 
orientation with sharp edges as is the required property of a trigon, and 
their position and shape varies with successive layers. They occur as the 
natural sequence of the enclosure of layers and have no relation to trigons. 
However, as described later, some features akin to trigon formation have 
been observed on the crystals. : 

Bunn and Emmett attributed growth from the edges of the face to 
overlapping of a layer from a neighbouring face. A layer reaching an 
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edge might wrap over on to the next face. This effect has been studied in 
relation to the present work and no justification for this view has been 
found. An instantaneous photograph might give this impression but it 
was always seen that in fact there was no relation between layers on 
neighbouring faces. The initiation of a layer at an edge might at one 
instant occur at the same time as one was being completed at the same 
edge on the adjacent face, but this did not always happen. 

The growth from centres appears in crystals which are no longer perfect 
and closed loops or single spiral growth, sometimes from more than one 
centre, may be seen. Figures 3, 4 (x 600) show a complex spiral system, 
and fig. 7 (300) a single spiral. Fig. 8 (300) shows the same 
spiral after a short time interval when it has changed to a double spiral. 
Periodic changes of this type were often observed. In fact, the centres of 
growth were not located precisely at the same point throughout, possibly 
because of multiple dislocations being present, distributed over a small 
area. Dissociation and recombination of layers in this region would give 
rise to these effects. The centre of the multiple spiral illustrates this 
point. The visible thick layers must build up from smaller layers by the 
process described by Frank (1949) of the effect of several sources in 
region. The actual centres of the dislocations where they emerge on to 
the surface are not revealed by the visible layers, as these are the result 
of the smaller layers bunching and interfering with one another, the 
process varying with time, and depending on the relative step heights and 
the separation of the centres. The origin of the dislocations causing 
growth from the centres of faces appears to be due to stacking faults 
occurring during the period of early rapid growth of the crystal, rather 
than the buckling mechanism which occurs in thin plate crystals. 

Figures 9-14 (300) show a phenomenon often appearing, that is, of 
one layer overtaking another which terminates in the face, joining with 
it as in fig. 10, and passing on leaving part of the joining layer terminating 
in the face in its place, fig. 11. This then curls around the point of termina- 
tion as it grows further, showing the behaviour of the layer tied toa disloca- 
tion source. This continues until a new layer crosses the region from 
the dominant growth centre at the top of the crystal face where there is 
a multiple spiral source. The tied position of the layer is not constant ; 
this again must result from the building up of layers from small sources, 
existing over small local area as already described. 


§5. Stacking Favutts In GrowrnH PHENOMENA 


Another feature frequently observed in naturally grown crystals is a 
crystallographically orientated ledge which acts as a centre of growth 
often of a spiral nature. Figures 15 (770), 16 (x 600) show such ledges 
with multiple layer growth. The directions of the ledge are always 
parallel to the edges of the octahedral face. The growth layers move 
freely over the face until they arrive at these lines, which they cannot 
cross but pile up against. The ledges must mark the line of a stacking 
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fault in the crystal since the layers cannot cross the line during growth. 
The fault may have only a single direction, i.e. be seen as a straight line 
terminating in the face across which growth is impossible, or have two or 
three joining line sections in the crystallographic directions as seen in 
the photographs. During growth the layers pile up against the line of the 
fault forming the visible step, and then travel round to the other side, 
where they are again held by the fault line and the step is cancelled as 
they pile up on this side of it. On one side the step is apparently grown 
in and on the other grown out. The original fault would be some fraction 
of a unit cell distance of misfit. 

Figures 17—22 (x 600) illustrate a single layer moving round such a 
ledge at different stages of growth. When the layer returns to the first 
side the fault line trace reappears. The inhibition of growth over these 
lines must be due to misfit along them. As growth proceeds with time the 
misfits or stacking fault is propagated through the crystal by the growth 
process described. These faults occasionally traverse the whole face, 
growth taking place freely on either side, with the layers originating from 
the fault line. These steps would be indistinguishable from slip lines and 
it seems feasable that the conjectured ‘slip’ lines observed on diamond 
octahedral faces recently (Emara 1956) are in fact fault traces of this 
nature which have been propagated through the crystal during growth. Some 
such propagation of faults might also explain the formation of trigons in 
diamond. The diamond structure probably gives rise to faulting in three 
planes of stacking fault arising from a point of original faulting, suggested 
by Halperin (1954), which would form during the further growth, as each 
layer would be inhibited at the trigon fault step which would result on the 
growing surface. This would become larger during the successive passage 
of growth layers. Growth would also be possible by nucleation inside the 
pits left on the surface by the holding up of the growth layers, but crossing 
of the fault planes would be impossible. Faulting could occur at different 
times and again within trigons already formed, as is often seen on many 
diamonds. 

Another effect observed in connection with the growth of crystallo- 
graphic fault ledges on the octahedral surfaces of lead nitrate is that these 
are not always stable at the points joining two directions of the ledge. 
One layer on passing this point will build up the ledge, and then, starting 
at the point of junction, break through and ultimately cross over and 
eliminate the ledge. One would imagine that subsequent layers would 
then pass freely over this region, the fault no longer impeding the growth, 
but this is not the case. Restriction occurs again initially with the 
following layer, with the subsequent breakthrough occurring after a 
short time delay, always beginning at the point where the two directions 
of the ledge meet. Each following layer follows the same pattern. 

Figures 23-28 (x 300) show this kind of breakthrough beginning at the 
joining point in fig. 24. The step has disappeared in fig. 26 and is being 


produced again in fig. 28. 
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§6. Errect or Impuriry on LAYER GROWTH 


The effect of impurity on the growth layers on lead nitrate crystals was 
studied with methylene blue ia solution. As is well known, the habit of 
the crystal is then changed to the form of slowest growth on the cube faces, 
as the impurity is presumably preferentially absorbed on these faces, 
retarding the growth. This is visually confirmed by the figs. 29, 30 
(x 600) which show the effect of impurity on the growth layers on the 
cube faces. With even little impurity present cube faces appear, but 
the layers on the octahedral faces are unaffected, and growth continues at 
the same rate and in the same manner as previously described. But the 
layers on the cube faces are strikingly restricted and growth is much slower 
than that on the octahedral faces. Cube faces may occasionally be seen 
in pure solution growth, when the layer form is very much as that of the 
early rapid growth described for the octahedral faces, with building up 
and rapid breaking down of the visible layers which are irregular in shape. 
There is no restriction of the layers. However, with methylene blue in 
solution the layers become thicker and growth is strongly inhibited in a 
direction parallel to one of the diagonals of the face, while still free in the 
direction of the other diagonal. This produces striated growth and a 
fast rate of layer growth along the direction of the striations and a slow 
rate perpendicular to them. The surface becomes stepped with small 
tilts at the striation boundaries, reported by Buckley (1951) as being (451) 
index faces from measurements on a large crystal. Figures 31, 32 (x 600) 
show an octahedral face on a crystal growing in solution with methylene 
blue present. The tilt edges of the striations on the cube faces are shown 
by the intersection line with the neighbouring face. It is seen that the 
position of the tilted region has changed in the time between the two 
photographs moving towards the corner of the cube face as growth 
proceeds. 
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ABSTRACT 


THE present paper investigates by experiment and theory the mechanisms 
governing work-hardening, work-softening and slipband formation in 
face-centred cubic metals. The rapid hardening stage (stage II) of the 
stress-strain curve is investigated by flow—stress measurements at different 
temperatures, combined tensile and torsion experiments, and by the 
observation of the length of slip lines as a function of prestrain. Stage [1 
of the stress-strain curve and the work-softening phenomena associated 
with it are investigated with the electron microscope, showing that the 
principal surface patterns (slip bands, fragmentation) are due to cross 
slip. It is argued that in stage II the slip distance is decreased con- 
tinually by the formation of Lomer—Cottrell dislocations. This also 
accounts for various observations on the end of the easy glide region. 
The temperature dependence of work -hardening in stage III is caused by 
screw dislocations circumventing the sessile Lomer—Cottrell dislocations 
by cross slip. The mechanisms through which cross slip causes glide 
band formation and fragmentation are discussed in some detail. 


§ 1. INTRODUCTION 


RECENT studies of work-hardening of metals have concentrated on single 
crystals of face-centred cubic metals. It has proved useful to distinguish 
three stages of the stress-strain curve which may be characterized as 
follows (fig. 1): 

Stage I is the so called ‘easy-glide’ region. The coefficient of work- 
hardening 0=dr/de is rather small, 6/G (G=shear modulus) being of 
the same order of magnitude as in typical hexagonal metals (Zn, Cd) 
at low temperatures. It is probably justified to consider the work- 
hardening mechanism in stage I as being analogous to the work-hardening 
in hexagonal metals. 

The coefficient of work-hardening of stage I is about one order of mag: 
nitude larger than that of stage I and rather independent of the impurity 
content ; stage IT seems to be characteristic for all face-centred cubic 
metals. The observation that the stress-strain curves of hexagonal metals 
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do not show a stage with a correspondingly large coefficient of work- 
hardening suggests that stage II is connected with the capability of face- 
centred cubic metals to slip on intersecting glide-planes. 

Stage III is characterized by a decrease of @ with increasing strain. 
From the temperature dependence of work-hardening it is obvious that 
in this stage the work-hardening process is partially balanced by a 
thermally activated process which will be called ‘dynamical recovery ’. 
It should be noted that this dynamical recovery proceeds during deforma- 
tion even at temperatures where no static recovery is observed. 

The main object of the present paper is an investigation of the 
mechanisms of work-hardening and dynamical recovery going on in stage II 
and stage III by experiment (work-hardening curves, combined tensile— 
torsion tests, surface observations by ordinary and electron microscope) 


Fig. 1 
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Schematic drawing of the typicai shear stress—shear strain curve of a face- 
centred cubic metal at two different temperatures 7, and T,>7'. 


and theory. We shall be but little concerned with stage I which is 
probably best investigated together with hexagonal metals such as zine 
and cadmium. The small coefficient of work-hardening, its dependence 
on the diameter of the crystal (Paterson 1955, Suzuki et al. 1956)+ and the 
appearance of very long (and rather faint) slip lines are best explained by 
assuming that in stage I the slip-distance of the dislocations is comparable 
with the diameter of the specimen (Mott 1952). We mention in passing 
that contrary to stage II the stress-strain curve in the easy-glide region is 
rather sensitive to structural irregularities (Diehl 1956 b). 


‘ A a comparison of the data obtained by different authors see Diehl 
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§ 2. Tae Mecuanism or Work-Harpenine In Stace II 


2.1. Survey 


The extent of the easy-glide region and therefore the strain €, at which 
stage I begins depends very markedly on the orientation of the tensile 
axis with respect to the crystallographic directions of the crystal. «, 
is smallest if the tensile axis is in a symmetric position. Several authors 
have suggested that the onset of rapid hardening at the beginning of 
stage IT is due to slip in secondary glide-systems} (Masing and Raffelsieper 
1950, Réhm and Dieh] 1952, Liicke and Lange 1952, Haasen and Leibfried 
1952). A rather striking argument in favour of such an explanation is the 
temperature dependence of ¢, (Seeger 1954 a). We shall return to the 
theory of this observation later (§ 4.2). 

A more detailed consideration of the mechanism through which slip 
in secondary glide-systems is able to increase the coefficient of work- 
hardening poses the following questions : 

(1) Is the increased 6 caused by a continuous admixture of slip in 
secondary glide-systems, or is the action of the secondary systems confined 
to the beginning of stage IT only, as proposed by Friedel (1955)? 

(2) Is the main effect of the secondary glide systems due to an increase 
in density of the dislocation forest which must be cut by the dislocations 
moving in the primary glide system (Paxton and Cottrell 1954) or is it 
of a more indirect nature due to the formation of obstacles which cannot 
be overcome by the dislocations in the primary glide-system with the aid 
of thermal energy ? 

These problems have been attacked experimentally by a combination 
of tensile and twisting experiments on Cu single crystals. We preferred 
copper to aluminium because copper shows at room temperature all 
stages of the stress-strain curve fully developed, and because extensive 
data on the stress-strain curve of the same material are available (Diehl 
1956 a). The main difference between the technique of Paxton and Cottrell 
(1954) and ours was that we avoided complications due to the elastic core 
during twisting by using tubular crystals. The external diameter of the 
erystals was 5-5mm; the thickness of the wall was 0:75mm. For 
details concerning the preparation of the specimens reference should be 
made to Rebstock (1956). 

Problem (1) was investigated mainly by measurements of work- 
hardening in tension after twisting (§ 2-2). The basic idea was to intro- 
duce by twisting ‘ artificially ’ additional dislocations in secondary glide- 
systems and to simulate thereby the anticipated effect of glide in secon- 
dary systems in stage IT. . 

Problem (2) was investigated by measurements of the (reversible) 
change in flow stress connected with a change in temperature (§ 2.3). 

ee 

+ The octahedral glide-system with the largest resolved shear stress 7 is 
calied ‘ primary slip-system’; all the other ones are denoted as “ secondary * 
slip-systems. 

Z2 


326 A. Seeger et al. on Work-Hardening and 


Both techniques were supplemented by observations of the surface 
appearance of deformed crystals (§ 2.4). 


. 2.2. Measurements of Work-Hardening 

Figure 2 gives some selected data on shear stress—shear strain curves, 
curve A corresponding to an ordinary tensile test. Curves B and C 
were obtained by tensile tests on crystals with the same crystallographic 
orientation as crystal A. At the strains indicated by arrows the specimens 
were unloaded, twisted without removal from the apparatus and reloaded 
in tension. It will be seen that a single twist increases the flow stress. 
There is no increase in the rate of hardening during the subsequent tensile 
test. Disregarding the behaviour of the stress-strain curve immediately 
after reloading in tension, it may be said that the principle effect of an 


Fig. 2 
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Tensile work-hardening curves of tubular copper crystals at room temperature : 
The origin of curves A is shifted by 4¢=0-1; Band C were shifted in 
addition to this in order to agree at e=c«, with 4. An arrow |, or 4 
indicates an interruption of the tensile test and a twist of the crystal, 
y being the surface shear of the twist. 


intermediate twist in stage Il of the work-hardening curve is a parallel 
displacement of the stress-strain curve towards higher shear stresses by 
an amount depending on the shear of the twist (see dashed part of B 
and C). 

Contrary to what should be expected on the basis of Friedel’s hypo- 
thesis, slip on secondary systems occurring only once does not 
result in an increase of the rate of hardening. Since however 
the jlow stress is increased by such a procedure it may be expected that 
repeated twists and therefore repeated action of secondary glide-systems 
de give an increase in the average rate of hardening. This is verified 
experimentally by curve D in fig. 2. After a shear strain -=0-2 in tension 
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crystal D was repeatedly twisted by the same amount y as indicated in 
fig. 2. The stress-strain curve obtained by the technique as discussed 
above does indeed show an increase of its mean slope after -=0-2. We 


conclude from these experiments that the increased rate of hardening 
in stage IT is due to continual slip in secondary glide-systems. 


2.3. Reversible Change of Flow-Stress with Temperature 

The following discussion will be based on the theory of flow-stress given 
elsewhere (Seeger 1954 a, b, c, 1955 a). The flow-stress of a pure metal 

can be divided into two contributions 7, and tr, according to 
PT Ger Tats eh Sa ae ae SA. ou, aman) 
For a given dislocation pattern 7, depends on temperature only indirectly 
through the temperature dependence of the shear modulus @ (or other 
elastic constants) and has its main origin in the long range stress fields of 
the dislocations of the primary glide-system. 7, is the contribution of 


Fig. 3 
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Division of the stress-strain curve of a copper crystal deformed in tension at 
room temperature into rg and ry. 


the dislocation forest threading the primary glide-plane. Leaving aside 
certain corrections, t, depends linearly on the apeolnte temperature 
T. Figure 3 gives the division of 7 into 75 and 7, as a function of sheet 
strain « for a single-crystal of copper at room temperature. These 
results were obtained (see Rebstock 1956) by measuring the Hove soroes 
after changing the temperature between T,=90°K and 7',=293'K, 
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and were evaluated using the theory of flow-stress as given for Cu 
(Seeger 1955 a) and the temperature dependence of the shear modulus 
as calculated from single crystal elastic constants (Overton and Gaffney 
1955). Figure 3 shows that whereas 7, and r, are of comparable magni- 
tude in stage I the relative contribution of 7, to the flow-stress is sub- 
stantially reduced in stage II. This difference between stage II and the 
beginning of plastic deformation is shown even more directly (with hardly 


Fig. 4 
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Change in flow-stress accompanying a change in temperature between 90°K 
and 293°k as a function of the flow-stress at 293°K for copper single 
crystals. Full circles were obtained in tension alone, empty marks after 
one or two intermediate twists. The hatched area indicates which 
experimental points would have been obtained if the increase in flow 
stress after a twist corresponding to a surface shear y=0-02 had been due 
to an increase in tg only. 


any theory intervening) in fig. 4, which gives the difference of the flow- 
Stress Tg) at 7';=—90°K and To, at 7',—293°K as measured by change-in- 
temperature tests, plotted against 7.5,. The data corresponding to the 
end of stage II agree with those derived from the experimental results 
of Adams and Cottrell (1955). This is rather satisfactory since these 
measurements relate mainly to stage III. For small strains we were 
unable to confirm the assertion of Adams and Cottrell (1955) that the 
ratio of the temperature dependent part of the flow stress to the total flow 
stress is independent of the prestrain. 
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From the results of the tests with change of temperature we conclude 
that the increased rate of hardening in stage IT is not due to a substantial 
increase in the density of the dislocation forest. This is not surprising 
since the observations on the change in crystal orientation during slip 
(Taylor and Elam 1925, v. Géler and Sachs 1929, Sachs and Weerts 1930) 
show that only a rather limited amount of slip in secondary glide systems 
occurs before the onset of double-slip. If slip in secondary glide systems is 
at all responsible for the increased hardening of stage II it must cause 
(indirectly through the formation of obstacles to dislocation movement 
which cannot be overcome with the aid of thermal energy) an increase in 
7. Whether this is so or not can be investigated by measuring the 
temperature dependence of flow-stress after twisting. 

Figure 4 contains also those results that were obtained in tension after 
twisting. They do not differ from those obtained without twisting. 
This implies that the division of the total flow-stress into the contributions 
by 7, and 7, remains practically unchanged by an intermediate twist, 
i.e. that the slip on secondary systems results mainly in an increase of r,. 
This holds even for large twists, since the data at 7,);,=3-2 kg/mm? 
(indicated by a square) was obtained after a twist corresponding to a 
(surface) shear of y=0-1. The other data were obtained with twists 
corresponding to shears between y=0-01 and y=0-02. The shaded area 
in fig. 4 indicates which experimental values would have been obtained 
if the observed increase of flow-stress by twists of y=0-02 were due to an 
increase in 7, only. . 

We may summarize the conclusions drawn in §§ 2.2 and 2.3 as follows : 
The rapid hardening in stage II and the increase of flow-stress after inter- 
mediate twists are both due mainly to an increase of t,. This is explained 
by the capacity of slip on secondary glide systems to form obstacles which 
hinder slip in the primary slip systems rather effectively. (The nature of 
these obstacles will be discussed in § 4.2.) The results presented in §2.2 
show that the formation of these obstacles must occur continually during 
stage II (and not only at the beginning of this stage). We expect there- 
fore a continual reduction of the slip distance during stage IT. This can 
be checked experimentally by surface observations. 


2.4 Length of Slip-Lines 

The discussion of this section will be based on the assumption that the 
surface markings of deformed face-centred. cubic metals correlate with. the 
slip and hardening mechanism in the bulk of the material. This seems 
to be well justified in view of the results of other investigations in which 
such a correlation could be established (Diehl ef al. 1955, Miler and 
Leibfried 1955, Seeger 1956 a). The following observations were obtained 
by dark field illumination with an ordinary microscope. Comparison 
with electron microscopy work showed that the ‘slip-lines visible 
under these conditions in stage II were arrangements of individual lines, 
the length of which is about half that visible in the optical microscope. 
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In order to observe the variation of slip-distance with strain we have to 
study what we call the ‘ active length of slip-lines ’/. This is the average 
length of slip-lines formed at acertain strain. In the present case where the 
active length of slip-lines decreases with increasing strain it is rather 
different from the average of the length of all slip-lines that were formed 
during the prestrain. For this reason we adopted the method of Blewitt 
et al. (1955) (in its simplest form first employed by Crussard 1945). It 
consists of polishing electrolytically the crystal after various prestrains «, 
thereby removing all traces of the previous deformation, and observing 
the surface pattern after a strain increment 4e. Figures 5 (a)-(f) are 
photographs obtained by this technique from the same crystal (yy»=49° ; 
\y=52° in the notation of Schmid and Boas 1935) after various amounts 
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Reciprocal of active length of slip lines / of copper crystals as a function of 
(shear—) prestrain in tensile tests at room temperature. Additional shear 
strains after polishing de=0-05. The orientations of the tensile axis 
of crystals A, #, and F are as indicated in the stereographic triangle. 


of prestrain in tension. Figure 5 (a) shows the very long slip-lines of 
stage I (prestrain «=0). Figures 5 (b) and 5 (c) were taken in stage IT 
with prestrains «=0-1 and «=0-2. Figures 5(d), 5 (e) and 5(f) were 
obtained in stage HI with prestrains of -=0-3 and «=0-6. They show 
the phenomenon of * fragmentation ’ (slip-lines appear to be interrupted, 
the fragments being displaced sideways), to which we shall return later 
(§ 3.3). 

It can be seen from these photographs that the active length of slip- 
lines does indeed decrease with increasing prestrain. Following Blewitt 
et al. (1955) we have plotted in fig. 6 1/1 (as observed on the top surface of 
the crystal after 4e—0-05) versus shear strain. In stage IT this results 
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in straight lines, the slope of which does not depend markedly on erystal 
orientation. At larger prestrains we note a somewhat stronger decrease 
of the active length of slip-lines with strain, the details depending on the 
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Work-hardening curves of the crystals 4, H, and F shown in fig. 6. 
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Reciprocal of active length of slip lines of a function of prestrain obtained with 
an additional strain of 4e—0-05 after one twist (B and C) or with two 
twists (D). The stress-strain curves are given in fig. 2. 


orientation of the crystals. A comparison with the stress-strain curves 
of these crystals (fig. 7) suggests that the deviation from the straight line 
starts at the beginning of stage TiI. In stage III fragmentation of slip 
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lines occurs (see § 3.3). In stage III we have taken / equal to the length 
of the individual fragments (compare figs. 5 (d) and 5 (e)), following Blewitt 
et al. (1955). The experimental results given by these authors resemble 
those of crystal (A), although their crystal had a somewhat different 
crystallographic orientation (indicated by BCR in the stereographic 
triangle in fig. 6). The rise of 1/1 with strain is not very marked for these 
orientations. This may explain why Blewitt e¢ al. (1955) described the 
relationship between 1/] and « as being linear over the whole range of 
stages IT and III. 

Figure 8 contains the results that were obtained after twisting and 
refers to the same crystals as the stress-strain curves in fig. 2. It will 
be seen that the effect of intermediate twists on 1// is analogous to that 
on the work-hardening curve. A single twist in stage II results in a shift 
of the graph to larger 1/l values; repeated twists however (curve D) 
give an increased slope of the 1// versus « plott. We may therefore say 
that the observation of the surface markings have confirmed the con- 
clusions drawn above on the basis of mechanical measurements. 


§ 3. DynamicaL REcOvERY AND WorK-SOFTENING IN Stace III 
3.1. Survey 


As was clearly shown by the results of Blewitt et al. (1955) on the varia- 
tion of the stress-strain curve of copper single crystals on temperature, and 
as may also be derived from the data of Andrade and Henderson (1951) 
on other face-centred cubic metals if proper allowance is made for the 
dependence of the shear stress-shear strain curve on orientation, the 
slightly temperature-dependent stage II is followed by stage III of the 
stress-strain curve which depends strongly on temperature. Increasing 
the temperature shifts the shear stress 7, (fig. 1), at which stage III 
replaces stage II, to smaller values. A straightforward interpretation 
is that in stage II the dislocated structure of face-centred cubic metals 
varies with plastic deformation but is independent of temperature. This 
has the consequence that 7, is also independent of temperature (if we 
neglect the temperature variation of theelastic constants). A temperature 
dependence of the slope @ of the stress-strain curve may arise through a 
variation of the density of the dislocation forest with strain (Seeger 1954 a, 
1954 b). 

The very small temperature variation of @ in stage II in copper crystals 
is in full agreement with the experimental result (see § 2.3) that the increase 


} For «20-2 this curve was obtained by polishing after the prestrain indicated 
in fig. 8, twisting by y—0-006, straining by 4e~0-02 in tension, twisting by 
y=0-006, straining again by 4e~0-02 in tension and observing the length of 
the slip lines in dark field illumination. By a suitable choice of the conditions 
of illumination it was possible to eliminate the effect of the slip-lines formed in 
torsion thereby retaining the slip-lines formed in tension only. For further 
details see Rebstock (1956), 
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in density of the dislocation forest is rather small y. Since at the transi- 
tion to stage III no abrupt change in the division of flow-stress into T@ 
and r, is observed, we conclude that the temperature dependence of stage 
IIT is due to the change in the dislocated structure with strain (and thereby 
also dr,/de) being temperature dependent. Since unloaded copper 
erystals do not show mechanical recovery at room temperature, this 
change in the dislocation pattern during plastic deformation must be 
due to the combined effects of temperature and applied stress. Such a 
change is of course irreversible and is thus distinguishable from the 
reversible change of 7, accompanying a change in temperature. The 
irreversible changes manifest themselves clearly in the phenomenon of 
‘ work-softening ’, studied in some detail on aluminium crystals by Cottrell 
and Stokes (1955). It seems to be generally agreed that the mechanism 
responsible for the work-softening phenomenon is the same as that causing 
the temperature dependence of 7, in stage III of the work-hardening 
curve. There are mainly two proposals for this mechanism : 

1, Diehl e¢ al. (1955) have proposed and seconded by observations on 
the formation of slip bands in face-centred cubic metals that stage [II is 
characterized by screw-dislocations surmounting obstacles in the primary 
glide system by cross slip. 

2. Friedel (1955) and Cottrell and Stokes (1955) have proposed that the 
temperature dependence is caused by dislocations breaking through 
obstacles on their glide-planes. 

The main object of this paragraph is to decide experimentally which of 
these proposals is the correct one. We shall show, mainly by use of the 
electron microscope, that the empirical results are in agreement with the 
cross slip hypothesis but not with the breakdown hypothesis. 


3.2. Formation of Slip Bands 

Before attacking the questions outlined in § 3.1 we shall briefly discuss 
the formation of slip bands. As was first shown by Heidenreich and 
Shockley (1948) and later studied in detail by other authors (Brown 1952, 
Wilsdorf and Kuhlmann-Wilsdorf 1952a,b) the strong ‘ slip-lines ’ 
visible on the surface of polished aluminium crystals deformed at room 
temperature can be resolved by the electron microscope into a cluster of 
individual slip-lines. Henceforth we shall refer to these individual lines, 
whether clustered into slip bands or occurring as isolated lines, as slip 
lines. 

For some time it was believed that these broad slip bands were typical 
for deformed face-centred cubic metals. Diehl et al. (1955) were able to 
demonstrate that they are typical for stage HI only. They showed that 
enema eg ee 


+ There may be metals (aluminium as a fa-e-centred cubic metal eee 
stacking-fault energy being possibly one of them) where beonuee of t - a 
latent hardening (Seeger 1956 c) the dislocation forest contributes signi bee y 
to the increase of the flow-stress with prestrain. This would ee rise to a 
temperature dependence of the work-hardening coefficient in stage LI. 
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aluminium crystals deformed at the temperature of liquid oxygen, and 
copper and nickel crystals deformed at or below room temperature do 
not exhibit clustered slip bands but a pattern of rather evenly spaced 
individual slip-lines of varying depth, provided the plastic deformation 
was not too extensive (fig. 9 and fig. 10). By polishing and strain- 
increment experiments slip band formation was shown to be a stage III 
phenomenon (fig. 11 and fig. 12). 

The correlation of the first appearance of slip bands with increasing 
strain and the beginning of stage II for various metals and different 
temperatures was so close in these experiments that it was very natural to 
conclude that the formation of slip bands is caused by the same process 
that is responsible for the temperature dependence of 7, as discussed 
above (Diehl et al. 1955). Cross slip had indeed been proposed earlier 
as a possible basic mechanism for slip band formation (Koehler 1952 a, b, 
Leibfried and Haasen 1954, Seeger 1955 b). 

The reason why on aluminium crystals deformed at room temperature 
(or above) slip bands form at rather early stages of deformation is as 
follows: At a given temperature 7,/G depends on the stacking-fault 
energy of a metal (see Seeger 1955 c) in such a way that it is the smaller the 
higher the stacking-fault energy of the metal under consideration 7. The 
stacking-fault energy of aluminium is so high that 7, at room temperature 
is of the same order of magnitude as 7,, stage Il of the work-hardening 
curve thereby reducing virtually to an inflexion point in the stress—strain 
curve. Aluminium crystals show therefore typical stage III behaviour 
already after rather small deformations at room temperature. 


3.3. Fragmentation of Slip Bands 


Blewitt et al. (1955) have demonstrated that the slip markings visible 
on copper crystals in dark field illumination in an optical microscope 
extensively deformed at room temperature show ‘fragmentation ’ (see 
e.g. figs. 5 (d) and 5 (e)). By a detailed investigation of this phenomenon 
we found that extensive fragmentation occurred only if the deformation 
was carried right into stage HII. In stage II no fragmentation was 
observed. The beginning of the fragmentation appears to be connected 
with the rise of the 1// versus « graph (compare fig. 6). 

Since the beginning of fragmentation seems to be related to the begin- 
ning of stage IIT { it was thought worth while to investigate fragmented 
slip bands by the electron microscope. The electron microscope employed 
in the experimental part of this paper was an Elmiskop I (manufactured 
by Siemens Berlin, maximum voltage 100 kv). Three typical photographs 
are shown in figs. 18, 14, and 15. They were obtained from a copper 


} The stress-strain curve for Al at —185°c in fig. 2 of (Diehl e¢ al. 1955) 
has been misdrawn ; in reality stage IIT begins at 7,=1-4 kg/mm?. 

+ This is in line with the observation of Blewitt et ai. (1955) that copper 
crystals deformed at 4:2°x show neither stage ITI nor fragmentation. 
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crystal that was polished after a prestrain of <-~0-6 at room temperature, 
strained in addition by an increment in shear strain of Ac~0-05, and 
treated with the same replica technique as used by Diehl et al. (1955). 
The replicas were shadowed with palladium in such a way that the pro- 
jection of the shadowing direction on the surface of the replica formed a 
rather acute angle with the traces of the primary glide system. 

Figure 13 shows that fragmentation patterns may indeed contain traces 
of one of the secondary slip-planes. In this case as well as in all similar 
ones these traces coincided with the trace of the cross slip {111} —plane 
as predicted from the crystallographic orientation of the crystal. Figure 14 
shows that between the adjacent ends of two fragments traces of both 
the primary and the cross-slip plane may occur in a kind of ‘ staircase 
pattern’. Figure 15 was included to demonstrate that sometimes 
only the primary slip-plane can be observed. Comparing this set of 
photographs there can be little doubt that these traces are also connected 
by cross slip, the cross slip lines however being so weak that they could 
not be made visible by our technique. A general idea of the strength of 
the cross slip lines between two fragments may be obtained by observing 
that the method of shadowing was such as to favour the cross slip traces. 
We conclude from this that the individual cross slip lines are rather weak 
as compared with the slip lines of the primary slip-system contained in the 
slip bands. Sometimes one gets the impression that the cross slip is 
distributed over a certain region rather than concentrated into individual 
lines (see e.g. fig. 16). 

Comparing the fragmentation of the type shown in fig. 15 with the 
clustering of slip-lines into bands, one gets the impression that there is an 
intimate relation between them. The mechanism by which they are 
formed could well be the same. The electron microscope techniques at 
present available give little hope that cross slip between the slip lines 
of a slip band can be made visible. Since however slip bands and frag- 
ments are of closely related appearance and since both of them begin to 
occur at the beginning of stage III we feel justified in concluding that 
the formation of slip bands is due to cross slip. This conclusion had been 
reached earlier on the basis of somewhat more indirect arguments (Diehl! 
et al. 1955, Seeger 1955 b). 

The general impression from both ordinary microscope and electron 
microscope photographs is that in stage III the dislocations overcome the 
obstacles in the primary glide-plane not by breaking through them as 
supposed by Friedel (1955) and Cottrell and Stokes (1955) (see the 
discussion in § 3.1) but by surmounting and circumventing them. This 
argument can be made quantitative on the basis of the results presented 
in fig. 6. As already mentioned the active length of slip lines /in stage III 
was taken equal to the length of those portions of slip bands that follow 
a unique trace of a slip plane, counting the fragments separately. On the 
basis of the breakdown mechanism / should start to increase at the begin- 
ning of stage III. The experimental results show that contrary to this 
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the onset of stage III is marked by an accelerated decrease of |. Even if 
it is admitted that the possibility of some of the weaker slip-lines invisible 
in stage II becoming visible through their sideways growth into slip bands 
might contribute to this acceleration it may be said that this observation 
contradicts flatly the breakdown hypothesis. 


3.4. Work-Softening 

It was already recognized by Cahn (1951) on the basis of investigations 
with the optical microscope that cross slip may be a prominent process in 
the deformation of aluminium single crystals at elevated temperatures. 
The surface patterns called cross slip by Cahn can be made very pro- 
nounced even at normal temperatures by work-softening of aluminium 
crystals (Cottrell and Stokes 1955, Kelly 1956). The work-softening 
procedure consists in prestraining a crystal at a lower temperature, until 


Fig. 17 


E 


Stress-strain curves of a face-centred cubic metal crystal showing the phenome- 
non of work-softening when changing the temperature of the tensile 
test from 7’, to a higher temperature 7',. The dotted line gives the 
stress-strain curve at 7’,. (Schematic.) 


a flow-stress 7 is reached which at a higher temperature would correspond 
tostage III. Ifthe crystal is unloaded, warmed up to this higher tempera- 
ture and strained again a stress-strain curve of the type shown in fig. 17 
is observed. The explanation is that the work-hardened state reached at 
the lower temperature is unstable under the combined action of the higher 
temperature and the applied stress 7. We have already mentioned in 
§ 3.1 that the mechanism by which the unstable part of +, is reduced in 
work-softening is thought to be the same as that responsible for the 
decrease of the rate of hardening in stage III. 

Figures 18 and 19, which were obtained on aluminium crystals pre- 
strained at the temperature 7',—90°K as shown schematically in fig. 17, 
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polished, and further strained at room temperature (7',), give the typical 
surface appearance of a work-softened crystal. At small magnifications 
the slip-bands appear wavy in the ordinary microscope (fig. 18). Larger 
magnification (fig. 19) shows that the slip-bands consist of straight pieces 
connected with what has been termed ‘ prominent cross slip’ by Cahn 
(1951)t. 

It must be stressed that this type of cross slip visible under the optical 
microscope only rarely follows the trace of the crossslip plane. Sometimes 
it doesnot appear to follow at all the trace of acrystallographic plane but to 
be of a non-crystallographic nature. In some cases Cahn was able to 
resolve such non-crystallographic cross slip into a zig-zag pattern by using 
high magnification of the ordinary microscope. 

Definite results on these problems which can be used to discriminate 
between various theories of work-softening can be obtained by the electron- 
microscope. Figure 20 shows an electron micrograph of ‘non- 
crystallographic ’ cross slip. It will be seen that it is resolved into a 
“ stair-case ’ of slip on the primary slip plane and on one of the secondary 
planes. The trace of the secondary slip plane was identified in this (and 
in all other cases investigated) as that of the cross slip plane. This is 
strikingly demonstrated in fig. 21, which is an electron micrograph of 
an aluminium crystal that showed double slip during work-softening. 
Slip bands were formed both by the primary slip system (dark traces) 
and by the conjugate slip system (white traces). The primary slip 
system and the conjugate slip system have their cross slip planes in 
common. The glide directions of the cross slip systems are different 
however. This is borne out in fig. 21. The cross slip traces of both types 
of slip lines are parallel. Since the glide directions in these traces and 
therefore also the nature of the surface steps are different they appear in 
different colours due to the shadowing technique. In the cross slip 
pattern of aluminium crystals work-softened at room temperature and 
not suitably oriented for double slip we never observed more than two 
different slip-planes, corresponding to the primary slip plane and the 
cross slip plane. 

Cross slip can be made visible in the optical microscope also on copper 
crystals, as shown in figs. 22 (a), (b). They were obtained by prestraining 
«=0-6 at the temperature of liquid oxygen, polishing, and work-softening 
at 300°c with an additional average strain of 4e=0-05. As was to be 
expected the cross slip is less marked than on aluminium crystals work- 
softened at room temperature. This is confirmed by the electron micro- 
graph fig. 23, which shows an appearance intermediate between the 
fragmentation of copper crystals deformed at room temperature (compare 


+ We think that the so-called * intimate cross slip ’ observed by Cahn (1951) 
on aluminium is essentially the same phenomenon as the fragmentation dis- 
cussed in §3.3. It was resolved by Wilsdorf and Kuhlmann-Wilsdorf (1952 b) 
with the electron microscope into individual lines of cross slip, which were also 
rather weak compared with the slip lines clustered into bands, 
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fig. 13 to 15) and aluminium crystals work-softened at room temperature 
(compare figs. 20 and 21). 

We were unable to find on work-softened face-centred cubic metal 
crystals those very long slip bands (or slip-lines) that can be observed by 
the optical microscope on plastically deformed single crystals of zinc and 
that would be expected if work-softening were due to a break down of 
barriers in the primary glide plane. The surface pattern was always 
determined by a wavy or fragmented appearance of slip bands which 
with the aid of the electron microscope could be resolved into a sequence 
of slip in the primary slip system and in the cross slip system. This 
shows that under the conditions of work-softening (i.e. high temperature 
and high stresses) it is rather easy for the dislocations to change over 
from their primary slip plane to their cross slip plane and back again to the 
primary plane. The observations on work-softened crystals therefore 
fit very well into the picture which we have drawn for the processes 
responsible for the dynamical recovery in stage ITI. 


Fig. 24 


Rectangular dislocation loops. L, and L, are slip distances of edge and screw 
dislocations. 


§ 4. THEORETICAL Discusston 
4.1. General Discussion of Work-Hardening 

The discussion of our experimental findings as given in §2 and §3 
was of a somewhat phaenomenological nature. For example, we have 
been talking of the effect of glide in secondary glide-systems on the propa- 
gation of slip in the primary glide system ete. Our present task is to 
translate these results into the language of dislocation theory and to fit 
them into a unified picture of the dislocation processes during the plastic 
deformation of face-centred cubic metal crystals. 
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; We have seen that the flow-stress 7(¢) of a metal crystal can be divided 
into two contributions ry and 7, according to eqn. (1). In general 7, 
is inversely proportional to the mean distance between dislocation lines 
penetrating the primary slip plane and therefore proportional to the square- 
root of the density of these dislocations (Seeger 1954 a, 1955 a). Tq is 
somewhat more difficult to interpret quantitatively. It may be said 
however that it is the larger the smaller the mean distance between 
dislocations in the primary glide system is. The coefficient of work- 
hardening due to 7g, which will be denoted by 0,=dr,/de, is related to 
the mean slip distance L of dislocations in the primary glide system. The 
strain increment is proportional to the product of the density of gliding 
dislocations and the slip-distance. The increase in flow-stress however is 
independent of the distance L slipped by the dislocations before they get 
stuck in the crystal. 06, is therefore small for large L and vice versa. 

In § 1 we have briefly mentioned the characteristic features of stage I. 
Both the very long and rather weak slip-lines and the small coefficient 
of work-hardening (in typical cases of pure metals 6/G being of the order 
2x10 to 10-*) are readily explained by assuming that the mean slip- 
distance of the dislocations in stage I is rather large (of the order of 
0:01 cm to 0-1 cm). 

To simplify the following discussion of slip-distances and work-hardening 
we consider the dislocation rings in the glide plane to have the shape of 
rectangular loops. The slip-distances of edges and screws in such a loop 
are denoted by LZ, and L,. The total number of such loops per unit 
volume is called NV. 

The observations reported in figs. 6 and 8 give the active length of 
slip lines / on the top surface of the crystal. We put them equal to twice 
the mean-slip distance L, of screw dislocations. Accordingly the active 
length of slip lines on the side surface will be identified with twice the 
slip-distance L, of the edge dislocations. The procedure of correlating the 
length of slip lines as observed on the surface of the crystal with the slip- 
distances of the dislocations in the bulk of the material is justified by the 
consistency of the results we shall obtain on the basis of this assumption 
presently. 

Our observations of J gave 1=0-20 cm at the beginning of stage II 
(fig. 6). This is in good agreement with the discussion of stage I given 
above. During stage II //2 (for which we now write L,) decreases with 
increasing strain according to the law 


Pig As (eazet se pe ak nee ee) 


e* is slightly smaller than «, and shows a similar dependence on crystal 
orientation as this quantity. A typical numerical value for a crystal in 
the central part of the orientation triangle is A,=4 x 10~¢ cm. 

The available observations on the length of slip lines on the side surface 
are less detailed than those on the top surface, since the first mentioned 
type of slip-lines is more difficult to observe. It can however be said that 
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they were always about twice as long as those on the top surface. We 
can therefore write 

Le=A eae") ae ee 
with A,~2A,. 

In copper crystals deformed at room temperature, stage II of the 
work-hardening curve can be well represented by a straight line (Blewitt 
et al. 1955, Diehl 1956 a). According to the discussion in § 2.3, 6 is not 
very different from 6, in this case. It seems to be a reasonable extra- 
polation from experiment to consider 6, to be constant in stage II. 
The constancy of 6g and the validity of the empirical law eqns. (2) and (3) 
can be related and understood with the following model of work- 
hardening (Seeger 1956 a). 

6, is considered to be due to stress fields around statistically distributed 
groups of a small number of dislocations. The amplitude of the stress 
field generated by the edges of the dislocation loops can be written as 

TOG (IN Lig) lee) & aes es tere a ( E 
NL, is the number of edge-dislocations of one sign per unit area and is 
therefore inversely proportional to the square of the mean distance 
between these edge dislocations. «, takes care of the numerical factors 
involved. It depends on the details of the geometrical arrangement of 
the dislocations, in particular on the effective number m of disloca- 
tions in the dislocation groups. This quantity is characterized by saying 
that the far reaching stress field around such a group is that of a single 
dislocation of strength nb. If tg is due to a distribution of individual 
dislocation lines (n=1) which is probably a good approximation to the 
situation in undeformed crystal (see Seeger 1954 a, c), «, should be of the 
order 1/5. In the present case we expect «,=1/5 n1/?. The infinitesimal 
increase de of the shear-strain due to dN dislocation loops per unit volume 
moving through the slip-distances L, and L, is given by 


de= beled: «.) 260) ocean 


When these loops get stuck in the crystal they will contribute to the stress 
field in the crystal] and therefore cause a change in tg. The work- 
hardening coefficient of this change is 


a Gb (dN dL, 
age + ane... 


If we eliminate ‘ and N from eqn. (6) with the aid of eqns. (2), (4), and 
(5), we obtain the differential equation 


dtq _ %°Gb Sac) en T@ 7 
de 2/1 ee Gs fee eee 85” Si tae (7) 
The general solution of this equation is 
C 
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Since the stress at the beginning of stage II must be finite we have to put 
the constant of integration O equal to zero. The physically admissable 


solution of eqn. (8) is a linear stress-strain curve with the coefficient of 
work-hardening 


PY b \12 
doa. G. (57) ; eg teed HG) 


The same type of treatment holds for screw-dislocations with suffixes 
1 and 2 interchanged. 

The observed coefficient of work-hardening of face-centred cubic metals 
is in stage II of the order 6g=4x10-G. Inserting the A,-value given 
above we find from eqn. (9) 

UGe= ey. SKC LO Pere Bais Aa be iw (10) 


We conclude that «, is of the order 1 and n therefore of the order 20 
to 257. The implication that n does not vary very much with strain in 
stage II is in line with the result of Blewitt e al. (1955) on the electrical 
resistivity due to dislocations as a function of flow-stress. The electrical 
resistivity is under very general assumptions proportional to the number 
of dislocations per unit areat. The experimental result is that it is 
proportional to the square of the flow-stress both in stage II and in stage 
Ill. This shows directly that eqn. (4) with «, fixed is applicable not only 
in stage II but also in stage III. A more complicated expression than 
eqn. (5) will have to be used for the strain-increment in stage III in order 
to allow for the contribution of the glide in the cross slip system. 


4.2. Formation of Obstacles to Slip in the Primary Glide-Plane : 
Lomer—Cottrell Dislocations 


From studies of the change of orientation with strain it is known that slip 
on secondary systems contributes very little to the glide in face-centred 
cubic crystals the orientation of which is in the interior of the stereo- 
graphic triangle. In stage II the number of dislocations in secondary 
glide-systems is therefore in general small compared with the total number. 
The large effect that those dislocations have on 6, must therefore be due 
to a ‘catalytic’ action. We have already (§ 2.3) pictured this as the 
formation of obstacles in the primary glide-plane which cannot be over- 
come with the help of thermal energy under the stress conditions of stage IT. 
Such obstacles would decrease the slip-distances L, and L, and thereby 
increase the coefficient of work-hardening. The only type of obstacle in 
face-centred cubic metals known at present which fulfils these require- 
ments is a special kind of sessile dislocation, the so called Lomer—Cottrell 
dislocation (Lomer 1951, Cottrell 1952). Lomer—Cottrell dislocations were 


This numerical estimate is not too reliable since the ratio between the active 
length of slip lines and the mean slip-distance may be somewhat different from 


two. ; > . 
t For recent discussions of the effect of dislocations on the electrical resistivity 


of monovalent metals see Seeger and Stehle (1956) and Seeger (1956 b). 
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used in the discussion of work-hardening by Mott as early as 1952 (Mott 
1952) and later by Leibfried and Haasen (1954). They play also an 
essential rdle in the views on work-hardening and work-softening held by 
Cottrell and Stokes (1955) and Friedel (1955). 

Lomer—Cottrell dislocations are formed along the intersection of two 
{111}-planes and lie therefore along (110)-directions. The detailed 
arrangement of partial dislocations and stacking-faults is as shown in 
fig. 25. (This is the most important type with which we shall solely 
be concerned. Further types which are less ‘ stable’ are discussed by 
Haasen (1953) and Friedel (1955).) The essential feature is that the 
stacking-fault ribbon ‘ bends over’ from one octahedral plane to another. 
The combination of partial dislocations is confined to both of these planes 
and can therefore neither glide nor climb away from the line of intersection 
of these planes. Since these sessile dislocations can be removed only 
with a considerable expenditure of energy they form effective obstacles 
to the glide motion of dislocations meeting them. 


Fig. 25 
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Lomer—Cottrell sessile dislocation, consisting of three partial edge dislocations 
connected by stacking fault ribbons. 


The Lomer—Cottrell dislocation shown in fig. 25 can be formed in two 
ways, namely either by the reaction of a dislocation of type } 1/2 [01T], 
(111) with one of type 1/2 [110], (111), or by a reaction between 1/2 [110], 
(111) and 1/2 [011], (111). The dislocations of such a pair attract each 
other through their long-range stress fields. The Lomer—Cottrell disloca- 
tion formed as the result of the dislocation reaction between the dislocations 


t+ We denote here dislocations by giving their Burgers vectors and their 
glide planes. They are understood to run parallel to the Lomer-Cottrell] 
dislocation. 
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of each of these pairs is ‘stable’ because energy is required to decompose 
It into the original dislocations and to separate these from each iy 
The following discussion will be based on the convention that ie 
primary glide system is 1/2 [101], (111). This is the octahedral system of 
maximum shear stress if the tensile direction lies in the interior of the 
triangle cross hatched with two perpendicular sets of lines (fig. 26), 


Fig. 26 
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Combinations of glide-systems capable of yielding Lomer-—Cottrell sessile 
dislocations. The direction of the Lomer—Cottrell dislocation is indicated 
by the direction of hatching. 


The range of orientations in which the shear stress is a maximum in the 
other glide-systems with glide plane (111) are also cross hatched. The 
dislocations of the glide-systems with glide-plane (111) can form Lomer— 
Cottrell dislocations with those glide-systems whose triangles of maximum 
shear stress are hatched by the same set of lines}. The direction of the 
Lomer-—Cottrell dislocation is uniquely related to the direction of hatching 
and is the same for full and for dotted lines as indicated in fig. 26. 


+ There are 12 octahedral glide systems and 24 triangles in the graph of fig. 26. 
Therefore two triangles of maximum shear stress belong two each of these 
glide-systems. We have hatched both triangles of such a pair only in the case 
where both are neighbours to the triangle belonging to the primary glide system. 
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Figure 26 is useful for a qualitative discussion of the relative resolved 
shear stress in different glide-systems. The smaller the distance between 
the’ point representing the tensile axis and one of the triangles belonging 
to a certain glide system, the larger is the resolved shear stress in this 
glide-system. For quantitative graphical representations see Diehl 
et al. (1954). 

For the generation of an appreciable number of Lomer—Cottrell dis- 
locations of a certain crystallographic direction it is prerequisite to have 
enough mobile dislocations of the types cooperating in their formation. 
This will be the case if the ‘ effective ’ resolved shear stress (applied shear 
stress minus the latent hardening in the slip system concerned) acting on 
these dislocations is larger than the critical shear stress. Since not very 
much is known about the latent hardening in stage I the following dis- 
cussion is qualitative only and does not allow for the differences in the 
latent hardenings of different glide-systems. 

Extending the idea first proposed by Réhm and Diehl (1952) we think 
that the transition from stage I to stage II begins when the applied stress 
becomes large enough to generate for the first time a considerable number 
of dislocations in secondary glide-systems which are able to form Lomer— 
Cottrell dislocations lying in the main glide-plane. This stress will in 
general, e.g. in copper, depend on temperature, since at least part of these 
dislocations have to cut the forest of dislocations in the primary glide- 
systems which has been built up during stage I. The resolved shear stress 
in the secondary systems corresponding to an appreciable production of 
dislocations increases with decreasing temperature. Since the slope of the 
easy glide region of the stress-strain curve is roughly independent of 
temperature this means that at a given orientation of the tensile axis the 
easy glide region should extend to larger strains «, at lower temperatures 
than at higher ones (Seeger 1954 a). This is exactly what is observed 
experimentally (see fig. 1). 

With the present picture it is possible to understand the variation of 
€, with the crystallographic orientation of the tensile axis, too. The 
orientation dependence of «, in the case of copper single crystals deformed 
at room temperature is given in fig. 27 taken from Diehl (1956 a). Its 
characteristic feature is that the largest values «, occur for orientations in 
the [110]-corner of the orientation triangle but excluding the [110]- 
direction and its immediate neighbourhood. This is contrary to the orien- 
tation dependence of the critical shear stress which takes its minimum 
values in the middle of the stereographic triangle (see fig. 27). Our 
explanation of the variation with crystal orientation is as follows: In the 
whole region adjacent to the plane (011) (which contains the [100]- 
direction and the [1]1]-direction) there is a tendency to form those 
Lomer—Cottrell dislocations characterized in fig. 26 by the hatching with 
vertical full lines. They require dislocations in only one of the secondary 
systems, namely the so-called conjugate system. In the immediate 
neighbourhood of the [110]-direction of the tensile axis there is a tendency 
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to form those Lomer—Cottrell dislocations characterized by horizontal 
hatching either with full or with dotted lines by the combination of disloca- 
tion pairs in secondary glide-systems. In whatever direction the orienta- 
tion of the tensile axis moves away from the [110]-pole into the interior 
of the stereographic triangle the tendency for the formation of this sort 
of Lomer—Cottrell dislocations will rapidly diminish since the resolved 
shear stress in at least one of the two glide-systems of each of the pairs 
will decrease. This accounts for the principal features of the experimental 
situation. 


Fig. 27 


Dependence of ¢, max/€, and 79/79 min on the orientation of the tensile axis of 
copper crystals deformed at room temperature. 


The end of the transition between stage I and stage II will presumably 
be reached when a sufficient number of Lomer—Cottrell dislocations in all 
of the three possible directions are formed. From there on the slope of 
the stress-strain curve is mainly determined by the rate of production of 
Lomer—Cottrell dislocations which must be such as to bring about the 
empirical law eqn. (2) for the variation of slip-distance with strain. 
These processes should be rather insensitive to the degree of perfection of 
the undeformed crystal, to its impurity content and to the diameter of 
the specimen. This agrees with the observations that for a given crystal 
orientation the work-hardening coefficient is very well reproducible in 
stage II. Its variation with crystal orientation is small compared with 
that of the length of stage I (Diehl 1956 a). This is understandable since 
so many secondary glide-systems cooperate in determining @ in stage IT 
that a small dependence on the orientation of the tensile axis should 
result. 

We have so far discussed the réle of Lomer—Cottrell dislocations in 
holding up the dislocations moving in the primary glide plane and in 
determining the slip-distance and thereby also the coefficient of work- 
hardening in stage II. As was pointed out by Mott (1952), Lomer-— 
Cottrell dislocations may significantly contribute to the stabilization of the 
deformed state against backward glide after unloading by forming 
obstacles in the rear of groups of dislocations piled up under the action 
of the applied stress. In our present picture this implies that the stability 
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of the deformed state is less in the easy glide region where no Lomer— 
Cottrell dislocations are available than in stage II. We have found a 
striking confirmation of this view in a recent paper by Buckley and 
Entwistle (1956). These authors show in their fig. 8 that the rate of 
increase of the Bauschinger effect (measured in terms of the “ Bauschinger 
glide strain ’ which is a measure for the instability of the deformed state 
under unloading or reversed loading) with increasing hardening is much 
larger in the easy glide region than in stage II and III of aluminium single 
crystals strained at room temperature. 


4.3. The Mechanism of Slip-Band Formation and Fragmentation 


The principal result of the discussion of experiments in § 3 was that the 
cross slip of screw dislocations is the governing process for the tempera- 
ture dependence of the work-hardening in stage III, for the phenomenon 
of work-softening, for the formation of slip-bands, and for the fragmenta- 
tion of slip-bands. The object of the present section is to consider some 
of the theoretical problems connected with these phenomena. 

As pointed out earlier screw dislocations undergo cross slip in order to 
circumvent and surpass obstacles which they are unable to break through. 
The discussion in § 4.2 suggests that these obstacles are Lomer—Cottrell 
dislocations. The Lomer—Cottrell dislocations holding up screw disloca- 
tions are those which cannot be formed by dislocations of the primary 
glide system but require the cooperation of dislocations in the other glide- 
systems with the primary glide plane. It is therefore essential (Seeger 
1956 a) not to restrict the discussion of Lomer—Cottrell dislocations to 
those formed by the dislocations in the primary system glide as was done 
by other authors. 

Next we will take up the related problems of the detailed mechanism of 
cross slip and of the stability of Lomer—Cottrell dislocations against 
‘collapsing’ under the combined action of stress and temperature. 
Dislocations generated during plastic deformation in face-centred cubic 
metals are extended ones. They can undergo cross slip only if the two 
partial dislocations of an extended screw dislocations are brought together 
over a certain length. The resulting complete screw-dislocation may then 
separate into partials in the cross slip plane. 

Lomer-—Cottrell dislocations can give way to the pressure of a group of 
dislocations piled up against them in several ways which have been 
discussed by Stroh (1956). For our present purpose it suffices to consider 
that way with the lowest energy barrier. Its first step is to reunite over a 
certain length the partial dislocations which make up the Lomer—Cottrell 
dislocation. The calculations of the energy barrier (activation energy) 
involved is very similar to the calculation of the energy barrier of cross 
slip and will have to follow the lines indicated in a simple case by Schoeck 
and Seeger (1955), Seeger (1954 d). Since more detailed calculations of 
this type are still under way we will confine ourselves to a qualitative 
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discussion of the relative magnitudes of the activation energies in both 
cases. The same stress component is involved in bringing together the 
partial dislocations of the screw dislocations and of the Lomer—Cottrell 
dislocations against which they are piled up. The lengths over which the 
recombination to complete dislocations has to be achieved depend on 
different components of the stress tensor but will not vary very much 
under typical conditions as shown by the calculations of Schoeck and 
Seeger (1954). The parameter which does determine the activation 
energies is the width of the stacking-fault ribbon between the partial 
dislocations involved. These widths are very much smaller in aluminium 
than they are in copper (or other low stacking fault energy metals) 
(Seeger and Schéck 1953). This implies that the activation energies for 
aluminium are much lower than those for copper under equivalent stress 
conditions. In a given material the width of the stacking-fault varies 
with the character of the complete dislocations. It is a minimum for 
screws and a maximum for edges. For Lomer-Cottrell dislocations it is 
somewhat smaller than for edges but still considerably larger than for 
screw dislocations. This has the consequence that the activation energy 
for cross slip is smaller than that for the collapse of Lomer—Cottrell 
dislocations, and that piled up groups of screw-dislocations will ‘ leak ’ 
by cross slip before they are able to break down the Lomer—Cottrell disloca- 
tions holding them up in the primary glide plane. The theoretical situa- 
tion is therefore in agreement with the facts deduced from experiment. 

The preliminary calculations so far available show that the temperature 
variation of 7. (being the resolved shear stress at which a noticeable rate 
of cross slip is reached) can be accounted for by the theory. The theory 
also explains fully the quantitative difference in the behaviour of 
aluminium and copper which is particularly obvious from the experimental 
results presented by Diehl et al. (1955). The empirical correlation with 
the magnitude of the stacking fault energy, though at that time based on 
much less complete experimental evidence than available now, was in 
fact the starting point for the work on slip bands presented here (Seeger 
1955 a). 

In §§ 3.2 and 3.3 we have presented the empirical evidence that the 
transfer of slip from one slip-line in a slip band to an adjacent one is due 
to cross slip. We shall now discuss possible mechanisms in some detail. 

It follows from the preceding discussion that in stage IIT screw disloca- 
tions are pushed out by cross slip from the head of dislocation groups 
piled up against Lomer-—Cottrell dislocations. The stress components at 
the head of such a pile up are proportional to the number of dislocations 
in the pile up (Eshelby ¢¢ al. 1951). The strongest and longest slip-lines 
will therefore leak first. (In quantitative considerations one must take 
into account, however, that the number of dislocations effective in a pile- 
up will in general be considerably smaller than the number of oman 
rings in the slip-line due to the * screening ; effect of dislocations in the 
environment of the slip-line.) The dislocations in the cross slip plane 
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will at first be pushed away from the group with a repulsive force varying 
roughly inversely proportional to the distance from the group. Later 
they will either get stuck in the work-hardened material, or they will meet 
another piled-up group of screw dislocations of the opposite sign. They 
will then annihilate with some of these, thereby giving the way free for 
more dislocations following the same way. Since the annihilated disloca- 
tions in the group just mentioned will in some cases be replaced from their 
source a weak slip line (say of the elementary structure type) may develop 
into a strong one. This mechanism is presumably stopped by the edges 
of the dislocation rings spreading in the cross slip plane. There are no 
dislocations groups with which they can annihilate ; they will therefore 
get stuck in the work-hardened material after a certain slip distance and 
create a pile-up the back stresses of which will eventually prevent further 
cross slipping. This accounts for the fact that the cross slip lines are 
rather weak as discussed in § 3.3. Since the stress magnification acting 
on dislocations which are not quite at the head of a pile-up is still large, 
these dislocations may undergo cross slip, too, thereby giving rise to more 
of these weak cross slip lines. This explains why one has frequently the 
impression that cross slip is spread out over an area rather than concen- 
trated in strong lines. 

The screw-dislocations which have not annihilated after having moved 
out of the primary glide plane by cross slip may under the action of shear 
stresses and temperature undergo once more cross slip back to the primary 
glide plane. Since by the argument given above the edges remaining in 
the cross slip plane will be held up somewhere the whole configuration 
forms a Frank—Read source in the primary glide plane and may by the 
Frank—Read mechanism generate a new slip line. This is (with modifica- 
tions) the mechanism of “‘ double cross slip ”’ first envisaged by Koehler 
(1952 a,b). It should be mentioned that also a mixture might occur 
of the two mechanisms which we have discussed: the dislocations may 
annihilate after cross slipping back to the primary glide plane. 

Depending on the scale on which these processes happen they may either 
lead to slip band formation or to fragmentation. If the new slip lines 
are formed within a few hundred Angstréms distance from the original 
one we observe clustering of slip lines into bands. If the dislocations slip 
farther away in the cross slip plane before creating a new line or making 
an existing one grow one obtains fragmentation. Both the ‘ double 
cross slip ’—and the ‘ annihilation ’ mechanism allow glide to proceed with 
a smaller increase in work-hardening than would correspond to stage II. 
They are therefore both able to account for the decrease of the work- 
hardening coefficient in stage III. The relative importance of the two 
mechanisms will depend on the conditions of stress and prestrain. If the 
material is very little work-hardened as is aluminium at the beginning of 
stage III at high enough temperatures double cross slip may well be 
important. Ina highly work-hardened material however the annihilation 
mechanism should be predominant. It accounts in particular for the 
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drop in flow-stress occurring during work-softening which in some way 
must be connected with the annihilation of dislocations which had con- 
tributed to the work-hardening. 

We conclude with a few remarks on the conditions at large strains. 
As can be seen from fig. 24 of Diehl e¢ al. (1955) at high enough strains 
there is no slip in the material between slip bands. Slip proceeds by 
adding more lamellae to slip bands (and possibly also by growth of slip 
lines within the bands). The question may arise whether a slip-line at the 
edge of a band may be so weak and contain so few dislocations that none of 
these can be expelled by cross slip, thereby stopping the sideways growth 
of the band. The general impression is that bands continue to grow side- 
ways with increasing strain. This is presumably due to the fact that the 
screw-dislocations in the interior of a band help to push out by cross slip 
the dislocations in the outer slip lines. 


§ 5. ConcLupING REMARKS 


We shall not summarize here the conclusions concerning the various 
dislocation mechanisms in the work-hardening and _ work-softening 
phenomena of face-centred cubic metals reached in this paper. We 
should rather like to stress the unifying réle played in these considerations 
by the stacking-fault energy. Dislocations generated during cold-work 
in face-centred cubic metals are extended ones, and their behaviour under 
the action of stress and temperature depends to a large measure on the 
width of the stacking-fault ribbons between the half dislocations. Until 
it was realized that aluminium, being a high stacking-fault energy metal, 
differs in this respect from other face-centred cubic metals such as copper, 
silver, gold and nickel it was impossible to integrate the seemingly con- 
flicting evidence on the temperature-dependence of flow-stress, on the 
temperature-dependence and the shape of the work-hardening curve, 
on surface markings and slip-band formation into a unified picture valid 
for all face-centred cubic metals. We hope to have been able to develop 
such a picture, which of course still requires refinement by future 
investigations. 
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INTRODUCTION AND SUMMARY 

It is the purpose of this note to show that the fluctuations in the 
instantaneous intensities at two arbitrary points in a stationary optical 
field are in general correlated and that this correlation is proportional 
to the square of the coherence function introduced into the theory of 
partially coherent fields by Blanc-Lapierre and Dumontet (1955) and 
independently by Wolf (1955). The result has a bearing on recent experi- 
ments by Hanbury Brown and Twiss (1956) in which fluctuations in the 
current output from two photoelectric cells situated close to each other 
and illuminated by an almost coherent beam of light were found to be 
correlated. 


§ 1. A THEOREM CONCERNING THE GAUSSIAN RANDOM PROCESS 
THE main result of the paper is an immediate consequence of a theorem 
on the two-dimensional stationary Gaussian random process, which we 
shall now establish. 

Consider two random functions X(t) and Y(t) of time ¢ and assume 
that the probability of finding X in the interval X, X+6X and Y in 
the interval Y, Y+é6Y at any given instant is P(X, Y)6X5Y, where P 
is the two-dimensional Gaussian distribution 


L ee 
= ee OX 4 Sees C(LD) 
icin) 277 \/(H11b22—H12") P| 2(HyiMe2—H 12)” 
Here p,,(i,j=1, 2) are the second moments defined by 
My1=(X"), M22=(Y?), Myp= (XY), OE este (h.2)" 


sharp brackets denoting average value. We need not distinguish 
between a time average and an ensemble average as we are dealing with 


an ergodic ensemble. reg." 
It follows that the average value of the product X?Y? is 


(X?Y2) = (aa X°y2P(X, Y)dXdY 
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To evaluate (1.3) we set 


X=o (=e, y=o,/(=) Romer fC): 


We also set 
Hys=—V(Hiae) CSP; - - - - + (1.6) 
then 
o=+/(tibee) Sing. - - - . + - (17) 
Substitution from (1.5)—(1.7) into (1.3) gives 
4o5 
ec 


x [7 ["* ertexp| —(e-+12+260 008 4) |dedn. (1.8) 


The integral on the right-hand side of (1.8) is well known ; its value is} 


scope a a[1+2 cos? d] 
22 Sez 2 pt) ea ee 
[| @rtexpl— (498+ 2&y 008 $)]dédn= eG 
[MaMa + 212")? 29? 
= 406 
(1.9) 
Hence (1.8) becomes, if we also substitute from (1.2) : 
CXS V2) == CX CYA) aE OC a eo eee) 
We also have the identity 
(XA Y 23 — (XA) CYA ae AX AA YA ek) 
where 
A X*=2X4#—{ X48) AY ta YES CY?) «  eeetict2) 
From (1.10) and (1.11) we finally obtain 
CX A V2 Ne? OXY ee eee (i213) 


The relation (1.13) implies that the fluctuations 4X? and 4 Y? are in 
general correlated and that this correlation is equal to twice the square 
of the correlation between X and Y. 


§ 2. CORRELATION IN INTENSITY FLUCTUATIONS AT Two ARBITRARY 
Ponts IN A STATIONARY OpTicaL FIELD 


We now consider a stationary optical field and fix our attention on 
the disturbance at two points P, and P,. For the sake of simplicity we 
shall regard the disturbance to be a real scalar, e.g. a rectangular com- 
ponent of the electric field vector. 

Let V(t) represent the disturbance at P, at time ¢ and V,(t+7) the 
disturbance at P, at time t+7. Keeping 7 fixed for the moment, we 
may regard V, and V, as the components of a two-dimensional vector 


} See, for example, Rice (1945), p. 69, 
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which is the sum of a large number of independent random vectors 
(contributions from the randomly emitting atomic radiators that form 
the sources of the field). Under these circumstances the conditions of 
the central limit theorem} of the calculus of probability are effectively 
satisfied and consequently the joint probability distribution P(V,, Vs) 
may be assumed to be Gaussian. Hence we may apply the theorem 
(1.13), with X(#)=V,(é), Y(t)=V,(t-+7). If we also set 


LHS ¥7(8), Plt) Vistlt), oT 
AI, (t)=I,(t)— (1,(4)), AI ,(t)=1,(t)— (I(t) me cs 
we obtain 
(AT, (t()AT(t-+-7))=2¢V,(t)Va(t+7))2 2. (2.2) 


In the limiting case when the points P, and P, coincide and at 
the same time 7=0, the discriminant p4,u..—p 2 in (1.1) is zero; 
nevertheless the formula (2.2) retains its validity also in this case, as 
can be seen by considering the probability distribution for X=V, alone. 
One then obtains (dropping the common suffix) 


C4 208 EURO OS eg an eat emer FE) 
The cross-correlation function 
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which occurs on the right-hand side of (2.2) and which is a generalization 
of the mutual intensity function of Zernike (1938) was introduced into 
the theory of partially coherent fields by Blanc—Lapierre and Dumontet 
(1955) and independently by the present writer (Wolf 1955). Morover, it 
was shown in the latter paper that in vacuo, the correlation function 
obeys rigorously the two wave equations 


1 67d 40(7) 
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where V,2 is the Laplacian operator with respect to the coordinates of the 
point P,. Thus although J is a quadratic function of the basic field 
variable (V), it behaves as the amplitude of a wave. By analogy with 
other wave theories one might have therefore suspected that J? should 
also possess a physical significance. Our formula (2.2) shows that J? is, 
in fact, proportional to the correlation between the intensity fluctuations at the 
points P, and P,, the intensity at P, being considered at tume 7 later than 
a Ps. 

To measure J,,(7) one may isolate the vibrations at P, and P,, for 
example by placing an opaque screen across the field, with two small 
openings at these points. From observations on the (time averaged) 
intensity some distance behind this screen, the values of J,.(7) may be 
determined (Wolf 1955, pp. 252-253). In an important recent publica- 
tion, Hanbury Brown and Twiss (1956) proposed a radically new method 
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for the measuring of correlations between light beams. Here photo- 
electric cells are placed at P, and P, and the correlation between the 
fluctuations in the current output from the two cells is determined. 
Hanbury Brown and Twiss (1956, 1957) found (see also Kahn 1957) that 
in certain cases this correlation is closely related to the correlation that 
exists between the amplitudes V, and V,. The principle of this method 
may at once be understood from the formula (2.2) : 

A photocell measures effectively the square of the electric vector, so 
that under ideal experimental conditions the left-hand side of (2.2) and 
hence the amplitude correlation J,,(7) may be determined by such an 
arrangement. In practice, one has not, of course, ideal square law 
detectors and ideal correlators, but instruments with finite resolving 
time (high-frequency cut-off). In consequence one cannot measure the 
correlation between the instantaneous values 4/, and AJ, but rather the 
correlation between certain weighted short-time averages AT and TMi 
Moreover, at very low intensities, the quantum nature of the absorption 
process may play an important part. Nevertheless, the basic principle 
of the method proposed by Hanbury Brown and Twiss is clearly brought 
out by the relation (2.2), which also suggests that this method is of very 
general validity. 
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SUMMARY 


The factors influencing the development of helical dislocations found 
in CaF, by Bontinck and Amelinckx are discussed. At temperatures 
near and below 700°c the principal electrolytic current in CaF, appears 
to be carried by the negative ions. Thus climb of dislocations is difficult 
to produce as a result of processes inducing electrolytic flow. It is 
proposed that the positive ions carry a larger component of the current 
near 1200°c and that maintenance of the electrolytic currents 
induced during additive coloration produce what might be termed mild 
climb. Under this climb dislocations having a strong screw component 
are induced to assume a helical or spiral prismatic form as proposed 
earlier by Seitz. The model of F-centres in CaF, and their equilibrium 
characteristics are discussed. It is pointed out that helical dislocations 
should be generated in many materials under a wide variety of circum- 
stances in which climb occurs and that such dislocations are, in principle 
at least, capable of generating whiskers when they occur sufficiently near 
the surface. Such whiskers would grow from the base and should possess 
many of the attributes found in the whiskers of this type which are 
observed in practice. It is proposed that controlled climb of any material 
may be produced with the use of an appropriate temperature gradient if 
the production or annihilation of lattice defects at the dislocations en- 
genders climb and if the production lies in an attainable temperature 
range. 


§ 1. IyTRODUCTION 


BonTINCK AND DEKEYSER (1956) have demonstrated that dislocation 
lines in calcium fluoride can be decorated by the colloids formed during 
the precipitation of excess metal in additively coloured crystals. A 
fraction of the dislocations decorated after additive coloration at 
higher temperatures, near 1200°c, possess helical character and merit 
particular discussion. The axes of the helical dislocations appear to be 
preferentially oriented in the (110) directions, which presumably are 
also the directions of the Burgers vectors of the dislocations. As pointed 
out by Bontinck and Amelinckx (1957), they seem to be one species of the 
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spiral prismatic dislocations proposed by Seitz (1950 a, b, 1952) in 1950 
which can be formed by the addition or subtraction of molecules (climb) 
from a dislocation which is primarily, or at least strongly, of screw character. 
The spacing between successive turns of the helix is large compared to 
the Burgers vector in the cases observed in this laboratory. The purpose 
of this note is to discuss some of the circumstances surrounding the 
formation of helical dislocations in CaF, and to point out possible attri- 
butes of such dislocations, particularly their ability to form whiskers 
under appropriate conditions. It should be emphasized that although 
whisker growth has actually been observed on calcium fluoride in this 
laboratory, the conditions under which whiskers appear have not yet 
been clearly associated with spiral dislocations. 


§ 2. ADDITIVE COLORATION OF CALCIUM FLUORIDE 


Mollwo (1934) has found that CaF, may be additively coloured by the 
addition of excess calcium under conditions which resemble very closely 
those under which the alkali halides may be coloured. Both electrolysis 
and addition from the vapour or bulk metal are effective. Actually CaF, 
appears to accept excess calcium far more readily than the alkali halides 
accept alkali metals. Mollwo has presented a figure (Mollwo, 1934, fig. 6) 
which indicates that the concentration of excess calcium varies in the 
manner A exp (—E/kT') with temperature. Here as usual, A and E are 
constants, k is Boltzmann’s constant and 7’ is the absolute temperature. 
The experimental values of the concentration of excess metal were deter- 
mined by various investigators in four different ways over four ranges of 
temperature extending from 300°c to the melting point at 1360°c. The 
constant H was found to be nearly the same as the heat of sublimation 
of elemental calcium (about 48 kg cal per mol), however the constant A 
was about 104 times larger than the corresponding constant in the 
expression for the density of atoms in the equilibrium vapour of the metal. 

In any event, the experiments show that deviations from stoichio- 
metry of the order of 1°%, or more can be obtained by additive means at 
appropriate temperatures near 1000°c. The decrease in mass density 
of the crystals during such addition suggests that, when atomically dis- 
persed, the stoichiometric excess is present in the form of fluorine-ion 
vacancies to which electrons are attached. In fact Mollwo used this 
assumption to determine the concentration of excess metal in the cor- 
responding range of temperature. Unfortunately, the lattice parameter 
measurements which should be made in conjunction with density 
measurements to test this model have not yet been carried out. As an 
alternative, it is possible that the excess metal is present interstitially. 
It is evident that the order of the reaction governing solution of the metal in 
the salt would be different in the two cases because of the divalent character 
of calcium and the monovalent character of fluorine. We shall discuss the 
nature of the colour centres further below after considering other evidence. 

When appropriately quenched from high temperatures, the additively 
coloured crystals exhibit two prominent absorption bands, provided the 
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concentration of added metal exceeds that of the miscellaneous impurities 
present, which seem to be of the order of 10 or more parts per million in 
the typically good crystals available at present. One of the bands, 
designated the « band by Mollwo, occurs near 3-3 ev (38504) and appears 
to be the analogue of the F band. The other, termed the B band, is near 
2-3 ev (5200) and seems to be the analogue of the R or M band as in the 
alkali halides, although the information available concerning its properties 
are meagre. Liity (1953) has shown that many of the bands observed by 
Smakula in x-rayed specimens are the result of combination of the « and 
8 centres with the foreign agents present in the best crystals available at 
present. 

Additive coloration has been produced in the crystals studied in this 
laboratory by heating with calcium and by heating either in the presence 
of hydrogen or sodium. Sodium has been used in the manner developed 
by Rexer (1932) and employed by Amelinckx (1956) in the case of sodium 
chloride. Presumably the treatment with hydrogen or with the alkali 
metal produces partial reduction and a relatively moderate density of 
excess metal. The extinction coefficient at the peak of the « band 
attains values of the order of several mm~! in quenched specimens, cor- 
responding to concentrations of centres between 10!” and 1018 per cm?. 
More intense darkening can be obtained with the use of calcium. The 
decoration of dislocations can be observed in specimens darkened at 700°c, 
however the helical dislocations have been observed only in specimens 
darkened at higher temperatures, in the vicinity of 1200°c. 

Other observations on x-ray darkened crystals appear in publications 
of Barile (1952), Smakula (1953, 1954), Schulman ef al. (1952), Urbanck 
(1951) has studied photoconductivity, whereas Deigen (1951) has given 
a theoretical treatment of double centres at high concentrations. 
Feofilov (1953) considered the F-centres to be anisotropic, situated in the 
crystals with a definite orientation. 

CaF, is a moderately good ionic conductor. For example, Smakula 
(1950) has found conductivities of the order of 10°! (ohmem)™? at 
100°c which are comparable to the conductivities of alkali halides of 
normal purity. The question of the sign of the carrier in the halides of 
the divalent metals has been a topic of general discussion for an extended 
period and it now appears that Ure (1956)+ has demonstrated that the 
negative ions carry the bulk of the current in the range from room tem- 
perature to 900°c, the disorder corresponding to Frenkel defects (i.e. 
vacancies and interstitial ions) in the fluorine lattice. This has been 
established by studying the electrical properties of crystals doped with 
NaF and with YF,. The first addition presumably enhances the density 
of F-ion vacancies and the second the density of interstitial F ions. 
Ure has found that the activation energy for migration of the vacancies 
is 20-L1-5 kg eal per mol at 200°c and decreases to about 12-41 kg cal per 
ee 

+ We are indebted to Professor Smakula for providing us with an account of 
this work, 
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mol at 600°c. The activation energy for migration of interstitials is 
38-48 kg cal per mol in the range between 640° and 920°c. He has found 
that the energy of formation of the defects is about 32-6 kg cal per mol in 
the same range of temperature, this being twice the activation energy 
for formation of the defects, presumably by a bimolecular process. 

It may be worth noting at this point that the positions occupied by 
the interstitial fluorine ions in the simple cubic lattice of the normal 
fluorine ions found in CaF, presumably are body-centred ones. Half of 
such positions are occupied by calcium ions. Thus the ideal free volume 
available is equiv alent to that of a calcium ion. The radius of a fluorine 
ion is about 30% larger than that of a calcium atom (about 1-334, com- 
pared to 1-06). Hence the fit is tight, but not unreasonably so. 

Actually, Ure’s results seem to be in excellent accord with the results 
of older work (Hauffe 1955) which may be summarized as follows. First 
the transport numbers of the negative constituents of the halides of 
divalent metals are usually close to unity (Tubandt 1932). PbI, seems 
to provide an exceptional case among the salts studied. Second, the 
conductivity of the salts can be enhanced (Gyulai 1931, Seith 1929) by the 
addition of the alkali halides. This is the case for example when NaCl 
is added to PbCl,. Presumably the alkali ion enters substitutionally and 
is matched by a halogen ion vacancy. This procedure has been used by 
Ure in the work quoted above. ‘Third and probably closer to the problem 
at hand, Zintl and Udgard (1939) have found that when appreciable 
quantities (33-59%) of YF, are added to CaF,, the Y ions enter substi- 
tutionally in the CaF, lattice and the extra F ions are located in inter- 
stitial positions which are probably the normally unoccupied body centred 
positions in the simple cubic lattice of fluorine ions. Zintl and Udgard 
have drawn the same conclusion for the system SrF,—LaF, and the system 
CaF,-ThF,. Ure has used the properties of YF, additions in his work. 

In view of the existing evidence, it seems safe to conclude that the most 
common defects affect only the halogen lattice and are of the Frenkel type. 
It may be noted in passing that the activation energy for migration of the 
carriers in Smakula’s specimens is about 0-9 ev which is close to the value 
Ure obtained for the halogen ion vacancies. It is not possible to state 
whether the carriers are vacancies or interstitial fluorine ions. 

It follows that the transport and structural work supports the view 
that the colour centres formed by adding excess calcium are halogen 
vacancies in which electrons are bound. In the range where the electro- 
lytic conductivity is essentially exclusively the result of the migration of 
halogen ions, the process of colouring additively by diffusion from the 
surface would appear to involve the simultaneous migration of electrons 
and halogen ion vacancies from the surface into the volume and of inter- 
stitial ions to the surface. The flow of halogen ion vacancies should be 
the most important of the two modes of migration of the fluorine ion if 
Ure’s activation energies are valid. The additional calcium ions and the 
halogen ions brought to the surface by the creation of halogen ion vacancies 
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there or by the migration of interstitial ions to the surface form new layers 
of the salt. It seems natural to suggest that this mechanism of additive 
colouring operates exclusively in the lower range of temperature near 700°C, 
Unless some of the calcium atoms diffuse into the crystal along dislocations, 
which emerge at the surface, and form F-centres within the erystal by 
giving up the valence electrons and capturing interstitial fluorine ions, SO as 
to induce climb of the dislocation, there is no good reason for expecting a 
change in the dislocation pattern during the colouring process. Ulti- 
mately, when the crystal is cooled and the F-centres form colloids the 
existing edge dislocations will be decorated. The condensation of 
F-centres may then cause climb of the dislocations, but we shall assume 
this effect is small. 

_ The specimens which are darkened at higher temperatures, near 1200°c, 
display helical dislocations and hence give evidence that climbing does 
occur at the higher temperature. It is possible that such climb is the 
result of the diffusion of calcium atoms along dislocations in the manner 
described above. However it seems much more likely that the lattice of 
the metal ion begins to become disordered at the higher temperature. 
A pair of interstitial fluorine ions may be converted into a calcium ion 
vacancy at a dislocation. Although this process is endothermic, it is 
evident that entropy will be gained because of the disordering of the 
Ca—ion lattice once interstitial halogen ions become abundant. Thus it 
seems entirely reasonable to suppose that a fraction of the lattice defects 
present near 1200°c are positive ion vacancies. 

Once positive ion vacancies become available as carriers, a part of the 
electron current from the surface may be neutralized by a flow of positive 
ion vacancies to the surface. This component of compensating current 
is particularly important because the electrons will have a tendency to 
become bound at the halogen ion vacancies and immobilize them in part. 
New positive and negative ion vacancies will be generated at dislocations 
(and other sources such as voids, if present) to replace those which are 
consumed by diffusion to the surface or by conversion to F-centres. Thus 
the dislocations are induced to climb as a result of the presence of a sub- 
normal concentration of vacancies of both types. Some climb presu- 
mable would occur merely as a consequence of heating and the production 
of the normal equilibrium density of positive ion defects at the higher 
temperatures, but the climb will be enhanced by additive coloration. 
During slow cooling, the positive ion vacancies which do not migrate to 
the surface or become involved in F-centres should return to the sources, 
which then act as sinks. In an uncoloured crystal, one might expect the 
climb of dislocations when they act as sources during heating to be 
balanced by the climb in the opposite direction when they act as sinks, 
although a redistribution of material from one part of the lattice to 
another could be associated with irreversible changes in. the dislocation 
pattern. On the other hand, one should always expect an irveversible 
climb during the process of additive coloration at higher temperatures if 
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+he mechanisnr postulated here in which a part of the electrolytic current 
countering that of the electrons is carried by positive ion vacancies is 
valid. The magnitude of the irreversible climb should depend upon the 
intensity of coloration. 

~ ‘Let us now return to the problem of the origin of the high density of 
colour centres which can be obtained in CaF. We shall assume that they 
consist of halogen ion vacancies which are neutralized by single electrons 
‘and that two such centres are formed for each excess calcium added to 
the crystal. We shall also assume that the cohesive energy of the electron 
to the F-centre is of the order of 1-5 ev, so that the fraction which are 
ionized under equilibrium conditions at 1200°c is not sufficient to affect 
the order of the reaction 


Ca (atom) = 2 -F-centres, 9.5 Fae ee) 


We shall adopt a reference state in which the free energy of a calcium 
atom in the metal is zero at absolute zero of temperature and is 7,,, a 
‘positive quantity, at the energy of formation of F-centres. The chemical 
potential of an atom in the vapour will be taken in the form 


T p= ented LOS (CIC le eos ie ane ee ae 
in which c, is the number of atoms per unit volume of the vapour and C 
is equal to (27 MkT/h?)3 in which M is the mass of the calcium atom, 
k is» Boltzmann’s constant, 7’ is the absolute temperature and /h is 


Planck’s constant. The corresponding chemical potential for a Ca atom 
dissolved in the crystal to form two F-centres will be taken in the form 


r=e,+2kT log (2c,/KNq) ei 2c eee) 


in which c, is the concentration of excess calcium atoms in the crystal of 
Cal, and « is a parameter, presumably of the order of 100 which takes 
account of the additional entropy gained by the crystal over and above 
the mixing entropy NV, is the density of fluorine sites in the crystal. The 
quantities «, and «, are the heats of formation associated with the calcium 
atoms in the vapour and Cal’,, respectively. 

We may determine the density of atoms in the vapour from the relation 
r,=r,, and the partition of calcium atoms between vapour and crystal from 
the condition 7,=r,. If we use these equations to determine the ratio 
c,/C,, we obtain 


Co[Cy=(N 9/20) exp {(2e,—e,—7,, /2kT}. 2 — 2 (4) 
If we accept Mollwo’s relation, according to which this ratio is approxi- 
nately 104 over a wide range of energy, we may conclude that 


eot2e,=96 kgcalpermol . . . . 4... (5) 
and 

(Nox/2C) exp (=n /2h TOS 2) ere) 
The first relation implies that the energy required to dissolve a calcium 
atom in Cal, from the vapour is about 48 kg cal per mol or about 2-1 ev 
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per atom. We shall consider this further below. The relation (6) seems 
to imply that « is considerably larger than 100. The quantity N» is 
8-0 x 10°? per cm’, whereas C is 2:8x 1022 at 1500°x. Since the expo- 
nential coefficient is less than unity, it follows that « should be of the order 
of 7 108 or more. 

The very large value of « found in this way presents somewhat of a 
mystery since it seems doubtful if the contribution to the internal entropy 
of the crystal, exclusive of the mixing term, can be associated with a factor 
so large. In this respect the situation would not be improved, but rather 
made worse, if we were to assume that the electrons which are normally 
bound to the F-centres become free during the process of additive color- 
ation. The partition factor corresponding to C above arising from free 
electrons is smaller by a factor (m*/M)*?, in which m* is the effective 
electron mass and M is the mass of a calcium atom, so that the contri- 
bution to crystal entropy from the electrons would be small at best. The 
principal change induced is that the relation (5) would be altered to «,=4e,. 

Actually, the formal situation would be improved if we were to assume 
that the units formed during additive coloration are combinations 
containing two halogen ion vacancies and two electrons, that is diatomic 
units analogous to those which Seitz has ascribed to R-centres in the alkali 
halides. In this case the reaction replacing (1) would be a monomolecular 
one ; Mollwo’s observed relation would imply that the heat of solution 
of a free calcium atom in CaF, is very small and that the parameter analo- 
gous to x is about 107. This value still seems large, but is not preposter- 
ously so. In lieu of more definitive experimental work, it seems wisest to 
regard the issue as an unsolved problem. 

It is interesting to discuss the energy of solution of a calcium atom in 
CaF, obtained in the foregoing way in more detail. Let us first regard 
the F-centres in the alkali halides as if formed from free atoms of alkali 
metal by a sequence in which the metal atoms are ionized, the positive 
ions are incorporated in the lattice in normal positive-ion sites and the 
electrons are substituted for halogen ions. For comparative purposes we 
shall assume that the electron behaves like a halogen ion. Under these 
circumstances the energy of formation of an F-centre from the neutral 
free metal, regarding endothermic energies as positive, would be #;—H, 
in which £, is the ionization energy of the atom and //; is the energy of 
dissociation of the lattice into free ions. The appropriate values for 
NaCl, KCl and LiF are given in the table, along with the observed values 
for the first three cases. The value for CaF, is also presented, H, being 
the ionization energy for removal of two electrons from Ca, and H, 
being again the energy of dissociation of the crystal into free ions. 
The ‘observed’ value for CaF, is that calculated from Mollwo’s obser- 
vations, quoted above by the method leading to eqn. (5). We note that 
the estimated values are too low by an amount between about 2-5 and 
3-0 ev in the first three cases, for which observed values are available. 
The discrepancies presumably originate in the fact that the electron 1s 
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broadly distributed and does not use the cavity associated with the 
F-centre as effectively as halogen ion would. The discrepancy seems to 
be about 7 ev per F-centre in the case of CaF, if we assume the bimolecular 
relation (1). Unfortunately a comparable value for LiF is not available. 
The discrepancy would not be altered much if we assume the centres are 
diatomic. On the whole, it appears that the electron associated with the 
F-centre in CaF’, may be distributed even less favourably than in the alkali 
halides. This may originate in the fact that the positive ions are arranged 
tetrahedrally about the F-centre and that, relative to the positive ions, 
the neighbouring negative ions are proportionally much closer to the 
vacancy in CaF, than in the alkali halides. The ratio of distances is 
2/4/3=1-16 in the first case and 1/2=1-41 in the second. Thus the wave 
functions of the F-centre electrons may vary much more rapidly with 
angle in CaF, than in the alkali halides. 


Estimated and observed values of the energy of solution of metallic atoms 
in salts to form F-centres. Endothermic values are regarded as 
positive (expressed in ev). (R6gener 1937, Mollwo 1933, 1934) 


Hy Hy Hy,-Hy, Observed 
NaCl 5-2 7-9 —2-7 0-0 
KCl 4-0) 7-1 —3-1 —0-] 
KBr 4:0 6-9 —2-9 —0-25 
Lif 5:4 10-4 —5-0 — 
CaF, 17-9 30 —1]2-] (+2-10) 


Mollwo quotes work of Wohler and Rodewald (1909) showing that 
essentially half the fluorine ions can be replaced by electrons. Presumably 
the resulting lattice should possess the zinc blende structure and the 
electrons should be distributed somewhat in the nature of a gas through 
the interstitial regions of the lattice. 


§3. CLIMB OF A GENERAL DISLOCATION WITH A BURGERS VECTOR EQUAL 
TO AN ALLOWED TRANSLATION 

We may summarize the discussion of §§2 and 3 by stating that the 
specimens of CaF, which have been additively coloured near 1200°c 
behave as if the dislocations present in them have undergone climb. In 
particular some of the dislocations have developed helical character. It 
is not certain, but highly probable, that the climb is associated with the 
generation of an appreciable density of vacancies in the calcium ion lattice 
as well as in the fluorine lattice. It also seems highly probable that the 
process of additive coloration has enhanced the climb, although some of 
it may be connected with the generation of lattice defects. 

It is evident that many types of reactions in solids which involve the 
lattice defects can produce climb. We shall restrict attention to the case 


Formation and Properties of Helical Dislocations 363 


of a dislocation whose Burgers vector is an allowed translational vector 
so that the climb does not generate a surface of misfit. This case seems 
to be of particular interest in CaF,. We shall also focus attention on 
dislocations which have a strong component of screw character, for these 
appear to be the natural sources of the helical dislocations observed in CaF,. 

Consider a straight dislocation line, running between two node points 
A and B, which possesses a strong screw component, that is which does 
not lie entirely in a plane PN normal to the Burgers vector. For the 
moment we shall not be concerned with the nature of the node points, but 
shall only assume that they are fixed. The line possesses a projection 
A’B in the plane normal to the Burgers vector passing through the node 
B (see fig. 1). Motion of the dislocation in the plane containing AA’B is 
possible without climb, being essentially a form of slip or glide. Actually 
such glide may require an activation energy unless the plane AA’B is a 
special one. 


5, 


---—-------+~------- 


A straight-line dislocation and its projection in the plane PN normal to the 
7 ‘ are 

Burgers vector b. A and B are the node points whereas A’ and B are 

their projections in the plane PN, which has been chosen to pass through 


B for convenience. 


Motion of the dislocation out of the plane AA‘B corresponds to climb 
and requires addition or subtraction of ions from the dislocation, the 
direction of displacement being opposite in the two cases, A direct 
measure of the amount of material added or subtracted is determined by 
the area F in the plane PN lying between A’B and the projection of the 
dislocation line on this plane (fig. 2), the volume of the added material 
being Fb, where 6 is the magnitude of the Burgers weanes b. In Sige ‘ 
which the projected curve between A’B in plane NP winds we a spiral, 
the total area under the multiple-valued curve must be employed in oe 
mining the amount of material involved in the climb. We shall oak 
lindrical sheet generated by the lines parallel 


articular interest in the cy eae 
L In the initial case in which the 


to b passing through the dislocation line. te 
dislocation line is straight, this sheet is simply the plane AA’B, however 
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the cylinder becomes curved once the dislocation climbs if the node points 
remain essentially fixed. Motion of the dislocation along the cylinder, 
parallel to b, can occur by glide*without the addition of material to the 
dislocation. Again, such prismatic glide or slip may require an activation 
energy. 

The energy per unit length of a dislocation depends upon its orientation 
relative to the Burgers vector, even in the isotropic case, and on its orien- 
tation relative to the general crystallographic axis as well in the aniso- 
tropic case. The second case has been discussed in detail for NaCl by 
Huntington ef al. (1955) for NaCl, in which the anisotropic results do not 
appear to be radically different from those for the isotropic case. 


bo 


Fig. 


The change induced in the dislocation shown in fig. 1 by climb. The node points 
have remained fixed. The dislocation line is curved in space; its pro- 
jection in the plane PN is also curved. The area F is proportional to the 
amount of material added or lost in climb. The dislocation can move 
along the cylinder generated by lines parallel to b, and passing through it, 
without the addition of material, provided the nodes remain fixed. This 
corresponds to prismatic slip. 


Consider a segment of a dislocation which possesses a radius of curvature 
of magnitude & directed normal to the Burgers vector, so that the curva- 
ture is maintained by a supersaturation of the agent (eg. Schottky defects) 
which cause the climb. It is readily found that the superconcentration 
required to maintain the local curvature is given by the relation (Bardeen 
and Herring 1951), 
lire kT log (N/N, 


b = 
R i) pe ‘ 4 . . ° ‘ 5 (7) 


in which N is the concentration of the agent inducing climb, N, is the 
equilibrium concentration for the given temperature, ¢ is the energy per 
unit length of the segment, which is assumed to be short compared to 
R, w is the increment of volume per deposited unit inducing climb and b 
is the magnitude of the Burgers vector. We shall focus attention, for 
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simplicity, on a type of ‘zero-order’ approximation in which e can be 
regarded as essentially constant. We then note that there is no restriction 
on the climb of the dislocation if the ratio N /No is sufficiently large that 
R,<AB. It is interesting to observe that the dislocation having a strong 
screw character possesses a great advantage for the process of climb in the 


case in which 
AB>R,>A’'B 


where R, is the critical radius for climb, equal to the reciprocal of the 
right-hand side of (7), if N is the actual concentration. In this case a 
pure edge dislocation of length A’B could not climb beyond a limiting 
point, whereas the mixed dislocation could climb indefinitely. 


Fig. 3 


A climb, more general than that shown in fig. 2, in which the dislocation has 
become helical in form. Its projection on the plane PN is a double 
spiral centred about the points A’ and B. The amount of material added 
is determined by the area under the multiple-valued double spiral in the 
plane PN. The dislocation can still undergo prismatic slip along the 
cylinder parallel to b. 

The ultimate shape of a climbing dislocation which has a strong screw 
component initially seems to be governed by the following factors : 

(1) Magnitude of the supersaturation. 

(2) Rate of arrival to the dislocation of the agent which determined 
climb. 

(3) Rate of redistribution of such material along the dislocation line. 

(4) Ease with which prismatic glide may occur along the cylinder 

determined by lines parallel to b. 

(5) The variation of ¢ with direction relative both to b and to the erystal- 
lographic axes. 

(6) The interaction of the various elements of the dislocation with one 
another through the elastic field. 
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It seems clear that if the material is redistributed along the dislocation 
at arate that is rapid relative to that of arrival, either because the diffusion 
coefficient along the dislocation is very rapid compared to that within the 
bulk material or because the supersaturation is very small, the dislocation 
will tend to assume the shape of lowest energy concomitant with the bound 
points and projected area. If « is a constant we might expect this to be 


Fig. 4 


(Db) 
/ 
i. PoOTections 17 P Varuna ara cs : p e 
The projections in PN corresponding to two interesting forms of helical dis- 
eles In fig. 4 (a) the points A’ and B are inside most of the loops 
0 the double spiral, whereas in fig. 4 (b) only one arm rotates about the 
projection of each node point, the two arms being connected by a single 
c ‘ bi oS 


OP. te. ah ee eae : 
loop. The dislocation shown in fig. 3 corresponds to a case such as that 
illustrated in fig. 4 (0). 
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the are of a circle. Lf N/N 4 is sufficiently small the limiting radius would 
be determined by (7) but the radius could be arbitrarily large if R,<AB. 
A form other than circular will occur in this case if depends strongly on 
orientation. 

On the other hand, if the material arrives at the dislocation at a sub- 
stantially greater rate than it is redistributed along the dislocation, as 
can occur if there is a relatively high superconcentration and if the rate 
of diffusion along the dislocation is not exceedingly high relative to the 
diffusion of the agent inducing climb in the bulk material. the amount of 
climb will tend to be the same over the entire length. The projection of 
the dislocation on the plane PN will ultimately become a double spiral 
with the points A’ and B serving as centres of each half of the double 
spiral (fig. 3). In three dimensions, the dislocation will become a form of 
helix in which the radius and pitch varies along the length in the general 
case. The winding of the helix about an axis parallel to the Burgers 
vector occurs preferentially near the nodes, that is near the centres of 
the double spiral, for one obtains the greatest angle of rotation there 
for a given increment of material, much as in spiral growth of crystals 
from the vapour or solution. 

It is readily seen that if prismatic slip is difficult, the helical dislocation 
will eventually develop lower pitch near the nodes than in intermediate 
regions, as additional turns are added. On the other hand, if prismatic 
glide is relatively easy, the arms of the helix will tend to move toward the 
centre from both ends as they are generated so that the pitch will become 
relatively uniform over the entire length. There are two interesting cases 
to consider: In one the central points A’ of B of the double spiral in the 
plane PN are well within the arms of the double spiral (fig. 4 (a)). In this 
case the outer loops are free to glide parallel to the Burgers vector in such 
a way as to form a relatively simple helix which winds around the original 
dislocation, provided prismatic glide is possible. In the second case 
(fig. 4(b)) most of the arms of the double spiral have a radius less than A’B 
relative to the appropriate centres of rotation. In such a case the loops, 
on undergoing prismatic glide, will form two helices which are inter- 
connected by a loop. It is evident that the two cases become equivalent 
if the radius of curvature of the loops grow larger than A’B. 

Presumably the interaction between the various elements of length of 
the entire dislocation will determine the amount of prismatic slip which 
occurs if glide is easy. In general it appears that neighbouring portions 
of successive turns will repel one another, much like edge dislocations of 
the same sign moving in the same slip plane. On the other hand, dia- 
metrically opposite sides of a given turn or successive turns should behave 
like dislocations of opposite sign and attract one another, at least as far as 
motion in the direction of the Burgers vector is concerned. The problem 
of the long range elastic interaction evidently is worthy of quantitative 
study. ; 

It seems worthwhile noting that as the helical disiageton winds and 
expands it automatically sweeps through a volume of space. This process 
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may aid the speed of growth by placing the dislocation into new regions 
where the supersaturation of the agent causing climb is large as the old 
regions become depleted. 

The helical dislocations which have been observed thus far in this 
laboratory evidently were formed from dislocations which were almost 
of pure screw type running in one of the (110) directions. It appears 
that during an important part of the period of growth of the helix the 
conditions favoured deposition along the length rather than redistribution 
along the dislocation line. Thus the regions near the nodes were able to 
make a number of windings. It is also evident that prismatic slip was 
possible during at least part of the history, for successive turns of the 
helices tend to be evenly spaced, as if slip occurred with comparative ease 
during the process of climb (Plate (a)). 


Fig. 5 


The conversion of a portion of a helical dislocation into single-valued prismatic 
loops. In the case shown, there were present originally two parallel 
dislocations primarily of screw character (horizontal dashed lines), the 
upper one being longer than the lower. It is assumed that the screws 
wind in opposite directions, so that they would have the same sign if 
connected by an edge at the left hand end. Both dislocations have 
developed helical character as a result of climb. The helices in the 
overlapping sections have combined to form single valued loops. 


to {ol 


The step-height of the helices is usually of the same order of magnitude 
as the diameter of the loops. Moreover a major fraction of ae turns 
appear to be uniform in diameter and spacing in many of the helices, 
although an abnormally large or abnormally small turn is not unusual. 
The general uniformity of many of the helices suggest that the prismatic 
glide occurred from the ends along with climb so that the turns formed 
earliest were gradually pushed toward the regions midway between the 
nodes where the superconcentration was fairly homogeneous and gradually 
dropped, perhaps because of the high density of sinks as represented by 
the helix itself. It is quite conceivable that winding up also takes place 
simultaneously at points of the screw dislocation other than the endpoints. 
A small local deviation from the pure serew orientation is all that is 
needed to give rise to winding up. This, with the help of prismatic elide 
would ensure an equal spacing of loops more rapidly. en 

The relatively homogeneous helix frequently terminates rather abruptly 
one or more relatively straight dislocations running off. There does nee 
seem to be evidence for segments near the ends which have diminishing 
radius of curvature and which were the most rapidly turning portions 
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during growth. It is possible that such regions existed during growth 
but ‘unwound * and disappeared after the superconcentration subsided 
because the curvature there was too high. 

Occasionally a helix terminates and is followed, along the same axis, 
which corresponds to the line of the original dislocation, by a sequence, 
of closed loops having about the diameter and spacing of the successive 
turns of the helix (fig. 5) (Plate (b)). It is not certain how such patterns 
are generated, however the result is what one would expect if the original 
dislocation had running parallel to it for a portion of its length a screw 
dislocation of opposite sign (fig. 5) which generated a helix winding in the 
opposite direction. Ultimately the two helices will meet as their diameters 
grow. The regions of contact will join with annihilation of the screw 
component and closed prismatic dislocations or loops will result. There 
seems to be no evident reason why the two screw dislocations of opposite 
sign should run so nearly parallel to one another. 

In connection with the preceding discussion, it may be worth noting 
that one of the simplest ways of promoting the irreversible climb of dis- 
locations would seem to be to establish a temperature gradient in the 
specimen at a temperature where the defects causing climb are readily 
generated at the dislocations. As a result of the thermal gradient, the 
defects should diffuse from one region of the specimen to another causing 
the dislocations in different parts to climb in opposite directions. This 
procedure might prove particularly fruitful in materials such as silicon 
and germanium in which the lattice defects have a high activation energy 
for formation. In such cases one might expect to observe relatively 
open helices such as those found in CaF. The techniques developed by 
Dash (1956) of the General Electric Research Laboratories regarding 
decorating dislocations in silicon with copper might prove useful in dis- 
playing the effect in that material. 


$4, THE GENERATION OF WHISKERS FROM HELICES 


As noted in §3, the helical dislocations found in CaF, probably were 
formed under circumstances in which the factors inducing climb were not 
as great as they might be because of the intervention of electrolytical 
currents associated with Frenkel disorder in the halogen lattice. 
Apparently only a small fraction of the defects which play a role in the 
additive coloration originate in a process requiring climb. This probably 
explains why the helices found have a step height that is large compared 
to the length of the Burgers vector. 

It is clear that the situation would be very different if a vigorous 
reaction that occurred at the surface of a crystal consumed lattice defects 
which could be generated only by climb in the interior of the specimen. 
In this event all sources would operate more intensely. In particular the 
helical dislocations of the type observed in this laboratory could become 
much more tightly wound. Under the proper circumstances such spiral 
prismatic dislocations could generate a type of whisker on the surface. 
We shall discuss some of the properties of such whiskers which may be 
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anticipated. It will become clear that they are of a type which would 
grow from a base within the solid and seem on the whole to be like the 
whiskers which have been grown on tin and other coatings of metals under 
specialized conditions, such as when in an appropriate reactive atmo- 
sphere or when squeezed in a vice. 

It should be emphasized that two general types of whiskers have 
been prepared and studied. The first type, which grows from the base 
was first made the object of careful physical study by Herring and Galt 
(1952) after observations of technological interest by Cobb (1946) and by 
Compton et al. (1951). They have been produced by mechanical strain 
by Fisher ef al. (1954). The mode of growth has been the object of 
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Schematic picture showing the development of a whisker from a spiral or helical 
prismatic dislocation. In fig. 6 (a) the upper horizontal full line repre- 
sents the surface of the specimen. The vertical line 4B corresponds to 
the screw dislocations whose Burgers vector b is shown. The horizontal 
dashed line above the full line representing the surface and running to the 
left from A represents the ridge on the surface associated with the screw 
dislocation. In fig. 6 (b), the screw has climbed into a helical form as a 
result of the addition of material. In 6 (c), some of the loops of the helix 
have been attracted to the surface, moving by prismatic slip, and have 
formed a cylindrical protrusion or whisker. The whisker can continue 
to grow by addition of loops from beneath the surface. In the case 
shown it is assumed that the prismatic slip of the helix to the surface has 
taken place in such a way that the ridge on the surface (horizontal 
dashed line) is cut and part is raised to the top of the whisker. In this 
event, the whisker would possess a screw along its axis (see fig. 7). 
Actually one would expect the screw to be unstable at the wpper end of 
the whisker in this idealized case and to slip out of the cylinder, at least 
at sufficiently high temperature, 


Formation and Properties of Helical Dislocations 371 


systematic study by Baker (unpublished); and by Levy and Kammerer 
(1955). The second type of whisker is grown. by deposition of material 
from the vapour phase or from solution and has been the object of study 
by many investigators in recent years (Sears 1953, Sears et al. 1954, 
Johnson and Amick 1954, Piper and Roth 1953, Gyulai 1954, Brenner 
and Morelock 1956, Brenner and Sears 1956). The mechanism of for- 
mation of the second type of whisker appear to be closely related to that 
of crystal growth from the vapour or from solution. Such whiskers 
will not be of primary interest to us here. Mechanisms for the growth of 
the first type of whisker have been proposed by Peach (1952), Eshelby 
(1953), Frank (1953) and Seitz (comment by Eshelby 1953). The mecha- 
nism of Peach does not provide an explanation of the fact they grow from 
the base, whereas the other three mechanisms do. Eshelby postulated 
an internal source of simple prismatic dislocations, such as that proposed 
by Bardeen and Herring (1951), and suggested that the prismatic dis- 
locations slip to the surface, being guided by a stress field. Frank’s 


Fig. 7 


(a) (b) 


A tightly wound helical or spiral prismatic dislocation for which the step height 
is equal to the Burgers vector (fig. 7 (a)) is equivalent to a stack of discs, 
or simple prismatic loops, through which a screw dislocation runs 
(fig. 7 (b)) when the helix protrudes from the surface. The screw is 
represented by the vertical dashed line running from B to A. 


+ Unpublished work undertaken in cooperation with J. S. Koehler. Baker 
carried out an extensive study of whiskers formed on coatings of tin and 
Pb-Sn alloy upon iron, on solid Pb-Sn alloy, on layers of cadmium eee 
upon copper and iron, on layers of Cd-Hg electroplated upon copper Coie 
with mercury, on electroplated zinc, and on a thin layer of AgNO, upon gl ass. 
Whiskers grew best on electroplated layers on the whole and not as a ou 
the thicker hot-dipped layers. A few of Baker’s numerous valuable conclusions 
are summarized in the text of the present paper. 
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mechanism proposes that the whiskers are formed by the repeated circular 
motion of the terminal end of a screw dislocation that cuts the surface. 
This mechanism is somewhat related to the processes described in this 
paper. Seitz proposed formation by the action of a helical, or spiral, 
prismatic dislocation along the lines discussed in more detail here. The 
authors are much indebted to Drs. Baker and Levy for the opportunity 
of reading unpublished work. 

We shall assume that the thermodynamic driving force inducing climb 
is very great so that climb is not halted merely by forces between dis- 
locations which act over the normal distances involved in networks, of 
the order of 10-4cm. We shall consider the case in which the climb adds 
material to the dislocation, that is, interstitial atoms condense or vacancies 
evaporate. 


Fig. 8 


An alternative direct method for forming a whisker without a screw on its axis 
from a helical dislocation. The surface of the specimens is shown along 
with the internal spiral prismatic dislocation and the ridge associated 
with the termination of a dislocation having the associated screw 
character on the surface. The ridge is straight on the left-hand side of 
the diagram but curves along the line corresponding to the projection of 
the helix on the surface and dwindles to zero at the point where the helix 
reaches the surface. As the helix moves to the surface by prismatic 
slip, the ridge extends around the loop on the surface corresponding to 
the projection of the helix. Once each turn a complete disc of crystal, 
one Burgers unit high, will be formed on the surface. These discs will 
form a whisker which does not have a screw dislocation. This process 
probably does not occur when whiskers are formed rapidly. 


The simplest case to consider in a geometrical way is that shown in 
fig. 6(a). Here a screw dislocation terminates at a ridge on the surface 
and at a node point inside the solid. The two terminal points may be 
regarded as equivalent to A and B in fig. 1. More exactly, the terminal 
end conjugate to B could be regarded as a point outside the lattice which 
is the image of B in the surface. If this dislocation is induced to climb 
into a helix by the development of loops at the bottom point, inside the 
solid, and by prismatic glide, the form could, at least in an idealized case 
achieve the state indicated schematically in fig. 6 (6). The neighbouring 
loops of the helix repel one another, but should be attracted to the surface 
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as a result of the image forces. Asa result, the system could achieve the 
state shown in fig. 6 (c). In this schematic case it has been assumed that 
in breaking through the surface the loops have cut the ridge and raised 
the segment running to the centre of the cylinder. Elementary reasoning 
shows that the protruding helix in which neighbouring turns are now 
separated by the Burgers vector, is equivalent to a stack of discs, or a rod, 
having a screw dislocation running down its centre (fig. 7). Thus in the 
idealized case shown, the whisker would possess a screw dislocation along 
its axis. Actually the screw could probably be very unstable in the case 
shown, for the ‘bend-point’ A in figs. 6 (c) and 7(b) should be attracted 
to the side of the cylinder and should pull the upper part of the screw 
dislocation with it. This may be expected to occur in the case illustrated 
when the cylinder has protruded by a distance of the order of or larger 
than its diameter. Thus the whisker would probably lose its screw dis- 
location in the course of growth. 


A part of a dislocation network which terminates at a surface. In the case 
illustrated, all the dislocations in the network are assumed to have screw 
character and hence may be converted into helical dislocations as a result 
of climb. Dislocation 1 is in an ideal position to form a whisker. The 
spiral prismatic dislocations formed from 2 and 4 are poorly oriented to 
break through the surface, however those generated from 3 and 5 could 
eventually force their way outward by compelling the node points I and 
II to move. In this case the initial whisker could develop two kinks as 
the two nodes emerge, assuming the growth remains continuous (see 
figs. 10 and 11). 


In a somewhat more idealized case, shown in fig. 8, the whisker would 
grow from start without a screw dislocation along its axis. In this case 
the terminal point of the helix at the surface rotates around the circle 
marking the intersection of the whisker with the original surface of the 
solid, making one revolution for each loop which protrudes. In effect, a 
complete disc of solid one Burgers unit high is formed at the base of the 
protruding cylinder for each revolution of the point where the helical 
dislocation emerges at the surface. 

The experimental investigations of whiskers grown from the base indi- 
cate that few conform to the ideal situation shown in these figures. Most, 
_ appear to have an irregular initial, or head ’ portion which is not tightly 
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attached. Frequently a large fraction are kinked, the kink angle being 
more or less well defined. Levy has observed cases, relatively rare, in 
which successive kinks have the same sign so that the whiskers have a 
spiral form. All investigators appear to agree that the kinks formed 
during the growth process and not after the whisker has emerged. 


Fig. 10 


The case in which dislocation 3 shown in fig. 9 has formed its own whisker 
(schematic). The bars projecting from the surface represent the whiskers, 
directed more or less along the Burgers vectors in the two case. 


Fig. 11 


Tn this case node I shown in fig. 9 has moved to the surface as a result of climb 
and the pressure exerted by the underlying prismatic dislocations. The 
whisker produced by dislocation 1 continues growth, but with a kink each 
time a node emerges. 


Actually, the mechanism of whisker growth based on the assumption 
that they originate in helical dislocations formed by the climb of those 
having Burgers character seems to have ample room for describing the 
observed phenomena if we take into account the fact that most of the 
dislocations near the surface are involved in network structures and do 
not have the simple form shown in fig. 6 (a). Consider the more realistic 
situation shown in fig. 9 in which a part of a dislocation network which 
terminates at the surface is shown. For illustrative purposes five dis- 
locations having mainly screw character are shown, number 1 terminating 
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at the surface. Dislocations 2, 4 and 5 link with other parts of the net- 
work not shown. All should be transformed to helices if the environment 
favours climb. The helix formed by dislocation 1 could produce a rela- 
tively simple whisker of the type described above, but growing diagonally 
from the surface in the direction of the Burgers vector. The direction of 
glide for dislocations 2 and 4 is poorly situated for formation of whiskers 
on the surface shown. On the other hand, the Burgers vectors for dis- 
locations 3 and 5 are disposed in such a way that the helices formed from 
them could produce whiskers if the helices could reach the surface. 
Because of the repulsive forces between successive turns of the helix, the 
helices will exert forces upon the node points I and II in the directions one 
might expect from compressed springs. If the forces between turns of 
the helix formed from dislocation 3 become strong enough they could 
glide to the surface in the direction of the Burgers vector and form a new 
whisker (fig. 10). On the other hand, it is also possible that the node I 
would be induced to climb to the surface as a result of the forces upon it 
(fig. 11), and the opportunity for climb inherent in the situation, so that 
the whisker formed from dislocation 1 would eventually stop growing and 
be joined to a whisker, oriented in a different direction formed from dis- 
location 3. The angle of kink should be closely related to the angle 
between the Burgers vectors of dislocations 1 and 3, although it need not 
be exactly the same if the whisker has not grown parallel to the Burgers 
vector on the surface. 


Fig. 12 


ees 


A dislocation pattern which may produce a Y whisker, with its base on the 
surface, of the type observed by Baker. The two whiskers produced by 
the upper dislocations join into one as the node point is forced to the 
surface as a result of climb and the action of the lower one. 


Similarly, the helix formed from dislocation 5 could eventually force 
node II to the surface and produce a second kink. Thus the ultimate 
kinked structure of the whisker would reflect, in a somewhat crude way, 
a part of the internal structure of the network. era. 

It is not clear to what extent the dislocations formed in this way would 
have a screw dislocation running down the centre of each segment 
between kinks. In the case illustrated, we might expect the part of the 
whisker originating in dislocation 1 to be relatively perfect since it 
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conforms somewhat to the situation shown in fig. 6(a). On the other 
hand the screw dislocations may or may not remain entrapped in the 
segments originating in dislocations 3 and 5, depending upon. somewhat 
accidental circumstances. 

If dislocation 1 in fig. 9 does not terminate on the surface, but is tied 
back into the network the helix formed from dislocation 3 would still 
attempt to break through to the surface. This could occur by a type of 
prismatic punching process in which a more or less cylindrical block of 
material between the node I and the surface would be forced out as a 
result of prismatic slip of the helix derived from dislocation 3 in the 
direction of its Burgers vector. Once the helix has reached the surface, 
the node point I would effectively be there also and further whisker 
growth would occur in what may be regarded as a normal way. It is 
possible that the head portions found on many whiskers correspond to the 
segments which are punched through prior to the establishment of con- 
ditions for development of a normal whisker. 

More complex situations than that described evidently are possible. 
Baker has found that pairs of whiskers can move together during growth 
and join to form a Y whose base is attached to the surface. A case of 
this kind could originate from the simple dislocation structure shown in 
fig. 12. In the case quoted by Baker the two upper arms of the Y were 
curved downward toward the surface. Bending of this type, and in fact 
bending of whiskers in general, is not difficult to understand if the nodes 
are induced to move so that the points at which the helices break through 
the surface shift during the course of formation of whiskers. 

Levy’s spirally kinked whiskers seem to require a more complex 
sequence of events. One must assume that a sequence of node points 
distributed in a more or less spiral manner inside the solid are forced to the 
surface in turn and that successive kink angles are in the same direction. 
Such cases are relatively uncommon and may be determined in the main 
by a somewhat chance sequence of events. 

It is clear that the basis for the explanation of the occurrence of kinks, 
i.e. the climbing of a node, does not depend strictly on the validity of 
the particular growth mechanism described here, but can be applied to 
any dislocation mechanism where the axis of the whisker is determined 
by the Burgers vector of the underlying dislocation 

Eshelby (1953) has proposed that the driving force for whiskers of the 
type found on tin plate is the result of internal tensions related to the 
desire of a convex portion of the metal surface to increase in area in order 
to form more oxide. This explanation of the source of climb seems to be 
in good accordance with the fact that whisker growth can also be stimu- 
lated by the action of mechanical stresses, as at the edge of a specimen 
mounted in a vice, for in the second. case it is evident that the dislocations 
in the region just outside the vice desires to accept material from the 
regions under pressure and should climb in the way needed to form 
whiskers. In general, the principal condition required for the formation 
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of whiskers, apart from the presence of appropriate dislocations, is an 
unstable energetic situation which can be brought nearer equilibrium by 
the migration of imperfections and the appropriate climb of dislocations. 

It is found that the whiskers require an incubation time before appearing. 
This presumably is the time needed for the helix to grow to a size at 
which the internal and image forces are sufficiently great to induce it to 
move to the surface. In general, we might expect the diameter at which 
this is achieved to be comparable to the depth of the source below the 
surface, that is to the linear spacing in the network. Baker has observed 
that the whiskers, although normally quite strong under bending or 
tensile stress, would flatten relatively easily if compressed normal to the 
axis of the whisker. This observation would bear further investigation 
since it seems to imply the presence of sources of dislocations possessing a 
Burgers vector normal to the whisker axis. The discussion given above 
would not appear to offer a simple explanation of the origin of such sources. 
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ABSTRACT 


Increases in the electrical resistivity of gold wires due to quenching 
have been measured and interpreted in terms of vacancies. The energy 
for vacancy formation was deduced to be 0-95 ev and the concentration 
of vacancies at the melting point 610-4. Annealing measurements 
between 70° and 130°c indicated that the movement energy of a vacancy 
was 0-68 ev and that it had on average a life of 108-10 jumps. 


§ 1. INTRODUCTION 


WE have recently given the results of experiments on the changes in 
electrical resistivity of platinum wires after quenching from high temper- 
atures and after annealing (Bradshaw and Pearson 1956). The increases 
in resistivity caused by quenching were interpreted as being due to 
vacancies. Since the published results of similar experiments on gold 
were somewhat conflicting (see table) we have also made measurements on 
gold wires using the same methods as before. 

The wires were from Johnson Matthey ‘Specpure’ material and were 
0-2 and 0-14 mm dia. They were heated by electric current to tempera- 
tures between 560° and 960°c and quenched into water at ~ 4°c. Their 
annealing characteristics were studied by immersing them in liquid 
baths at temperatures between 70° and 130°c; to ensure complete 
removal of vacancies at the end of an anneal they were heated briefly to 
~ 400°c. The changes of resistance following quenching and at appro- 
priate stages during annealing were found by measurements in liquid 
nitrogen. Potential leads 0-02 mm dia. were welded on for this purpose. 
The temperatures of the wires just before quenching were obtained by 
measuring the hot resistance and using resistance/temperature data 
provided by the N.P.L. (Powell, private communication). The gold wires 
were welded on to platinum leads whose diameters were chosen so as to 
make the temperature gradient between the potential leads on the gold 
negligible at the quenching temperature. 
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§ 2. QUENCHING 

The values of Apo, the increase in resistivity resulting from quenching 
from different temperatures, were plotted on the assumption that the 
usual relation would hold, 

v=A exp (—4A,/kT) 

where v is the vacancy concentration, 4H, is the energy required to form 
a vacancy and A is the entropy factor. Figure 1 shows a plot of 4 Po 
against 1/7’ for the two wire sizes. Each point is the mean of about six 


Fig. 1 
T (°C) 


600 650 700 750 800 850 900 950 
0.05 


0-0l 


AP, (out cm) 
0) 
a 


© 0-14 MM DIA 
x O0:20MM DIA 


0.00] 


12 " 10 9 8 
4 
ie (aon) 
Increase of resistivity by quenching. 


observations. Using Abelés’ (1955) estimate for the resistivity of 
vacancies (1 atm °% corresponds to 1-45 1 ohm cm) the values of A from 
the lower temperature linear portions of the graph are found, for the 0-2 
and 0-14 mm wires, to be 1-9 and 3-1. The values of 4H, were 0-92 and 
0:98 ev. 

Assuming that the same annealing process for vacancies operated at 
the higher temperatures during quenching as during the annealing 
experiments (see below) we have calculated, as in the paper on. platinum, 
the transition temperature above which one might expect to lose 
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vacancies on quenching owing to their high mobility. The initial 
quenching rates were estimated from the previous work on platinum to 
be 30°/msec (0-2 mm) and 60°/msec (0-14 mm). The transition temper- 
atures were calculated to be 1/(10-2 10-4) and Li(9-9peL05*) Ke tis 
interesting that both lines appear to be straight to temperatures higher 
than the calculated transition temperatures, although one would expect 
them to start curving some way below this. (A similar result may be 
deduced from the figures given by Bauerle et al. 1956.) The calculated 
transition temperatures would have been compatible with what was 
observed if 4H, had been ~ 0-57 ev ; this is outside our experimental 
scatter. The best evidence for adequately rapid quenching and hence 
correct values of 4H, is the straightness of the lines in the figure at lower 
temperatures and the agreement between the results for different wire 
sizes. The shift between the lines is thought to be due to slight errors in 
estimating the temperature. The mean figures for formation energy 
and entropy factor are given in the table. We suggest that the total 
error in 4H, is probably + 0-lev. Also given with the results previously 
published for gold are values of A calculated from the data provided. 
We agree reasonably with the figures of Bauerle e¢ al. and would predict a 
concentration of 6x 10~* vacancies at the melting point. 


AH, ev A AH, ev nB 
Kauffmann and Koehler 1-2 — 0-4 +0-14 — 
(1952) 
Kauffmann and Koehler 1-28 +0-03 12 0-68 +0-03 1:0 x 10-8 
(1955) 
Bauerle et al. (1956) 1:02+0:06 | 5-1 0-66+0-06 | 0-3x10-8§ 
Lazarev and 0:79 0:2 0-52 4-0 x 10-U 
Ovcharenko (1955) 
Present work 0:95 +0-1 2-4 0-68 +0-05 1:0 x 10-8 


§ 3. ANNEALING 


Most anneals were carried out on wires quenched from temperatures 
in the region of 850°c and, as with platinum, the decay of 4p, the increase 
in resistivity, with time happened in a single stage. There were measur- 
able differences in the annealing rates even for the same wire with the 
same Ap, and annealing temperature and, for different wires, character- 
istic differences in the shape of the decay curves were also. observed 
(see fig. 2). However, the general curve shape usually approximated to 
that of a second order reaction (down to ~ 15%, App remaining). To see 
whether varying 4p, had any effect, a few anneals were Spee at the Bane 
temperature on the same wires quenched from ~710°o and ~ 960°c. 
For a wire quenched from ~ 960°c the annealing curve shape was about 
the same as that for wires quenched from ~ 850°c but for wires quenched 
from ~ 710°c the initial slope was only about half that for wires quenched 
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from the higher temperatures and the curve shape approached that of a 


third order reaction (fig. 3). 
Fig. 2 
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If n is the probability per jump that a vacancy is lost, and if despite 
the shape of the annealing curve we interpret annealing in terms of 
dl 
i “ =12 nv,B exp (—AH,/kT) 
where vp is the atomic frequency (assumed here to be 101% per sec), B is 
the entropy factor for vacancy movement (probably in the region between 
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1 and 10) and JAH, is the activation energy for movement, then the 
differences in annealing can be regarded as being due to differences in n 
at any instantf. 

4H, was obtained by altering the annealing temperature to a new 
value when 4p had decayed to ~ 44p, and noting the change of annealing 
rate. Values of 4H, obtained by straightforward complete anneals at 
different temperatures in the cases where n appeared to be stabilized did 
not differ significantly. The mean value of 4H, for both wire sizes was 
0-68 ev; it had an S.D. of 0:02 ev. The total error limits were perhaps 
+0:05ev. The mean value of nB was 1:010-® (for quenches from 
~ 850°c and when 4p = $4p,), though this figure should only be taken as 
giving an order of magnitude. These figures are compared in the table 
with those of other workers (nB calculated from their figures). The 
grain sizes varied from those equal to the full wire diameter down to 
~ 0-03 mm diameter. Sizes of this order imply that the grain boundaries 
and surfaces were not important as sinks for vacancies since assuming 
B is between 1 and 10 they would lead to values of n considerably smaller 
than the mean value of n obtained experimentally. 

Despite the shapes of the annealing curves approximating to second 
or third order a value of n lying between 10-* and 10-® when v was 10-4 
would appear to preclude the idea of annealing by vacancies pairing or 
clustering. As with platinum the picture of dislocations providing the 
sinks seems to be more plausible, though the observed annealing behaviour 
indicates that the process is not a simple one. 

The predicted value of 1-63 ev for self-diffusion activation energy 
agrees reasonably with the most recent experimental value of 1-70 ev 
(Okkerse 1956). 
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+ The factor 12 arises since there are 12 atoms round the vacancy which 
can jump into it. We are grateful to Dr. B. 8. Berry for pointing out that 
the factor was omitted in a previous paper on Pt. This however does not 
affect the conclusions therein except that n for Pt should have been ~1023; 


not ~10-°. 
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ABSTRACT 


A study of the factors controlling the coefficient of volume expansion of 
a pure metal, in a previous paper, has led to the conclusion that there are 
two isomorphous states, differing only in volume, of any given crystal 
structure. In this paper the relative stabilities of these phases is dis- 
cussed and the possibility of isomorphic transformations is examined. A 
satisfactory description of a paradoxical phase change observed by 
Bridgman in the case of caesium under high pressures is achieved. An 
isomorphic transition observed in cerium appears to be another example 
of this type of phase change which, under suitable conditions imposed by 
the external variables, pressure and temperature, and possibly composition 
in the case of alloys, is possible in all metals. 


$1. INTRODUCTION 


RECENTLY an analysis of the factors which control the thermal expansion 
of a pure substance has been carried out (Varley 1956) in an earlier paper 
hereinafter referred to as loc. cit. It has been found that the equation 
relating the volume coefficients of expansion «, to the electronic and 
vibrational characteristics of a metal (analogous to Griineisen’s relation) 
is quadratic in « ; a«?+ba+c=0. Of the two possible roots for « one is 
positive and the other negative and the root with the larger modulus 
corresponds to the state of lower free energy. The argument for this is 
based on the fact that (C,—C,) is proportional to «? where C,,, C, are the 
total specific heats at constant pressure and volume respectively. There 
is no restriction upon the sign of « so that the stable state can have a 
positive or negative coefficient of volume expansion, dependent upon 
the details of the electronic and lattice vibrations of the solid. 

If, then, one considers a plot of all possible values of the Gibbs free 
energy, G and not only those values obeying the equation of state, as a 
function of volume at constant temperature and pressure, the curve will 
show two minima separated by a maximum. At 0°K the minima and 
maximum will coalesce into a single minimum. The maximum will 
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always occur at a volume V,(P) which is the volume of the substance at 
0°K for a given pressure P. The free energy at the maximum is then the 
free energy of the substance at constant volume V,(P) and temperature 7’, 
That this point must be a maximum can be seen from the argument 
that heating a crystal at constant pressure from absolute zero will result 
in a change of volume and an associated increase in entropy ; and since 
C',—C, is proportional to «? and is greater than zero, G must have a 
lower value than if the substance were held at constant volume, V(P), 
whatever the sign of the change in volume. 

There will be a family of curves of G as a function of V at constant 7 
and P, one for each value of the pressure. The loci of the minima give 
the equilibrium relation between G and P at constant 7. Again, a plot 
of G as a function of 7 at constant P will give a trace of the values of G 
for the two minima which occur in the equilibrium relation between @ 
and V. In other words, the loci of the minima define the equation of 
state of the substance. 

As the pressure upon the substance is increased at constant temperature 
one will expect, by Le Chatelier’s principle, that the state with the smaller 
volume will tend to become more stable. Thus, for a substance with a 
positive coefficient of volume expansion, the state (minimum) with the 
lower volume will have a higher free energy than that with the larger 
volume, at low pressures. As the pressure on the substance is increased 
the family of curves of G as a function of V will be sheared in such a way 
that ultimately the state (minimum) with the smaller volume will have 
the lower free energy. 

For a substance with a negative coefficient of volume expansion similar 
arguments apply with the sign of the pressure reversed, e.g. for a crystal 
in a polycrystalline solid under internal stresses giving rise to negative 
pressures. 

The two states (minima) will have equal free energies when the co- 
efficients of volume expansion are of equal and opposite sign, that is 
when the coefficient of the linear term in the quadratic expression for « 
is zero. If this condition can be achieved by a suitable variation of 
pressure and/or temperature an isomorphic phase change will occur. 


§2. PHasz CHANGES UNDER PRESSURE IN CAESIUM 


It can be shown, loc. cit., that for a solid in which the valence electrons 
may be considered as distributed in a free electron spectrum (Fermi gas), 
the coefficient b, of the linear term is 


Fee AMPED ct 2039 side sy abe eae 


where C’, and C,, are the lattice and electronic specific heats at constant 
volume, respectively, y,=—(0 In 6/0 In V)7, B= — (0 In 6/d In T),, and @ is 
the characteristic temperature defining the thermal entropy of the 


lattice (approximating to the Debye-@). 
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Consider now the metal Caesium. For an alkali metal such as this 
C.<C, at normal temperatures so that the condition for equal free 
energies reduces to y,-+0. Now 6, the characteristic temperature is 
proportional to V1/6 K1/2 for a Debye solid with, in addition, Poisson’s 
ratio independent of volume (Slater 1939). To take into account the 
variation of Poisson’s ratio with volume it will be assumed that this 
quantity varies through K and that 


Occ VEER Rees a te) el eee een a) 
where 7 is a constant. 
Then sy, = (9 In 6/8 In V)p—=—}—n(O In K/@In V)p. . ~ (3) 


For caesium y,—1-29 at normal pressure and temperature. Also, from 
Bridgman’s measurements (1949, p. 180) it is found that K for Cs is 
proportional to V, approximately 
K=K,+m(V,—V)/V, ee ee ee 

where K, and V, are the bulk modulus and volume at normal pressure 
and m is the rate of change of K with V/V,. Using these results for Cs 
it is found that m=+10K,, K,=16-3 kg cem~?. Substitution of this data 
into (3) and using the relation dK/dV=—m/V, gives the value 0-146 
for n. 

The volume at which the two phases have the same free energy can now 
be obtained by substituting the condition y,—0 into (3). This gives 

ViVi=0-59- 

Experimentally Bridgman (1949, p. 181) finds two transformations under 
pressure in caesium, one occurring at a pressure of 23 000 kg cm~? at a 
relative volume V/V )=—0-6284 and a second one at 45000 kg cm? at 
a relative volume V/V,=0-498. One of these transformations (at the 
lower pressure) is attributed to a polymorphic transformation from the 
body-centred cubic structure to the face-centred cubic form. The 
other is a paradox. The above arguments remove this paradox and lead 
one to expect an isomorphous phase change from a structure with a 
relatively large atomic volume to one with a smaller atomic volume. 
Two phases with the same crystal structure but with different atomic 
volumes are always possible, and with a suitable adjustment of external 
variables (P and 7’) a transformation from one to the other of these can 
be brought about in suitable cases. The above calculations suggest that 
the transformation at the lower pressure in caesium is associated with 
this isomorphic transition while the high pressure transformation is due 
to a polymorphic transition to a more closely packed structure. 

Referring back to the quadratic expression for « it will be seen that at 
the transition point the coefficients of expansion «, have the values 
+(c/a)'/?. There will be a discontinuous volume change, therefore, of 


T 
27, |«,] a7’ where V, is the volume of the substance at the transition 
0 


pressure and 0°xK. 
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Assuming, once again, that the valence electrons behave as a Fermi 
gas, it can be shown, loc. cit. that the coefficient C in the equation for « 


is given by An(ep) 1 
=c,( dep dita | 


where n(¢;) is the density of states per unit volume at the Fermi surface. 
Alsoa=KVJ, loc. cit.,and hence A V ,, the volume change at the transition 


point, is 
fl 27,2/N2 
AV —=2V,/ epee a |e: 
0 be, K 


Evaluation for 7=300°K with V=0-59V,, V,;=71 cm? and again 
using the experimental data of Bridgman (1949, p. 180) for K, gives 
AV ,/V,=1-9;x 10-8 which is to be compared with the corresponding 
experimental quantity 610-3. The theoretical calculation is thus of the 
right order of magnitude and supports the suggestion that it is the low 
pressure transformation which is the isomorphic transition since the experi- 
mental relative volume change at the higher pressure transformation is 
5:6 10-*. However, more detailed calculations may indicate the converse 
to be true. 


§3. IsomorPHIc PHASE CHANGES IN CERIUM 


An isomorphic transformation will always be expected if the external 
variables, P and 7’, can be varied in such a manner as to make the linear 
coefficient 6, in the quadratic expression for « pass through zero. The 
transition in caesium has already been discussed. There is another 
example in ceriumy. 

It has been observed (cf. Lawson and Tang 1949) that under high 
pressure and at low temperatures cerium undergoes a transition from one 
face-centred cubic structure to another with a smaller lattice parameter. 
At the transition volume the coefficient 6 in the quadratic expression for 
a% must be zero. 

Now 0 is given by, loc. cit., 


Cry, 2a ne) te (as) (1 a } 
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where N is the number of electrons per unit volume in partially filled 
bands and n(e,,) is the density of states at the Fermi energy, «p. 

Experimentally, Bridgman (1949, p. 168) finds that for cerium the 
change of compressibility y, with volume is negative, and since this 
term is a major factor in y, then y,<.9. The coefficient of expansion of 
cerium at normal pressure is positive so that 6>0 at large volume and 
low pressure. 

If, as a result of the compressibility decreasing with increasing volume, 
y,, becomes more negative on compression then, since the coefficient of 
C,>0 at normal pressures, } will for some value of V’ pass through zero. 
An isomorphic transition will then occur to a phase with a smaller atomic 
volume. 


+ The writer is indebted to Dr. L. Gutman for drawing his attention to this. 
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Again, if the pressure is held constant and the temperature reduced 
the coefficient b may again pass from positive to negative values. This 
is possible if the coefficient of C, is large enough. At very low temper- 
atures, where the electronic term in 6 predominates the structure should 
transform back again to a continuation of its room temperature form. 


§4. ConcLUDING REMARKS 


Isomorphic transitions are first order phase changes and will be accomp- 
anied by a latent heat of transformation. This is evident since there is a 
discontinuous change in volume and hence in internal energy. 

Thus whenever the coefficient of volume expansion changes sign, as it 
must do when the coefficient 6, of the linear term in the quadratic expres- 
sion for « changes sign, an isomorphic transformation will occur with a 
discontinuous change in crystal volume. 

In practice there may be a gradual transition from one structure to the 
other, a complete transformation at any given temperature being hindered 
by the generation of internal strains. 
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SUMMARY 
The magnetic susceptibility has been measured at 77° and 300°x 
of a series of copper—zine alloys, ranging in composition from 0—47-6 atm %, 
zine, and of one silver-cadmium alloy containing 30-5 atm °% cadmium. 
For the copper—zine system the absolute susceptibility ranges from 
—0-085 x 10° e.m.u./g for pure copper to —0-191 x 10-6 e.m.u./g for 47-6 %/, 
zinc. No significant difference was found between the values for a given 
alloy at the two temperatures of measurement. 
The significance of the results is discussed in relation to a theory 
recently proposed by Henry and Rogers. 


§ 1. INTRODUCTION 


TuIs paper deals with the measurement and interpretation of the magnetic 
susceptibility of brass as a function of zinc content over the face-centred 
cubic, «, range and at a single composition in the more restricted, body- 
centred cubic, f, range. <A single measurement on the analogous 
silver-cadmium solid solution is also reported. Two temperatures of 
measurement were used, 77° and 300°K. 

The copper-zine system was chosen because of the relatively large 
amount of experimental information on various properties already 
available. Such a concentration of effort on one system should make 
easier the task of interpreting the behaviour of any one property and of 
formulating a general theory to account for the effects of alloying in 
a solid solution of this type. 

When the work was started the variation of susceptibility, x, through 
the « and £ ranges was known approximately through measurements by 
Weber (1920) and by Endo (1925). These showed that x in these ranges 
decreases with zine concentration though the actual values obtained are 
probably not reliable. Endo also reported values for the y range which 
are larger than those found (Smith 1935, Marcus 1949) by several later 
observers (e.g. —1-02 instead of —0-78 x 10-6 e.m.u./g at 67 atm % zinc). 
Neither worker made any check on the possible errors arising from 
ferromagnetic impurities. 

Since the present investigation was completed the results have been 
published of an extensive series of measurements by Henry and Rogers 


+ Communicated by the Authors. 


2D2 


390 B. G. Childs and J. Penfold on the 


(1956) of the room temperature susceptibility of four copper based 
solutions including copper-zinc. They put forward an interpretation 
of the decrease in x with concentration based on a suggestion by Friedel 
(1952). This is that in calculating the contribution to the diamagnetism 
the solute atoms can be considered as being singly ionized and free from 
interaction with other ions. Henry and Rogers require further that the 
slope of the electronic density of states curve for copper should be negative, 
the decrease in density of states at the Fermi level and hence in the 
Pauli paramagnetism as the electron to atom ratio is increased being 
the chief reason for the decrease in y. The discussion in the present 
paper centres on showing that this latter assumption is probably not 
valid and that some other explanation must be found for the decrease 
in the magnetic susceptibility. 


§ 2. MeTHoD OF MEASUREMENT 


The method of measuring susceptibility was that due to Faraday. 
A specimen placed in a region of non-uniform magnetic field experiences 
a force in, say, the X-direction proportional to the susceptibility and 
to the average, over the specimen, of the quantity H(dH/dX), where H 
is the magnetic field at X. In the present case it was arranged that the 
force acted principally in the vertical direction, the horizontal Y and Z 
components being very small. 

The apparatus (see fig. 1) thus consisted essentially of two sections, 
an electromagnet for producing fields of the required intensity and 
intensity distribution, and a sensitive balance for measuring the force 
exerted on the specimen. This force was of the order of one dyne in the 
present tests. 


2.1. Magnet and Field Calibration 


The magnet was designed to give a maximum field of 15 000 oersted 
under continuous operation, the power consumption being 50 amp at 
100 volt. Cooling was provided by a forced air draught in the case of 
the coils and by water for the magnet core. The pole pieces were 
machined to a cylindrical profile r°/? cos 36/2=constant (r being the 
radius vector and @ the angle to the vertical), calculated by Davy (1942) 
to give an extended region of constant force. The size of this region 
was limited in practice by errors in machining and by the necessarily 
finite cross sectional area of the magnet core. The minimum gap of the 
polepieces was set at 1 in. and their width and maximum vertical height 
at 4 in. 

The magnet was mounted on a trolley rumning on rails which allowed 
it to be removed from the rest of the apparatus when assembling the 
specimen. 

The power supply was from a mercury arc rectifier, the current being 
adjusted manually with a bank of resistances in series with the magnet. 
The order of precisicn required in measurement and control of the current 
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to give an error of less than 0-1° in susceptibility was, for example, 
‘0-01 amp at 20 amps. 

The calibration of field strength H, and [H(dH /dX)],—the ‘ force 
constant ’ /’,—as functions of current i, was carried out by two methods. 
The first of these was the conventional search coil and fluxmeter method, 
the usual precautions being taken to avoid the cutting of lines of flux 
by parts of the circuit other than the coil. The coil—the dimensions of 
which, 4mm long by 5mm diameter, were accurately measured—was 
mounted on a brass holder which allowed it to be rotated and its position 
adjusted and measured to 0-1 mm in three dimensions. The fluxmeter 
itself was calibrated by reversing a known current through a standard 
mutual inductance in series with it. 


Schematic Arrangement of Apparatus. 
M, Travelling microscope; A, Rectangular aperture 5 O, Slit source ; 
B, Sucksmith balance; D, Damping member; 8, Specimen; P, Pole 
pieces ; C, Calibrating device. 


This field calibration was carried out for five currents 13, 16, 20, 26 
and 40 amp which were selected after a rough preliminary measurement 
to give approximately equal increments of 1/H (see §3). It was made 
at intervals of about 5mm over a volume 2 cm high by 2 cm deep by 
lem in width over which the force constant was a maximum and 
approximately independent of position. The accuracy was limited 
to about 1% in H on account of the small size of search coil required. 

In the second calibration method F’; was measured directly as a function 
of position and current in terms of the force on a small piece of tantalum, 
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the susceptibility of which had previously been measured accurately by 
Hoare and Walling (1951) using the Gouy method. By integration of the 
FX curves, values were found for the variation with x of 4H?, 4H? 
being the difference between H? at the point 2 and at the centre of the 
calibrated region. These values were fitted to the directly measured 
curve of H? against X, thus allowing its accuracy to be improved. Mean 
force constants for specimens of any length lying within the limits of the 
calibrated region could then be calculated from this curve using the 
relation F,—1(H,2—H,?/X,—X,). A small further correction (up to 
about 2°) was applied to these values to allow for the effect on F’; of the 
differences in diameter between the specimens and the standard. This 
also was found experimentally. 


2.2. Balance 


A Sucksmith (1929) ring balance was constructed to measure the forces 
on the specimens. The ring was of beryllium copper strip 0-15 mm thick 
by 4 mm wide by 14 cm in diameter and had mirrors attached at positions 
of maximum sensitivity 50° above the horizontal diameter. From the 
bottom of the ring was suspended a silica rod 1 mm diameter by 40 cm 
long which ended, outside the constant force region, in a hook. The 
specimens were hung from this hook on a loop of fine gold wire. 

Vibration of the ring as a whole was prevented by attaching oil damping 
members to points half way up the ring at the ends of the horizontal 
diameter. The upper quadrants of the ring which carried the mirrors 
were still free to vibrate in their normal modes, but vibrations of this 
frequency were prevented from reaching the ring by mounting the metal 
frame which carried the balance and optical train on suitably designed 
sponge rubber pads. The vibrational amplitude was thus reduced by 
a factor of 1000. 

The deflection of the ring was measured with a travelling microscope 
in terms of the angular displacement of a beam of light reflected at the 
two mirrors. To get the required accuracy of force measurement the 
position of the beam had to be observed to an accuracy 0-002 mm or 
better. For this purpose a fine structure of about 0-01 mm in width 
was introduced into the image formed at the microscope. This was 
done by interposing a rectangular aperture in the path of the beam 
between a horizontal illuminated slit source and a focusing lens, both of 
which were on the opposite side of the balance from the microscope. 
The image of the slit, formed in the focal plane of the microscope, was 
thus superimposed with a diffraction pattern of fine black lines. It was 
found possible to set the cross wires on one of these lines to a precision of 
+0-001 mm. 


The balance was calibrated using a series of known wire weights. 
The calibration was found to be linear within 0-2°% over a range of 40 dyne, 
the sensitivity being about 4dyne/mm for a microscope to balance 
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distance of 2m. A further calibration could be made during an experiment 

by dropping a steel ball of known weight from a small electromagnet on 
to a pan attached to the balance. This electromagnet was normally 
switched off, the ball being held in contact with it, off the balance, by 
remanent magnetization. 


2.3. Ancillary Equipment 


The temperature of the specimen could be varied between 23° and 300°K 
by passing a stream of hydrogen or nitrogen gas from a reservoir of the 
boiling liquid through a heat exchanger coil wound around the tube which 
enclosed the specimen. Thermal insulation was provided by a Dewar 
flask which fitted between the poles of the magnet. In this particular 
investigation only two temperatures of measurement were employed, 
liquid nitrogen temperature 77°K and room temperature 298°. 

The specimen temperature was measured to an accuracy of 1°K using 
a calibrated copper—constantan thermocouple soldered to the outside 
of the heat exchanger tube level with the specimen. Preliminary 
calibration experiments determined the conditions under which the 
temperature recorded by this thermocouple was within 1°k of the 
temperature of a dummy specimen measured by a thermocouple in 
contact with it. 

To prevent condensation on the specimen at the low temperatures and 
to reduce temperature fluctuations of the balance the specimen and 
balance enclosures were evacuated. It was necessary to keep the balance 
temperature as constant as possible because of the temperature sensitivity 
of the deflection which amounted to about 1mm/°c. This was chiefly 
due to an additional elongation of the ring under its static load of about 
10 g as a result of the variation with temperature of the elastic modulus. 
Accordingly the room temperature was controlled to --0-1°c over the 
period of an experiment, using a multijunction thermocouple as a sensitive 
element. The thermal capacity of the balance enclosure and the insula- 
tion provided by the vacuum reduced still further the temperature 
fluctuation at the balance and, providing the balance was in equilibrium 
with the room before pumping down, the resulting zero drift was relatively 
small and could be allowed for. 


§ 3. OPERATING PROCEDURE 


The specimens used were for the most part of low susceptibility. 
They were made in the form of cylinders about 4mm diameter and 
1-5-1-:8 em in length. A hole was drilled along a diameter near one end 
of the cylinders to enable them to be hung from the balance. 

In the case of the brass specimens appropriate amounts of 99-995% 
pure copper and zinc containing about 10 p.p.m. hydrogen and 50 p.p.m. 
of oxygen, were melted in a high frequency induction furnace. The 
ingots obtained were sectioned and examined metallographically for 
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traces of segregation, samples also being taken from the top and bottom 
to check that the concentration was uniform. The individual specimens 
were machined from pieces cut from the centres of the ingots. 

After manufacture the specimens were weighed, after a preliminary 
acid treatment to remove any traces of iron contamination from the 
surface. They were then suspended from the balance in a measured 
position, small weights being added or subtracted at the balance to keep 
the total load approximately constant and avoid grossly distorting the 
ring. The optical system was then adjusted to obtain the required sharp 
and achromatic diffraction image at the microscope. Care was taken 
to avoid the mirrors acting as additional stops and introducing a subsidiary 
diffraction pattern. The balance lid, heat exchanger tube and Dewar 
flask were next put on and the magnet moved into position. The 
apparatus was then evacuated, first to a pressure of 1 mm Hg and then, 
after the specimen had reached the required temperature for the experi- 
ment, to 10-3 mm Hg. 

At the start of a run a number of readings of the zero position of the 
image were made as a function of time over about 2min. The magnet 
current was then switched on and adjusted to 13 amp, particular care 
being taken to avoid overshooting the current and introducing field 
errors due to hysteresis. Readings of the deflected position were taken 
as a function of time over the next 2 min, after which the magnet current 
was switched off, a small current being passed momentarily in the opposite 
direction to demagnetize. 

A further series of zero readings were taken for the next 2 min after 
which the magnet current was switched on at the next value (16 amp). 
In this way a series of readings of undeflected position of the zero were 
obtained as a function of time, interspersed by deflected readings 
corresponding to the five standard currents. At the end of the series 
the balance was recalibrated. 

By plotting the zero readings against time, interpolated zero readings 
were obtained corresponding to each deflected reading, and hence the 
actual deflections for each magnet current. After correcting these for 
the deflections due to the suspension, corresponding values of the apparent 
mass susceptibility (Xapp) for each current were calculated using the 
force constants, balance sensitivity and mass of the specimen. 

To correct Xapp for errors due to ferromagnetic impurities, which 
introduce a field dependent term, the method of Honda (1910) and 
Owen (1912) was used in which, providing the impurity is saturated, 


xX — X -[- a 

t app 

rue H,, 
o being the saturation magnetic moment of the impurities per g of 
specimen and #H/, the average field. 

Xtrue is thus found by plotting Xapp against 1/H, and extrapolating to 
Vibe) 
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§ 4. ResuLts 


As a check on the accuracy of the calibrations the susceptibility was 
first measured of 99-995°%, pure copper, silver, gold, aluminium and 
platinum. Table 1 gives the values found for these compared with 
apparently reliable values found by various other observers. Since our 
figures are relative to tantalum as 0-849 10-6 e.m.u./g the value found 


Table 1. The Susceptibility of Various Pure Metals 


Metal | Specimen | Temp. Xm &-m.u./g Xm found by other observers at room temperature 
Number °c xaos e.m.u./g x 108 
Cu i 22 | —0-084-+.. 0-002 | —0-086 (D) | —0-085 (HRe) | —0-085 (A) | —0-086 (HRo) 
2 22 | —0-084-+.0-002 
3 20 | —0-084+0-002 
—196 | —0-085+.0-002 
Ag 1 22 | —0-177+0-002 | —0-183 (D)]| —0-185 (V) —0-176 (A) | —0-182 (HRo) 
2 22 | —0-179+0-002 
3 22 | —0-177+0-002 
Au 24 | —0-140+ 0-002 | —0-142 (D)} —0-145 (V) —0-140 (A) | —0-142 (HRo) 
Al 22 | +0-615+0-005 | +0-608 (D)}| +0-62 (HRe) 
2 )-967--0-005 | -Lo- 
Pt 23 | +0-967+0-005 eS 0-971 (HW) 
Ee 
Ta 0-864 (D) 0-849 (HW) 
References (D) Donoghue, J. J., 1953, NAA-SR-117. 
(HRe) Hutchinson, T. 8., and Reekie, J., 1948, Phys. Rev., 78, 517. 
(V) Vogt, E., 1932, Ann. d. Physik, 14, 1. 
(A) Auer, H., Riedl, E., and Seeman, H. J., 1934, Zetts. Phys., 92, 291. 
(HW) Hoare, F. E., and Walling, J. C., (1951). 
(HRo) Henry, W. G., and Rogers, J. L. (1956). 


by Donoghue for tantalum is included. Agreement is generally satis- 
factory except in the case of platinum where Donoghue’s value is 
appreciably higher than our own. This may be connected with the 
particular sample used by Donoghue since he found appreciably different 
values for two apparently identical specimens. 

The silver results are interesting since they demonstrate what errors 
can be introduced by quite small amounts of impurities in ferromagnetic 
form and also the reliability of the Honda-~Owen method of correction. 
Figure 2 shows the apparent susceptibility as a function of 1/H for the 
specimens Ag], Ag2 and Ag3. These were prepared from a single rod of 
spectrographically pure silver supplied by Messrs. Johnson Matthey. 
Agl was machined directly from the rod and showed only a moderate 
degree of ferromagnetic impurity. Ag2 was Agl annealed in vacuo 
for 1 hr at 900°c. The slope of the curve, proportional to the saturation 
magnetic moment of the impurities, is greater for this specimen. Ag3 was 
remelted in vacuo and cooled slowly. Here the slope of the curve is 
greater still, the values of apparent susceptibility for the low fields being 
actually positive instead of negative as for pure Ag. 

In each case any surface contamination of the specimen was removed 
by boiling it in concentrated HCl while in the case of Ag2 the specimen 
was removed once from the apparatus to be reboiled in HCI and a second 
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time to have its surface removed by electropolishing. This treatment 
had no effect on the apparent susceptibility. 

Despite the great range of slopes the values of Xtrue taken from the 
intercepts of the curves at 1/H=0 coincide within the experimental 
error. The fact that the lines are straight indicates that the impurities 
were saturated at the lowest field (5500 oersted). Analysis of the samples 
after measurement showed in each case the principal ferromagnetic 
impurity to be 15 p.p.m. of iron. The presence of nickel and cobalt 
could not be detected, the limits of sensitivity being 1 and 5 p.p.m. 


Fig. 2 
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Apparent Susceptibility—1/H for Three Silver Specimens. 


respectively. The difference in the ferromagnetic behaviour of the three 
samples is presumably to be attributed to the varying degrees of precipita- 
tion of the iron content on going from Agl where the metal was chill 
cast and heavily cold worked to Ag3 where the metal was cooled slowly 
from the melting point. Metallographic examination of Ag3 at magni- 
fication of 1300 showed no indisputable signs of a precipitate which could 
have been responsible for the ferromagnetism though this is not necessarily 
surprising since, as shown by Reincke (1954) in the case of thin iron 
films permanent magnetization can occur at a thickness of only 11 A. 

The exceptionally high ferromagnetic impurity content of Ag2 and 
Ag3 characterized by the large slopes of the yapp-1/H curves did not occur 
for the other metal and alloy specimens. The slopes for these were 
roughly the same as for Agl1. 
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Table 2. The Mass and Atomic Susceptibilities of Brass 


og ote eS Soa A eee TE Se a ea A 


Concentration —Xm * 10° e.m.u./g axa 10~* e.m.u./g 

atomic % zine at atom at 
298°K 77°K 298°K 
0-084 5-41 5-34 
0-090 5-67 eo 
0-0925 5:89 5-89 
0-097 6-06 6-18 
0-109 6-90 6-97 
0-12] 7-81 75 
0-135 8-71 8-65 
0-139 8-75 8-9 


0-190 


Estimated accuracy of y,, is + 0-0015 x 10-8. 


Fig. 3 
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The variation of the mass x, and atomic x, susceptibilities of ie 


copper-—zine solution at 77° and 298°K are shown in fig. 3 and table ‘ 
i i 1X j rements, each 
Each point shown is the average of six separate measur 


corrected for ferromagnetic impurity. 
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Over the range of « solid solution from 0-32% zinc the alloys become 
increasingly diamagnetic the values of y,, ranging from —0-084 x 10-8 
for pure copper to —0-138 x 10-6 e.m.u./g for 32 atm % zine. The value 
for B brass, 47-6 atm % zinc, is lower again at —0-190x 10-6 e.m.u./g. 
Within the experimental error, estimated at about +0-015x 10-6, the 
values of the susceptibility at 77° and 298°K are the same. Also 
reproduced in fig. 3 are the values found by Henry and Rogers. Their 
values for the « range are systematically more diamagnetic than our 
own, the discrepancy at the higher concentrations being about 6% of 
our value. There seems no obvious explanation for this. 

In connection with another investigation, a single measurement was 
also made on a specimen in the analogous silver-cadmium system at 
30-5atm °% cadmium. The values found were —0-215 and 
—0-217+0-002 10-8 e.m.u./g for y,, at 77° and 295°K respectively, 
compared with —0-177 at 295°K for pure silver. In atomic units this 
represents a decrease of 4:6 x 10-6 e.m.u./g atom between 0 and 30 atm % 
compared with a decrease of 3-2x10-%e.m.u./g atom over the same 
concentration interval for copper—zinc. 


§ 5. Discussion 


The total magnetic susceptibility of an alloy is usually discussed 
theoretically in terms of separate and distinct contributions made by 
various groups of electrons. Distinction may be made in the case, say, 
of the copper—zine system between contributions from electrons 


(a) up to and including the (3d)? shell of the Cu* ions ; 
(>) up to and including the (3d)? shell of the Zn*+* ions ; 
(c) distributed around the Zn** ions to screen their single excess charge. 


These three groups give rise to an orbital diamagnetism of the type 
first treated classically by Langevin. 

The 4s electrons from the copper and zinc atoms which enter the 
conduction band give 


(d) a paramagnetic effect (Pauli paramagnetism) as a result of the 
partial nett alignment of their spins. This contribution, is at 0°K 
proportional to n(¢) the density of states at the Fermi level. It may 
be further enhanced by exchange interaction effects which, if positive, 
tend to increase any existing inequality in the numbers of + and — spins. 


(e) A diamagnetic effect due to the curvature, imposed by the field, 
on the paths of these electrons as they move through the lattice. 


(In this treatment any orbital contribution to the diamagnetism from 
the copper 4s and second zinc 4s electrons has been ignored). 

No one of these contributions can be calculated exactly (although some 
separation is possible in practice using information from other types of 
experiment) so that, in cases such as copper-zine where the contributions 
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are believed to be all of the same order of magnitude, no great accuracy 
can be expected in the final theoretical value for the absolute susceptibility 
at any particular concentration. 

Variations in susceptibility with concentration can be expected to arise 
in varying degrees from the following : 


(i) The replacement of Cut ions by smaller Zn++ ions having a smaller 
diamagnetism. 

(ii) The change in the individual contributions of Cu+ ions and Zn++ 
ions as, with the expansion of the lattice on alloying, the ions move further 
apart and their 3d shells (responsible for about 85°% of the diamagnetism) 
interact less. 

(iii) The total diamagnetic contribution from the screening 4s electrons 
which will increase with zine concentration. 

(iv) Changes in the shape of the n(#) curve and hence in the Fermi level 
and n(¢) for a given conduction electron concentration due to (a) the 
overall expansion of the copper lattice and (b) by local perturbations of 
the periodic potential. 

(v) Filling of the band as the concentration of conduction electrons 
increases with zinc concentration. 

(vi) Variation with lattice parameter and wave function symmetry of 
exchange interaction effects. 

(vii) Variation of the conduction electron diamagnetism as the first 
Brillouin zone is filled. 


The problem of deciding which of these are the major effects and of 
estimating their size is no easy one. 

In their treatment of this problem Henry and Rogers follow Friedel 
(1952) in assuming that one important term governing the slope of the 
susceptibility concentration curves for their copper base alloys is the 
resultant change in orbital diamagnetism of (i) and (iii). The distribution 
of screening charge is assumed to be the same as that of a (4s) 
electron on an isolated Zn++ ion. Stoner’s (1929) formula for ionic 
diamagnetism is used with the charge distribution for Cut, Zn*++ and 
Zn(4s) calculated by the self consistent field method. This orbital term 
accounts for about 31—4 of the observed negative slope found by Henry 
and Rogers for copper-zinc, copper-gallium, copper-germanium and 
copper-arsenic, but provides a small positive contribution in the case 
of copper-aluminium where the nett slope is also negative. 

The major contribution to the slope is assumed to come from a decrease 
in the Pauli paramagnetism (governed by (iv) and (v)), a change in the 
shape of n(Z) solely with lattice spacing being assumed which is the same 
as that for free electrons. For the n(#) curve for copper, Henry and 
Rogers take the curve postulated by Jones (1937 b) to account for the 
« and f solubility ranges in copper—zine and similar alloys: in this the 
Fermi level is close to a maximum beyond which n(#) decreases roughly 
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linearly with electron/atom ratio. Henry and Rogers assume that ¢ 
is actually past the maximum so that addition of further electrons on 
adding zinc (or on Friedel’s model the subtraction of states from the bottom 
of the band) produces a linear decrease in n(¢) with concentration. This 
term is the one in their view primarily responsible for the decrease in 
susceptibility; other possible effects such as (ii), (iv 6) and (vi) are 
ignored. 

The agreement cited by Henry and Rogers between theory and 
experiment for the slopes of their five solid solutions is apparently 
excellent (cols. 2 and 3, table 3) and at first sight the assumptions made 
seem to be amply justified. On closer examination, however, there 
seems considerable doubt as to whether the basic assumption of a negative 
slope of n(H#) for copper is in fact correct. This point will now be 
examined in some detail. 


Table 3. Calculated and Observed Slopes Derived from Henry and Rogers 


(1) alloy system; (2) slopes observed by Henry and Rogers; (3) quoted 
‘calculated’ slopes; (4) ‘ calculated’ slopes for zero slope of n(#) for Cu. 
The slopes of the curves of susceptibility against concentration are in units 
of 10-* e.m.u./g atm. 


In his calculation of n(H) Jones proposed an expression for EH as a 
function of k,, k, and k, to hold in the immediate neighbourhood of the 
k,-direction (taken as one of the 8 equivalent (111) directions). This 
was arrived at purely intuitively as being the simplest expression which 
satisfied certain necessary conditions at the zone boundary and centre. 
n(H) was obtained from this by the further arbitrary approximations of 
assuming the E(k) relation to hold over a circular cone of solid angle 
37 in k space and of replacing the actual Fermi surface by the surfaces 
associated with 8 such cones, one for each octahedral face of the first 
Brillouin zone. Thus Jones’ n(#) curve cannot be regarded as sufficiently 
firmly based for a numerical calculation and the quantitative agreement 
which he finds for the « and f solubility limits may well be fortuitous. 

What evidence there is for the actual n(H) curve for copper suggests 
in fact that the Fermi level is well on the low energy side of a maximum 
so that n(E) increases initially with #, or in other words that the Fermi 
surface has not touched the zone boundary. 
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Thus for example Cauchois (1953) in an analysis of the copper K and 
L x-ray absorption spectra finds an absorption maximum for the L III 
band at about 1-0 ev beyond the absorption edge. Since this corresponds 
to transitions to states of primarily s symmetry in the band and the 
s character decreases and p increases with increase of energy, this must 
imply that n(#) itself is increasing with EH up to this point. 

A maximum in the & absorption band at 2-7 ev indicates that n(#) 
probably has a maximum, with a fair admixture of p character, at about 
this energy. Beyond 2-7 ev n(Z) is believed to drop to a minimum before 
rising again, corresponding to overlap of the first and second zones. 
This experimentally deduced behaviour agrees well qualitatively with 
the theoretical n(#) curve for copper derived by Rudberg and Slater 
(1936) from Krutter’s (1935) quantum mechanical calculations of EZ as 
a function of k. Rudberg and Slater predict that n(Z) rises to a maximum 
about 1-5ev beyond ¢ with a minimum at about 3-5 ev and a second 
maximum at 5 evt. 

If n(¢) in the alloy depended, apart from a slight correction for changes 
in lattice parameter, simply on filling the copper conduction band, then 
it seems most likely that it should increase with concentration and not 
decrease as required by Henry and Rogers. The situation is not as 
straightforward as this however since as shown by Bearden and Friedman 
(1940) the edges of the copper and zinc K absorption bands for copper—zine 
alloys do not shift appreciably relative to the fine structure of the bands 
as the composition is increased. This would seem to imply that in fact 
little movement in €¢ occurs relative to the first maximum in n(#) for 
copper the actual slope of n(#) not being the factor principally determining 
n(C) in the alloys. 

Fortunately direct information about n(¢) is now available from 
measurements by Smith (unpublished work in progress) of the electronic 
specific heat, Cz, in the copper-zine system. Smith finds y(=C;,/7) 
which is proportional to n(Z) to show a significant increase with zinc 
content from a value of 1:65-0:05 x 10~ cal/g/atm/deg? for pure copper 
to a value of 2-00+0-05 x 10~4 cal/g/atm/deg? for 32 atm % zinc, a 
behaviour quite at variance with that on which Henry and Rogers base 
their calculations. Thus it seems likely that »(¢) and hence the para- 
magnetic term in the total susceptibility for copper-zinc should increase 
initially with zinc content. Even if it only remains constant however 
the agreement obtained on Henry and Roger's treatment is completely 
destroyed (see column 4, table 3). t 

We must accordingly anticipate a diamagnetic effect additional to those 
which Henry and Rogers have considered and one, moreover, probably 
ee ee 

+ A similar calculation of E(k) made by Howarth (1953) for the (111) and 


(200) directions in phase space also predicts that n(#) should increase initially, 
the Fermi surface being 1-2 ev below the energy at which the surface first 


touches the zone boundary. | 
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related in some way to the valency of the solute so as to give the observed 
increase, with solute atom valency, in the negative slopes of the 
susceptibility/concentration curves. Such a contribution could come 
from the conduction electrons since, as Krutter shows, these cannot be 
considered as free, the energy surfaces being far from spherical. 

In the case of bismuth and the y brasses Jones (1934, 1937 b) has shown 
that the relatively much larger diamagnetism can be accounted for by 
supposing that the energy surfaces of the electrons overlapping into the 
second zone are highly ellipsoidal. 

Unfortunately it is difficult to make a reliable estimate of the conduction 
electron diamagnetism for the « solutions. Until this can be done 
however further speculation as to the reason for the drop in susceptibility 
would not appear profitable. 


§ 6. CONCLUSION 


The decrease in magnetic susceptibility with concentration found by 
Henry and Rogers for the copper—zinc, copper—arsenic, copper—gallium 
and copper-germanium solid solutions and by ourselves for copper—zine 
does not appear to be attributable to a decrease in Pauli paramagnetism 
consequent on filling an electron band derived from the copper conduction 
band. Neither, allowing for the inaccuracies in the calculation, does it 
seem likely that the additional diamagnetism of the solute ions—assumed 
singly ionized—can be the major factor. No satisfactory explanation 
can be suggested at present for the decrease. It may, however, be due 
at least in part to an enhancement of the conduction electron 
diamagnetism. 
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SUMMARY 


The energy of a gas of fermions of spin }, interacting with a hard sphere 
repulsion, is obtained by means of Brueckner’s approach to the many- 
body problem. The single particle potential is computed in the effective 
mass approximation, and the value of the maximum density for which 
this approximation is valid, is derived. By means of a modification of the 
effective mass approximation, in which the effective mass becomes spin- 
dependent, it is shown that above a certain critical density, the gas 
becomes ferromagnetic. 


§ 1. INTRODUCTION 


THE method which has been developed by Brueckner and his co-workers 
(Brueckner and Levinson 1955) for the derivation of the properties of a 
nucleus in its ground state and in states of not too high an excitation is 
applicable to any many-body problem in which the particles interact by 
means of two-body forces only, even when these are very strong, and in 
which it is sufficient to treat two particle correlations precisely without 
a precise treatment of correlations of three or more particles. In parti- 
cular, the method should be applicable to the theory of liquid *He. 
Since we are dealing with a system obeying Fermi—Dirac statistics, the 
arguments of Brueckner (1955 b) and of Bethe (1956) concerning the 
smallness of the cluster terms should hold, although the three particle 
cluster must, of course, be evaluated in order to estimate the error in the 
energy. In fact the theory is applicable to this problem in a much 
cleaner fashion than to the nuclear one, since in the former case, the only 
source of error arises from the cluster term, whereas the terms inversely 
proportional to the number of particles constitute no source of error 
whatever, and there is strict conservation of momentum throughout the 
problem. Furthermore, in the nuclear case, a precise treatment involves 
the guessing and making self consistent of a single particle potential 
matrix which is, in general, not diagonal in any simple representation 
like position or momentum, and the definition of matrix elements of various 
operators with respect to Slater determinants which are not built up from 
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plane waves, but from functions whose form must be changed in the 
course of the guessing in the self consistent calculation, Although 
Brueckner (1955 a) and Skyrme (1956) have carried out calculations 
taking nuclear surface effects into account, no calculation which takes 
into account the full complication of the finite nucleus case has as yet 
been carried out. 

In this paper, only results for the idealized case of a gas of hard spheres 
with spin } and Fermi—Dirac statistics are discussed. Those for liquid 
*He will be presented in a subsequent paper. In § 2, an expression for 
the single particle potential in the effective mass approximation (e.m.a.) 
is written down. In § 3, the value of the effective mass is derived in the 
low density limit, and an expansion of the energy as a function of density 
is given. In§ 4, the effective mass and the energy are calculated at higher 
densities, and the limit of validity of the e.m.a. is clearly set forth. In § 5, 
the spin-dependent effective mass approximation (s.d.e.m.a.) is introduced 
and the single particle potentials are obtained. The effective masses and 
energy in the low density limit are given in § 6, while in § 7 is calculated 
the value of the density above which the antiferromagnetic configuration 
is no longer the least energetic. In § 8, the characteristics of the gas in 
the ferromagnetic region of densities at absolute zero are given, as well as 
a model which leads to an overestimate of the Curie temperature. 


§ 2. SINGLE PaRTICLE POTENTIAL IN THE EFFECTIVE Mass 
APPROXIMATION 


Let the interaction between a pair of particles be denoted by v;,;, while 
the corresponding reaction matrix is denoted by G,;.. The (non- 
symmetrized) matrix elements of v will be denoted by 


(n,'n;' pyran) = | b(n,’ r,)b*(n,'; 1 OR 9 CC r yb (n;; r;) dr, dr, (1) 


where (n,; r,;), etc. denote the solutions of the Schrodinger equation with 
the single particle potential. The (n,;r,;) are plane waves in our case, 
since the single particle potential V is diagonal in the momentum re- 
presentation for an infinite medium. In writing the matrix element 
(n,'n;'|v,;\n,n,), it is understood that all particles not explicitly written 
down are in the same state as in the chosen configuration (Bethe 1956). 
(When a definite particle is said to be in a definite one particle state, a 
convention like picking out the main diagonal of a Slater determinent is 
being used.) Of course, no single particle state is ever multiply occupied. 


Tn general, let Derg, OF Uae lop aeons 4 (2) 
for a triplet or singlet pair, so that 
V:;=4(80,+0,) +2(0,—%,)(5; 0; sme «a SeuntOD 


Z2E2 
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Then the reaction matrix for two particles with spins pointing in the 
same direction is given by the integral equation : + 
(k,’ky |@ gr] ik) = (Ie, [51 — Pay), 1kk;) 
—(2r)-# f dk,"(k,'k,’ |b —Py)odk,"k;”) 
2 We 
Ei(k,")+E(k;")—E(k,)—H(k;) 
where P,,, is the Majorana exchange operator. That for spins pointing in 
opposite directions is given by 
(k,’k,’|Gsir71k:k;)=(k,’k,’|3(1—P.y)v, |k,k;) +(k,k,’ 21 + Pa)0,|kk;) 
eh eg | Bk,"[(k,'k,'[3(1—P3,)¥,1k;"k;”) + (k,'k;’ [21 +P 3,)0,|k,"k,")] 
dz: 
X tL PCL He Pic Ie) 
E\(k,")+ E(k;")—E(k,)—(k;) 
In eqns. (4) and (5), k,;’+k,’=k,”+k,;’=k;+k,, the single particle states 
being denoted by the wave number k of their plane waves. The integral 
over k,” is such that neither k,” nor k;’(=k,+k,—k,”) is occupied in the 
chosen configuration, while k,;, k; are assumed to be the wave numbers of 
particles 7 and j in the chosen configuration. H(k) is the single particle 
energy for wave vector k, and we have assumed that the excitation energy 
of (k,"k,"k,°. . .) with respect to (k,;k,;k,°. . .) can be expressed entirely in 
terms of one particle energies, as in eqns. (3)—(8) of Bethe 1956. Goldstone 
(see Bethe 1956) has pointed out that this is not in general true, but we 
shall make use of it as a simplifying assumption throughout this paper. 
We take v,=v,=v, and we assume that v is an ordinary (Wigner) force. 
The single particle potential in the chosen configuration is given by the 
integral over the diagonal elements of the reaction matrix : 


V(k)=(2ry8 [ Pk, (kiks|2G 47+ Gisrrlkk,) tee cee), 40) 


where the integrals run over all k; in the chosen configuration. The 
factor 2 in front of G,,;; in (6) arises from the possibility of having both 
spins up or both spins down. 

In the e.m.a. one assumes that V(k;) can be written down as 


(k,"k;”|Gisrlkk;), iy (4) 


(k,"k,"|@sr7|k,k,). By ks (5) 


V(k,) == V (O)}asb 02 sce. ei eC) 
where an effective mass M* is defined by the equation 
1 1 


since H(k,) is then given by 
E(k,)=V(0)+k/2M*. 2 oh an en Eo) 
Let us now take for the chosen configuration that in which all states 
with wave numbers from 0 to a certain maximum (Fermi) wave number 


+ In this paper, only neutral fermions are considered. There are 4 instead of 
2 different reaction matrices if two charges (isotopic spin orientations) are also 
allowed, 
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ky» are occupied by two particles of opposite spin, while all single particle 
states with wave numbers k>k, are empty. Then the integral in eqns. 
(4) and (5) runs only over those k,” for which k;" >ky and |k,;”—k,—k,|>ky. 
Let us simplify the problem by neglecting this restriction and integrating 
over all k,” instead. The effects of this replacement will be discussed 
below. 

The reaction matrix is now that for the interaction of two particles of 
mass M* each. For a hard sphere interaction of radius a, i.e. V;=+0 
for r;;<a and 0 for r;;>a, the solution is well known to be 

Ne? ee t Stine ST issjalkis'@) 

ee weigh) a6 OUT te) 

where k,;=3|k;—k,|, &;,’=4|k,’—k,’'|, and @ is the angle between k,—k, 
and k;’—k,’. For the diagonal matrix elements, we may write 


P,(cos 8) (10) 


— 47; | 
(k;k,|G,;|k;k;)= Us kez, 2» (20+ 1) tan 8; : fe 3 . (11) 
where 
tan 8,=(—1)"1 Siraj (K;,0)/J_y ajo (hiya) 2 6. (12) 


6, is the phase shift for orbital angular momentum / for a hard sphere of 
radius a. (Eqn. (11) applies to any central potential.) 
Hence, eqn. (6) becomes, for a central potential writing V(k,;) for V(k,) : 


ee ne | rie Tea 21 1h, an (13) 
(e=— aaa)” [fas gag, ELF Vii tam Bik) — 
where 

ete Ly ee ge ew ae ee (LO) 
6, is the phase shift for orbital angular momentum / caused by vt. By 
means of a straightforward transformation of variables, (and writing x 
for k,;) this may be written as : 

8 1/2ky+kp) 
Viki=— sy | 

for k;<ky, while for k;>k,, we get 


dis - uP, (x) > (20-1) fi ile) = (LS) 


“=U 


paw hla )¥ (A+VS, 8:(z) (15) 
Ve — 75 | ssn Gt (© ,O,(%).  . 
2 for 0<a<4h(kp—h) 
H a = SS SOV aN 2 ; 
eg tte 2) kp b<ac(hhp th). . . . (06) 


ka 


+ Riesenfeld and Watson have shown that the integral eqns. (4) and (5) should, 
in a more rigorous treatment, contain a matrix element of the reaction matrix 
as a term in the denominator of the propagator, and that therefore the factor 
tan 8, in (13) should be replaced by 6;. This treatment corresponds to en- 
closing the entire system in a very large but finite box and imposing the suitable 
boundary condition of wave functions going to zero on the surface of the box. 
In what follows, tan 8, will be systematically replaced by 9. The author is 
very much indebted to Dr. W. Wada for bringing this point to his attention. 
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Oe ee $(k— lp) <2<4 (kooks os lonrla) 


In the chosen configuration, we do not need eqn. (15), since k;<ky, but in 
5, this case will be required and it has therefore been derived here. 
For the case of a hard sphere interaction, (14) becomes 


8 pUAkg+kp) 
Vieje= a de. aP (x) 52+ hi 
1 va ce. Se ents 
as fe ) J _14/2(%@) ae 


After this paper had been written, the author read a preprint of a paper 
by Bethe and Goldstone, who criticize the computation of the effective 
mass by considering only the range 0<k<k,. The reasoning used in the 
present paper is based on the arguments at the end of § ITI of Bethe (1956). 
Since V(k) is not a well defined quantity for k>k,, no attempt has been 
made in this paper to make it self consistent for k>k,. Bethe and 
Goldstone find that their criteria for the validity of the e.m.a. in which 
the behaviour of V(k) for k>k, (suitably defined in a certain manner) is 
also investigated, are such to allow the e.m.a. to hold for very high 
wave numbers k. However, the ‘ effective mass’ for high & is not the 
same as the ‘ effective mass ’ for low & in general, so that this cannot be 
spoken of as validity of the e.m.a. to high k, except in a very rough manner. 
The functional dependence of eqn. (7) in the range 0<k<k, is still a 
necessary condition for the e.m.a. and the conditions derived in § 4 of 
this paper cannot be relaxed. It is true, however, that M* is increased 
by not neglecting the Pauli principle in the intermediate states of eqns. 
(4) and (5), and that this will therefore lead to an increased maximum 
density for the range of validity of the e.m.a. 


§ 3. Errective Mass anp Tota ENERGY IN THE Low Densrry Limrr 
The chosen configuration, in which each single particle state with wave 

number /<k, contains two particles, while all those with k>k, are empty, 

corresponds to a gas at the absolute zero of temperature and of density 
N 1 

Sas ee ae ee ad ee on a 
OF 8a 9) 


where N is the total number of particles and Q is the volume. Hence- 
forth, the hard sphere radius will be taken as the unit of length, so that 
a=1, (h=1 throughout this paper.) For sufficiently low density, ice. 
low kp, and k<k,, it is possible to drop all J values except the first few in 
eqn. (18), and to expand the remaining Bessel functions in a power series. 
Let only /=0, 1, 2 be retained, and replace 


> (2/+1)f, arctan | (- 1) aaa Tivol jit, |by —x — 3x34 at (20) 


1=0)1,2 hire 1/2 wv) 
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then there results 
anM*V(k jes Ens : i 
Siernel es” +(gbeP+ 7 ky ) 


19 
—(sag be" a rosa te” Pas igo te *k) verge) 


Now the validity of eqn. (11) (or of its Riesenfeld~Watson modification) 
for the reaction matrix depends on eqn. (7) being satisfied. If this is 
not so, then the procedure is not self consistent, and a more complicated 
assumption for the k- dependence of V(k) must be made. Equation (21) 
shows that in the limit of low densities (k,a<1), eqn. (7) is indeed well 
satisfied. To the extent that the terms of degree 7 in ky, k(k<k,) on the 
right-hand side of (21) can be neglected, the effective mass approximation 
is now completely justified. Hence, 

aM*V(k 1 

= =(5h+ g a le® \+q ptr... (22) 
Define c= M*/M and z=M/M*. Also define U(k)=7MV(k)/8. Then it 
follows from definition (8) of effective mass that 


16 
7=1+ Fa UH) U0) \ . . . . . . (23) 
16 
or o=1— To U(h)—U)}. < heh Baap tomers (24) 
Hence, 
ee Se aa, ple et 6 8) 
7 


This means that for small densities, the effective mass is slightly less than 
the real mass. It is interesting to note that if the hard sphere interaction 
had the p-wave part (J=1) suppressed, the result to this order of approxi- 


mation would have been 
Gates tre Ce eee a te tee er 


while if the s-wave part (/=0) were suppressed, eqn. (25) would still have 
been obtained. Hence, the quantity 1—o is different from zero at low 
densities chiefly because of p-wave interaction. The d-wave does not 
contribute to 1—o in the approximation (22), although in the approxi- 
mation (21), it contributes to the k* term in V(k) which causes the e.m.a. 
to break down. (The limit of validity of the approximation is discussed 
in § 4.) 

The total energy, or the energy per particle E, may now be derived in 
the low density limit. 


ky 
| * 12 db{k?/2M+4V(b)} 


D | tik 
Q 


where the factor } in front of V(x) serves to prevent the potential energy 


om) 
“ 


E— 
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between each pair of particles from being counted twice. Equations (27), 
(7), (8) lead to 


E—1V(0) pas (1 tS) Berar ae eC) 
It should be noted that the leading term dropped in (25) is of order ky? 
| o=1—~ ky +O(kp) Med ae ek, emo) 

so that (28), (25), and (22) ae the oe 
MEA hy? ky?+ ro hp? + gghe® +O(kp’). . - (28) 


The neglect of iy Pauli sere in intermediate states in the integrals 
in eqns. (4) and (5) should obviously be least serious in the low density 
limit. It turns out that if this approximation were not made, the first 
two terms of eqn. (29), proportional to n?/3 and n, would not be changed. 
This calculation will be contained in a subsequent paper. 


§ 4, Limit of VALIDITY OF THE EFFECTIVE Mass APPROXIMATION 


For k, approximately equal to 1-9, the e.m.a. breaks down because 
o/ky?{U(kp)— U(0)} becomes so large (viz. >7/16) that o passes through 
0 and becomes negative, according to eqn. (24) with k=k,. Clearly, 
o> 0 is a necessary condition for validity, since if <0, eqns. (18) and (28) 
lead to a negative total energy, arising from an interaction which is 
entirely repulsive. Since k,=1-9 corresponds to a density of hard spheres 
below the density of close packing, the range k,>1-9 is not to be con- 
sidered as unphysical. This absurd result simply means that for k, >~1-9, 
eqn. (7) does not hold at all for V(k), O0<k<ky, and that 

V(k)= V (0) OK? +b k* ee et ae (80) 
or some other form for V(/) must be guessed, with b, 6’ then being adjusted 
in place of /* in the course of the self consistent calculation. 

Even if o>0, it is still necessary that the o obtained from eqn. (23 6) be 
roughly independent of k for 0<k<k,. This condition is well satisfied 
throughout most of the range of k, for which co remains >0, as table 1 
shows. 

Figures 1 and 2 give the effective to actual mass ratio o and the value 
of ME (where 4=1, a=1) for various values of k,. At the lower end of 
the range eqns. (25) and (29) are reasonably accurate. The results were 
calculated by replacing 


> (2l+1)f, arctan [ (cy pee 


J _ 1-1/0 x) 
by an odd nalyneiial in w of degree 9 which gives a good fit over the 
whole range of «in question. Its coefficients were not obtained from any 
power series expansion of the are tangent of the ratio of the two Bessel 
functions. 
Inspection of table 1 shows that for k»,<1-6, the effective mass approxi- 
mation for hard spheres is reasonably self consistent. For ky21-9, it 
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Ratio of the effective mass to the actual mass o=M/*/M as a function of the 
Fermi wave number ky. 


Table 1. 


and k=4k, near the limit of validity of e.m.a. 


kp k=kp 


1-4 0-490 
1-6 0-317 
1-8 0-116 


Values of o(=M*/M) as obtained from eqn. (24) with k=k, 
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is absurd, since o becomes negative. For a hard sphere core with an 
attractive potential outside it, as in *He, the effective mass approxi- 
mation will be valid to a somewhat higher value of k,a than for the 
present case (a=radius of hard core part of the potential). 


Fig. 2 


2.0 
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Ma°E/h* (E=energy per particle) as a function of the Fermi wave number ky 


§ 5. Spin DerenpDENT Errective Mass APPROXIMATION 
The chosen configuration of § 2 is now replaced by another one, in which 
each single particles state with wave number k<k,,_, is occupied by two 
particles of opposite spin; each state with kiy-y<k<kyy.) occupied by 
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one particle, whose spin is aligned in the same direction for all singly 
occupied states; while all states with k >kiy4) are empty. We shall 
refer to the direction of the magnetic moments of the particles in the 
singly occupied states as the ‘ up ’ or ‘ plus’ direction, and to the opposite 
direction as the ‘ down’ or ‘ minus’ direction. 

The energy of the new chosen cenfiguration will now be worked out in 
the spin dependent effective mass approximation (s.d.e.m.a.) in which 
the single particle energy of an ‘ up ’ particle (one whose magnetic moment 
points up) is given by 


E(k)=k?/2M*(+)+V(-+ ; 0) Smee MOL) 
while the energy of a ‘ down ’ particle is 
H_(k)=k?/2M*(—)+V(— 0). Sw eS SC(81 8) 


M*(+-) and M*(—) being the two effective masses. 
Equation (13) (with the Riesenfeld-Watson modification) must now be 
generalized to the following form : 


rhA(+) 6 Ao odd 
V(t sk)=—Qny2| || de, rag D2) Silks) 


Lh ke(—) ( Pl. ie (QI-21)0(E,,) * (32 a) 
ot ie > i 


l 


' ki —) : 47 odd oe eh 
V(—; &,)=—(20)3 E | | Oks TE, > 2021+ 1) bulk) 


kay all 
—(2n)-*| pe | | ak, ey Ee (ene Ot) 
0 vs sti eis 

M,* is the harmonic mean 2J/*(+- YM (—)/ {ME *(+-) + M(— )} of 1*(+) 
and M*(—). Itis for the second term on the right-hand side of eqn. (32 a) 
that eqn. (15) was derived, since in (32a), O<k,;<hy), 80 that h, exceeds 
ky) for some of its range. It then follows from eqns. (14) and (15) that 
V(+-) and V(—) may be written as follows : 


If k< ky); 1/(k-+k Peal 
8 eth 4) nie 

: 5 od ee. P.(+ ; #) > -2(21+-1) 3,@) 
/2(k +key all : 

Ps 4 oe Bes 2) YU+1) (0). (83) 
TLL ¢ 0) 

Lf k>ky_) 

8 71/2(k-+k p+) care ae 
ee . a 2 P(-+; k)> 2(21-+1) 5,(x) 
V(+; Deer) dx « P(+ ; 2, ( i 
aM 9* J 1/2(k—ky(—)) 7 


where 
Plt ; v)=2 when O0<a<B(ky4)—h) 
2 (Lk—ax)*}/ka 
when Hhkp.)—h)<a<(thiyay th) . - - + (35) 
and Q,(— ; x) is obtained from @,(x) of (17) by replacing ky by Ky_). 
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If n(-++) and n(—) denote the densities (number of particles per unit 
volume) of ‘ up’ and ‘ down’ particles respectively, then 
N(+) 1 
DE) =O = gar 


(36) 
and 


1 
— 6n2 (kaa thr’): . . ° . . (37) 


§ 6. ErrectivE Masses anp ENERGY IN THE Low Densrry Limit FoR 
S.D.E.M.A. 


Let the phase shifts in (33), (34) be given by their values for a hard- 
sphere interaction (of unit radius), and expand up to and including terms 
in 2°. Then 


odd 
> 2(2h-- 1) 3;(f)==—= 2a ee ee oo a, 
1 
and 
all 
(Ziad 6,(¢) == (ae) se Petes 6) 
Then eqns. (33), (34) give 
8 I. 1 it 
Vilziss i) a (3 kaze t+ 50 “ney + 48 key h? ) 
8 i il aoe 


Define o(—)=1/7(—) = M*(—)/M and o(+)=1/7(+)=M*(+)/M. Then 
eqns. (31), aes ae Say two simultaneous equations 


a ee ~ (1-re( (+)/o())ka=F=1—o(+) me eg Cee), 


which have ae sees 


3 3 
2hra thy, 
37 


In eqn. (41), O(k,°) means that the leading terms dropped in deriving (39) 
and (40) are of degree 5 in ky), ky_) together. Equation (41) shows 
that if n(+-)>n(—), then o(+-)<o(—), i. e. the effective mass of a particle 
with magnetic moment in the direction of the majority is lower than the 
‘minority effective mass ’. 
If E is the average energy per particle, then 
Fe) Mere) 
ME= ef ee ae (HI) FAVE 5B) mf er. 
0 
NY (=e ‘ 
= [ k2 dk [(k?/2M)-+4V(— my uedk. . (42) 


o,=1— =O (po). ae eee: 
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Equations (42), (36), (39), (41) then yield the following result : 
mew 2 Rear thay?’ , 2 ker ekm_? 
10 ky pthky_) 37 ky tke 
J 1 
ae yy Pein a 2k the») + (hers Ske) + ebm 4)°)} 


a 
Me acaenes Rs COR) er = ee, (48) 


which of course reduces to eqn. (29) when ky 4)=hy_y= kp. 

The first term in (43), which is quadratic in the k,’s, is the contribution 
of the kinetic energies of the individual particles. The next term, which 
is cubic in the k,’s, arises from the term in (39) which is cubic in the ky’s, 
and which comes entirely from the s-wave interaction of the hard spheres. 
By eqn. (37), a fixed density of particles implies that ky, ,)3+k n(—)° 18 fixed. 
At sufficiently low densities, the first term in (43) dominates the others, 
and for fixed ky, 3+kp_)8, it is a minimum if ky.)=ky,_). Hence, at 
low densities, the antiferromagnetic configuration, for which as many 
spins (and therefore, magnetic moments) are in the ‘ down ’ as in the ‘ up’ 
direction, has the lowest energy, and is the configuration in which the gas is 
in equilibrium at temperature 7=0. For an ideal Fermi—Dirac gas, 
this is, of course, true for arbitrary density. 

Atasomewhat higher density, the second term of (43) becomes important. 
Its effect is to favour inequality of ky.) and k,_;. The physical inter- 
pretation of this is quite simple. This term is the leading term (in powers 
of ky), ky_)) of the energy due to s-wave interaction of the hard spheres, 
which at not too high densities is certainly more important than higher 
values of /, and it can exist only between a pair of particles with opposite 
spins, not between a pair with the same spin, in accordance with the 
Pauli principle. Hence, this interaction tends to line up spins (and so, 
magnetic moments) of near by particles and to produce ferromagnetism. 

Now, fix ky? +hy_?=2a, and minimize $(y)=ME(ky.>+ky_)*) as 
a function of y=}(ky4)3—hky_)?). We shall use (43) for WE, neglecting 
O(k,’), and assuming that |y|<« for the range of y in question. (Indeed, 
we know a priori that at low enough densities, y—0 minimizes ¢(y), and 
we wish to investigate at what density this ceases to be the only relative 
extremum.) For this purpose, it is necessary to expand :j(y) in powers of 
y?, up to and including y*. The result is 


3 5/8 By? 5 yt Lagi eg Ie 109 
oy)= po (1455+ aaa ghrieges (Coc Uh tiage eee Vitis Sg 
2 
33 (0° —y"). tc ppt Lee, (22) 
For « low enough, dys(y)/dy has only one real root, viz. y=0, corresponding 


to the above mentioned minimum. For higher %» dyb(y)/dy has three, 
real roots, y=0 corresponding to a local maximum in 2(y) (instability of 
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the antiferromagnetic configuration at 7'=0) and two minima y= +Yo, 


where 
27 Oe V5} 2 9 4/38 
i ae be sey jit Sp SMe eee, thee ED 
Yo x 5) (1 aes arpa ee ) ( ) 


yy becomes real as «1/8 becomes greater than ~3-0 which, by § 4, is far 
beyond the range of validity of the e.m.a. In § 7, the minimum density 
for ferromagnetism at 7'=0 is calculated without the low density approxi- 
mation (which (45) shows to be inapplicable), and it turns out that the 
limiting density is that for which k,—1-674 rather than the above estimate 
k,»=3-0. This lower value is within the range of validity of the e.m.a., 
although table 1 shows that the e.m.a. is beginning to fail at a somewhat 
higher density. 


§ 7. Mintmum DENSITY FOR FERROMAGNETISM AT 7’=0 


Define 
M 
Ut; = 7-V 
Assume that ky.) and ky_) are approximately equal, with ky4)>kp_); 
so that 0<D<C®* where 


(ESR) on. ee 


ky p=OLD. . : . . . . . . (47) 
Define 
D _ kp kro? 


4= = c Py TER! ays Ge tas ehh ties 
C8 kp) +key? 


(48) 
The average energy per particle is to be considered as a function of d for 
fixed C. It is an even function of d, and since d<1, it may be written 
in the form 


ME=A+Bd?+0(d!) . 2... 2 1 (49) 


where A, B are independent of d but are functions of C. From the 
discussion of § 6, it follows that B is positive for low values of C, while for 
higher values of C, it becomes negative. In fact, we know from pre- 
ceding considerations that in the limit of small C, B is given by 


i 1 Len fel 
eB aay gir Ye 15 7 
Be ee Oe UC no 


but that this does not hold to a high enough value of C to establish the 
point where B changes its sign and becomes negative. The lowest 
positive zero of B(C’) gives the density above which ferromagnetism sets 
atl =0: 

The effective masses now depend on d, and one expands the ratios 
7(-+)=M/M*(+) as follows : 


T(-+-)=To9+d7,,+d?77,.+0(d5) . re Ret DF bear (51 a) 
7(—)=T9+d7_,+d?7_.+O(d*). Brie es tea) 
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Furthermore, it follows from (47), (48) that 


I I 
sane se +0(@)}, i ee) 
Define 
odd 
)=— ae (27+ 1) 5,(a) HTS, ct bnew Psiay nea 
and 
all 
b(e)=— & (+1) 5,(x). ee ERY) 


Substituting the above expressions into eqns. (33) and (34) and keeping 
in mind that both the integrands and the limits of integration are 
d-dependent and that only quantities which are of cubic and higher degree 
in d may be neglected, it turns out that 

C= — Lc ( . 3 * r; 5 . = is 54 


T4973. 


Writing 7, for 7,, and 7, for 7,5, this gives 


Te toy UT tO Te. oo mw one (OS) 
T_=T)—d7,+ 077, . . . . . . . . (55 b) 
T9=4(74 +7) =T 99+" 70. Jere Hemera WOO C) 


Having established eqns. (54), it suffices to deal with the function U(— ; k). 
The following results are then obtained : 


1/2C 
LU 30)=tf der af B(0)+$@)} 
1/2C 
+af =A dx .xP(x)— iB eae 2 (40) $000) | 
+e vf dhe . 2{D(0)-+-$(@)}+ 35 O2PUC) 
+ 7055 [—B(IC)—AUO)HICH'EC) HICH'CY} |, = (06) 
C x 
U(—; be) 7m| dx. (1-7) . {D(x)+4(x)} 
0 
1 C 2D(x) x(a) 
+a{ — 3 reo] ee sel, ded C0 (a a yagi } 


: ae eee 2 Too (" eee 
+d? 18 TyC7P(C an 9 ToC : dud(x)— 9 a : a6 % p(x) 
; . 0 3 fn © 


2 To0 2 eel : In. 22@(a 
= | dix wP(x)+ 3G ote (a) 


traf de(1— a) [PC + $e) f er: (OL) 
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Now, from eqn. (23 a), it follows that 
16 
(T9947 +-4?7,—1) sarmmsgy iakbrs ) ky_))— U(— ; O)} (58) 
'F(-) 


which permits the 7», 7, and 7, which make the s.d.e.m.a. self- 
consistent to be deduced. The results are : 


6 a4 
re 4 1= 8 (26) re a Be (59) 
aC? al (ae 
ee 
aC? 
tre Set Ger—9+ Foro) ZO"wEO) |}. (60) 


C a0? (OE, C- 
va{ Te wate LT nl 5-F edo | 
Z 


+r] Tt OOOH E-F 


C2 
+ (a0) +400)—-IeHGC)—208'AC) |P 6D 


where «, 8, y, «, 7, u,v are the following functions of C. (In eqns. (59), 
(60), (61), it is understood that the argument of each one is C). 


(0)=- da . (1-7). {D(a)+d(z)} . . . (62a) 
1/20 

pO)=| dx .a{O(x)+o(x)} . . . . . . (626) 

0)=[ de .x(1— 7) Oe ig Ree ne OG Oe, 
lene 

eC) == . 02 

(C) Clee eeu Ea ken a ey 
°C 6 

ntC)=[ de .2(=— 2) ote) . A Sadat 8 (PRS 
°C C ¢ 

MO)=I de (5 — Fb) Be Se ana (O2EF) 
1/2C 

HO)=| AG ED) Rie ee, eee oie a ere (O20) 

A and B of eqn. ae then ae Pa values : 
=A 04° See = O' 75-1) oa mnt any EGEY 


1 Ou 


B= 024 = = {U2 2(0)— U_,(0)} +O {nt 3 ™1+ - 5 (70D . (64) 
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1 
where U_,(0)=—27,»— g T0P{P(3C)—4(3C)}. et OO) 
1 
U_,(0)=27,8+ 3 7,0? B(5C) 


1 
+ 75 C rool —2P(3C) —26(30) +C8'(3C)+C4'(4C)}. (655) 


For a hard sphere interaction, the following are good approximations for 
® and ¢ over the range 0<a<2 


P(x) = 2a°— 0-96282x5+-0-2900x7—0-0319929 . . (66a) 
$(a)=x+ 0-860323—0-14202°-+0-0154527. . . . (666) 
Fig. 3 


B(C) as a function of C, where B is defined by eqn. (49). 


The function B(C) is plotted in fig. 3, and has a zero at C=1-674. Since 
this is still somewhat within the range of validity of the e.m.a., it is 
concluded that at absolute zero, a gas of Fermi—Dirac particles of spin }, 
interacting by a hard sphere potential of radius a, becomes ferromagnetic 
if the density 

(1-674) 


n> ot eet es ee 


§ 8. ENERGY IN THE FERROMAGNETIC REGION OF DENSITY 
Equations (42), (36), and (31) lead to the following expression for the 
energy per particle E in the s.d.e.m.a. 
3 keh’, * hemi U (+3 0) +h U(—3 0) 


Jiieeetee 2 | 
10 kyyptky_y 7 kp +k 
3 ; haa tut) — 1}+ky_) {7(—) : 1}  —_ere (68) 
20 kaye thm) 
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For C fixed at a value slightly greater than 1-674, ME as an even function 
of d has a local maximum at d=0 and local minima at d=+d, where 
dy<1. As OC becomes larger, dy increases towards 1. At a value of C 
above that for which dj=1, the local minima disappear, and d=+1 
become absolute minima. It is only this latter case, for which E has its 
minima at d=+1, which corresponds to the familiar kind of ferro- 
magnetism. The behaviour in the other case is of an intermediate sort, 
in which the energy is lowered as |d| is increased and thus the magneti- 
zation tends to be increased, provided 0 < |d|<d, (ferromagnetic behaviour) 
while for d)<|d|<1, the energy is lowered as |d| is decreased, and thus 
the magnetization tends to be decreased (antiferromagnetic behaviour). 
This will be referred to as ‘ intermediate ferromagnetism ’. 


Table 2 


1-40 0 1-9257 2-0391 2-0391 
1-40 0-2 1-9348 2-0435 2-0193 
1-60 0 4-2379 3°1522 3°1522 
1-60 0-2 4-2439 3:1406 3°1270 
1-68 0 6-3771 4-1519 4-1519 
1-68 0-05 6-3770 4-1524 4-1477 
1-68 0-2 6-3752 4-1320 4-1101 
1-68 0-5 6-3623 3°9990 3:°8975 
1-68 1-0 6-4160 3°5316 

1-75 0 10-106 5:8764 5-8764 
1-75 0-5 9-871 5:5484 5-4002 
1-75 0-75 9-470 5-0614 4-7154 
1-75 1-0 8-949 4-3948 


Table 2 gives the value of the energy per particle E and of the reciprocal 
effective mass ratios r(-++-)=M/M*(+-) and 7(—)=M/M*(—) for various 
values of C and d. It is seen that for C<1-674 (C fixed), E increases as 
|d| increases, while for C' slightly above 1-674, E has a local maximum at 
d—0, reaches minima at certain values d=-+d,, and increases as |d| 
increases from dy to 1. For somewhat larger O, E is a monotonic de- 
creasing function of |d]. 

Since the densities for which ‘intermediate ferromagnetism’ and 
ferromagnetism occur lie near the limit of validity of the s.d.e.m.a., a 
detailed study of this phenomenon must use a more complicated wave 
number dependence for V(+ ; k), V(—; &), and must also exclude 
A states for the intermediate states in the integrals of eqns. (4) 
and (5). 
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The most obvious application of these results to an actual physical 
system is a study of the paramagnetic susceptibility of ?He, and an 
investigation as to whether or not ?He becomes ferromagnetic at a 
sufficiently low temperature. L. Goldstein and M. Goldstein (1950) have 
concluded that at 7'=0, the antiferromagnetic configuration (d=0) is 
more stable than the totally ferromagnetic one (d=-41). Their deriva- 
tion involves the use of a different set of approximations from those of 
this paper, viz. the use of plane waves, neglect of spatial correlations 
between the atoms, and the use of two types of semi-empirical inter- 
atomic potentials with cut-offs at physically plausible values to avoid 
divergencies. Furthermore, they do not investigate ‘intermediate ferro- 
magnetism’ at 7=0. The attractive part of the *He interaction increases 
the minimum density at which d=0 is most stable with respect to its 
value in the hard sphere case (taking the same radius for the hard sphere 
part of the *He interaction as in the pure hard sphere case). 

A gross overestimate for the Curie temperature of a gas of hard sphere 
fermions is given by 

h2 

Ma? 
where B is the function of C defined in eqn. (49). For M=5:0x10-*4 g 
and a=2-6 x 10-8 em, this overestimate becomes 

ieee B AS estes ott te (c0) 
Equations (69) and (70) were derived on the extremely crude assumption 
that even for 7'>0, all the filled wave numbers for ‘ up’ particles lie 
below all the empty ones in each member of the canonical ensemble, 
and similarly for the ‘down’ particles. However, it was, of course, 
assumed that the number of ‘ up’ and ‘down’ particles in a member of 
the ensemble differs, in general, from the numbers at 7’=0. 


fe 2B (69) 


Cries =a ne 
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CORRESPONDENCE 


The Effect of Free Electrons on Lattice Conduction at 
High Temperatures + 


By R. Srratron 
Metropolitan Vickers Electrical Co. Ltd., Manchester, 17 


[Received December 4, 1956] 


Tue thermal conductivity of a semiconductor is equal to the sum of two 
parts which are due to the transport of energy by current carriers (free 
electrons or holes) and thermal lattice vibrations (phonons) respectively. 
The lattice contribution will be controlled by inter-phonon, phonon— 
lattice defect and phonon-electron (or hole) scattering. The last process 
bas been considered by Ziman (1956) for lattice temperatures (7’) below 
the Debye characteristic temperature (6){. We will give a simple exten- 
sion to the theory for the high temperature range (7'>0@). 

Now a simple argument indicates that whereas at low temperatures 
phonon-electron collisions may be important, this will not be the case at 
high temperatures. Thus according to Ziman (eqn. (7)) it follows from 
the conservation of energy that the minimum wave number k, of an 
electron which is capable of scattering a phonon with wave number 
p (emission or absorption) is given by 


ho= [S000 fh. [aoe de meme) 
where m* is the effective mass of the electron and v is the velocity of the 
phonon. This leads to a minimum electron energy 
peaTithg? Tanto 
Oo 2m*~ 8m*v2 
where w=p/v is the angular frequency of the phonon. Now when 
7’ <6, most of the phonons have an energy of about 2k7' so that 


Fed he ae dt m*y 

It follows that when 7'~m*v?/k(~1°K) most of the electrons are 
capable of scattering the phonons. With increasing 7’ only the higher 
energy electrons will come into the picture until finally, when 7’ is greater 


+ Communicated by the Author. 
t The restriction on the temperature is implied in Ziman’s paper (eqn. (19)) 
since the integration over the lattice frequencies in the expansion for the mean 


A ae is extended to an infinite upper limit which is only permissible if 


+im*e®—fie . . 1S. (2) 
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than 6 and the most important phonon energy is k0, the minimum electron 
energy required becomes 


RESUS me we eo hd 


which is of the order of 1 electron volt. 

It will thus be expected that electron—phonon scattering is relatively 
unimportant at high temperatures contrary to the results of ter Haar 
and Neaves (1956) who have however not taken account of the fact that 
a minimum energy H, for scattering exists. (This will be referred to in 
detail after eqn. (10).) 

The detailed calculation for the thermal conductivity can be carried 
out by taking over Ziman’s expression for the phonon relaxation time = 
(eqn. (9) in Ziman’s paper and corrigendum) applied to the case of high 
temperatures, i.e. 

1 mp kT 1 
t(w) = 4rh4p © 4m*v?* 1+ exp [—L/kT + h2w2/8m*v?] 

Here ¢ is the arbitrary Fermi energy of the electron gas, p is the mass 
density of the lattice and C is an electron—lattice interaction constant 
whose magnitude is of the order of a few electron volts. 

Substituting for z7(w) into the expression for the thermal conductivity 
(cf. Klemens (1951), eqn. (24) for the general expression) namely 


ke peal 7.(w) 


ae 5 
Gr") G8; 


(5) 


CS pean twee. () 


if 7’ >6 where j=1, refers to the longitudinal direction of polarization, 
j=2 and 3 to the transverse, gives 


Php 1 re ) ex (£,/kT)—1 7 
1 ee eee ee Ah) Se Reh 
~ 6am™ Seal +exp (— kT (Bo|kT) 

If (CP) /kT >1, then « is a constant in agreement with the result for 
metals given by Wilson (1953, eqn. 9.82.2) in the temperature range 
AST if 

For a semiconductor, provided that (H)—{)/k7T >1, eqn. (7) can be 
reduced to 

kipe* I exp { — SRDS 
(a a Se aE ee nt) 
™ 6rm*? ei ian 


The unknown quantities C; can be eliminated from eqn. (8) by intro- 
ducing the expression for the electron mobility (cf. Wilson (1953), eqns. 


9.36.3 and 9.36.4) 
4 e porhrke? 


leek 9 
= 3 Qam*hT)¥2 * 4n(6n*)28m*C2T (9) 


where 8 is the lattice constant and C?=C+C?+0;°. Now in the 
elementary theory only C, (longitudinal waves) is non-vanishing which 
would make « infinite since transverse waves (C’, and C3) are not scattered. 
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A more exact discussion however (cf. Wilson 1953, p. 279) indicates that 
there will always be some interaction between the transverse waves and 
electrons; it has been suggested (Makinson 1938) that it is better to 
assume that C,~C,~C3. 

From eqns. (8) and (9) 


(10) 


m*\12 pT _ C2 exp{(Hy—O kT} 
ade 10-28 ( —— bat ged 3 a 
ea ee eS eo) ayer 


cals/sec em°K 


where p is in practical units (em2V—! sec-1), m is the free electron mass 
and if C;=C,—C, then XC?/C,?=9. The result by ter Haar and Neaves 
(1956, eqn. (171)) differs from the above by the factor (H)/kT’)~1 exp (H,/kT) 
and applies only to non-degenerate semiconductors. The omitted factor 
is however a very large quantity so that the thermal resistance (1/«) 
due to electron phonon collisions is very small compared to the experi- 
mental values (of the order of several cmsec°c/cals for covalent semi- 
conductors at room temperature) which must therefore be explained in 
terms of the other phonon scattering mechanisms. 
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REVIEWS OF BOOKS 


The Life of Arthur Stanley Eddington. By A. Vinert Dovauas. (London and 
- New York : Thomas Nelson & Sons Ltd., 1956.) [Pp. xiv+208.] Price 25s. 


Tuts admirably written biography gives an account, with a minimum of 
technicalities, of the work of one of the most distinguished astronomers and 
astrophysicists of his age and draws a portrait of an attractive personality. 

Born in 1882, Eddington entered Trinity College, Cambridge, as a Scholar 
in 1902 and became a Fellow of the College in 1907. From 1906 to 1913 he was 
Chief Assistant to the Astronomer Royal at Greenwich, and from 1913 to 1944— 
the year of his death—Plumian Professor of Astronomy at Cambridge. From 
1914 he combined with his Professorship the Directorship of the University 
Observatory. He was elected to a Fellowship of the Royal Society in 1914, 
was knighted in 1930 and was awarded the O.M. in 1938. 

His first original work was on the motion and distribution of stars, and in 
1914 he published his book entitled ‘‘ Stellar Movements and the Structure of 
the Universe ”’. 

Shortly after this, Eddington was ‘““led to make a close scrutiny ”’ of Einstein’s 
‘theory of relativity ; the results of this study were published in a ‘“‘ Report 
on the Relativity Theory of Gravitation’ (1918). This was followed five 
years later by his “ great, now classical, exposition’’, “‘ The Mathematical 
Theory of Relativity ”’. 

His skill as a practical astronomer and his grasp of the principles of relativity 
made inevitable his inclusion among the team of observers sent to Sobral and 
Principe for the solar eclipse of 1919. One of the objects of the expedition 
was to test the Einstein theory of the gravitational deflection of light, the 
Einstein value of the deflection being double that calculated on the classical 
theory. Mr. Douglas tells the story with a keen sense of drama : 

“In Dyson’s study at Greenwich on the last evening before sailing, he and 
Eddington were discussing the amount of deflection calculated on the classical 
Newtonian theory and the Einstein value, about twice as great, which Eddington 
fully expected to obtain. Cottingham who was to accompany Eddington 
asked ‘ What will it mean if we get double the Einstein deflection ?’ to which 
Sir Frank replied ‘ Then Eddington will go mad and you will have to come 
home alone’. Three months later at Principe, when he had made his first 
plate reduction, Eddington turned to his companion with the words, 
‘ Cottingham, you won’t have to go home alone ’.”’ 

From relativity Eddington returned to stellar physics. His book on the 
“Internal Constitution of the Stars” was published in 1926 and, says 
Mr. Douglas, “ immediately became and has remained the standard treatise 
on this new field of astrophysics which he himself had so largely opened up and 

eloped ”’. 

alae the last years of his life Eddington occupied himself with the 
preparation of a work on ‘‘ Fundamental Theory ”’, lett unfinished at his death. 
Here he brought together relativity and quantum theory—and especially 
Dirac’s work in the latter field—and endeavoured to derive from theoretical 
premises numerical values for the natural constants. Of this work, full of 
mathematical and other difficulties, Douglas remarks that while most physicists 
have “dismissed it from serious consideration ’’, others have found in it 
starting points for further researches. “ Unfinished it stands as a challenge 
to future manipulators of the mathematical symbolism of physics . . .”. 

A bibliography lists 13 books written by Eddington and nearly 160 papers 
communicated to learned societies and other publications, C.F.M, 
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Thermal Power from Nuclear Reactors. By A.S. THompson and O. E. RopGErs. 
(Chapman and Hall.) Price $7.25. 


Tue authors of this new book state that their main object is to present an 
analysis of the engineering of nuclear reactors, with special emphasis on the 
thermal aspects of reactor design. This is undoubtedly an aim to be 
recommended, since though there are now many excellent books on nuclear 
reactor science, there has long been a need for an adequate work dealing 
with the special problems of heat transfer and heat removal. 

The first two chapters are presented as an introduction to the required 
reactor physics and reactor theory, and give such background as is necessary 
for an understanding of the principles involved in the nuclear physics design 
of reactors. The interactions of the elementary particles, diffusion theory 
and the calculation of critical mass are discussed, and a particularly interesting 
treatment of the numerical solution of the reactor equations is also presented. 
In the following chapter an excellent analysis is given of reactor kinetics, in 
which the interdependence of reactor power and temperature, through the 
temperature coefficients, is well brought out. The authors show that a satis- 
factory solution to such problems is equally dependent on an adequate 
calculation of component temperature variations and neutron flux changes. 

After rather brief treatment of reactor shielding problems and desirable 
reactor material properties, the remainder of the book is devoted to thermal 
aspects of reactor design. This part opens with a most useful chapter dealing with 
thermal stress problems which arise in reactor design. This is a subject rather 
neglected in previous reactor engineering text books, and it is pleasing to see 
the authors developing working theories and showing how they may be applied 
in various examples. Certain heat conduction solutions are also given in this 
section, as being necessary for calculation of the temperature distributions 
giving rise to thermal stress. In following chapters the calculation of pressure 
losses, the working heat transfer correlations and boiling heat transfer are 
dealt with. It is felt that some discussion on the necessity for considering 
compressible flow with gas cooling and with heat addition could usefully have 
been added, and also that the conditions for application of the heat transfer 
correlations should have been made clearer and emphasised. There is now 
considerable information available on boiling heat transfer and although it 
is a rather specialised topic, more of this might perhaps have been usefully 
included. 

The book concludes with a chapter on thermal power cycles giving the 
fundamentals of both steam and gas turbine cycle analysis. 

Considering the book as a whole, it is felt that it should be a very useful 
introduction to the special problems of reactor engineering, especially to those 
working in other engineering fields who would like a link up with their well 
established practices. D.V.W. 
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First row: Figs. 1, 2, 3, 4. 
Second row: Figs. 5, 6,.7, 8. 
Third row: Figs. 9,10, 11, 12. 
Fourth row: Figs. 13, 14, 15, 16. 


Order of magnification given in the text. 
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‘irst row: Figs. 17, 18, 19, 20. 
Second row: Figs. 21, 22, 23, 24. 
Third row s) Wigss 25,26) 27.028) 
Fourth row: Figs. 29, 30, 31, 32. 

Order of magnification given in the text. 
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Fic. 5 
ig. 5 


Slip-lines on a copper crystal deformed at room temperature by additional 
shear strains 4e—0-05 after various prestrains e. (a) prestrain e—0; 
(6) e=0-1; (c) e=0-2; (d) e=0:3; (e) e=0-6; (f) «e=0-6 (a-e) dark 
field illumination ; (jf) bright field illumination, same site as (e). All 
photographs x 225. 


A. SEEGER et al. Phil, Mag: Sera 8; Vole 206i = 


Slip markings on the top surface of an Slip markings on the top surface of a copper 
aluminium crystal deformed at room crystal deformed at room temp- 
temperature showing a typical slip erature. No clustering of slip-lines 
band. into bands. Shear strain ¢«=0-:3. 
Shear strain «=0-3. Electron micro- Electron micrograph » 7500. 


graph x 7500. 


Fig. 11 Fig. 12 


Slip line pattern of an aluminium crystal Clustering of slip lines into bands on an 
deformed at 90°K in stage II. Pre- aluminium crystal deformed at 90°K 
strained up to a flow-stress 7=0-5 in stage III.  Prestrained e¢==0-3, 
kp/mm®, polished, strained in addition polished, and strained in addition by 
up to a flow-stress 7=1:3. kp/mm? de=0-10. — Electron micrograph x 


corresponding to a strain increment 7500. 
Ae=—0-07. Electron micrograph » 

=e 

7500. 
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Fig. 13 . Fie. 14 


Electron micrographs of fragmented slip-bands on the surface of a copper 
crystal prestrained «=0-6 at room-temperature. Additional strain after 
polishing 4e=0-05, «9000. Same crystal as figs. 5 (¢) and (/). 


Fig. 18 Big, 19 


Aluminium crystal, prestrained at 90°K Same site as fig. 18, x 3465. 
€=0°35, polished, and work-softened 
at room temperature. Bright field 
illumination, x 90. 
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Fig. 20 Fig. 21 


Same crystal as figs. 18 and 19, electron Same crystal as figs. 18, 19, 20, showing 
micrograph, x 5000. traces of the primary slip plane, 
the conjugate slip plane, and the 

cross slip plane, 5000. 


(a) () 
Copper crystal, prestrained «=0-6 at 90°K, polished, and work-softened at 
300°c, x 200. (a) Dark field illumination ; (b) Same site in bright field 
illumination. 


Same crystal as fig. 22, electron micrograph, x 8000. 
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(a) 
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(a) General view of helical dislocations, 
(6) Helical dislocation with closed loops. 
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